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Chapter 1 
Introduction 



1.1 Overview 

The origin of baryon asymmetry (matter-antimatter asymmetry) in the present universe 
is one of the fundamental puzzles in particle physics as well as in cosmology. Assuming 
an inflationary phase in the early universe, any initial baryon asymmetry is diluted and 
becomes essentially zero during the inflation. Therefore, the observed baryon asymmetry 
should be generated dynamically after the inflation. Such a dynamical generation of 
baryon asymmetry (baryogenesis) is possible if the following three conditions are satisfied; 
(i) baryon-number violation, (ii) C- and CP- violations, and (iii) departure from thermal 
equilibrium, which are known as the Sakharov's three conditions |]]. 

The original idea of the baryogenesis in the grand unified theory (GUT) || satisfies 
all of these conditions |§, and provides an elegant interplay between particle physics and 
cosmology. Actually, the GUT [|J unifies baryons and leptons in the same gauge multiplets 
and predicts the existence of baryon number violation, the delayed decay of superheavy 
GUT particle satisfies the out-of-equilibrium condition, and its CP- violating asymmetric 
decay into baryons and antibaryons can offer a net baryon asymmetry. 

On the other hand, the electroweak gauge theory itself violates the baryon asymmetry 
by a quantum anomaly M . This effect is highly suppressed by a large exponential factor at 
zero temperature, and hence the stability of the proton is practically guaranteed. However, 
it was found that in the early universe at high temperatures above the electroweak scale, 
the baryon-number violating interaction is not suppressed and even be rapid enough to be 
in thermal equilibrium [El. An important point here is that this baryon number violating 
( "sphaleron" ) process violates lepton (L) number as well as baryon (B) number, and 
it conserves a linear combination of them, B — L. Therefore, if the baryon asymmetry 
is produced in (B — L)-conserving processes, like those in the SU(5) GUT, it would be 
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washed out by the sphaleron effect before the electroweak phase transition, and hence no 
baryon asymmetry can remain until the stage of big-bang nucleosynthesis to generate the 
light nuclei. 

The existence of the sphaleron effect opens a new possibility of baryogenesis, the gen- 
eration of baryon asymmetry at electroweak phase transition || [7] . This mechanism does 
not use any (B — L)-violation, and the baryon asymmetry is generated by the sphaleron- 
induced (B + L)-violating process. (Hence, same amount of baryons and leptons are 
produced.) As other baryogenesis scenarios have, this mechanism also has a double-edged 
behavior, i.e., the produced baryon (and lepton) asymmetry tends to be washed out after- 
wards by the sphaleron process itself. In order to avoid the erasure of the produced baryon 
asymmetry, the electroweak phase transition must be strongly first-order. However, in the 
case of the minimal standard model with one Higgs doublet, it turns out that, for Higgs 
masses allowed by current experimental bounds, the electroweak transition is too weakly 
first order or just a smooth transition, which excludes the electroweak baryogenesis in the 
standard model. Although the extension of the standard model to the supersymmetric 
version may cure this difficulty, there remains only a small parameter range ||. 

Therefore, it seems natural to consider that the baryon asymmetry was generated 
before the electroweak phase transition. As long as we consider a baryogenesis above the 
electroweak scale, some (B — L)-violating interaction is mandatory, since the sphaleron 
process washes out any baryon asymmetry in the universe unless there exists nonzero 
B — L asymmetry. (Namely, the first one of the Sakharov's conditions, (i) I?-violation, 
should be replaced with (i)' (B — L)-violation.) This means, on the other hand, as first 
suggested by Fukugita and Yanagida |9|], that the lepton-number violation is enough for 
baryogenesis and explicit baryon-number violation is not necessarily required. Actually, 
there now exists an implication of lepton-number (L) violation, whereas there has been 
discovered no evidence for baryon-number (B) violations: that is, the neutrino oscillation. 

Neutrino oscillation [HJ.Q especially the strong evidence for the atmospheric neu- 
trino oscillation reported in 1998 by the Super-Kamiokande Collaboration |T0| , 



is one 



of the greatest discoveries in particle physics after the success of the standard model. 
The experimental data strongly suggest that the neutrinos have small but finite masses, 
impelling the standard model of particle physics to be modified. 

Such tiny masses of light neutrinos are naturally understood if we introduce heavy 
right-handed Majorana neutrinos to the standard model, in terms of so-called seesaw 
mechanism [O]. Because the right-handed neutrinos N are singlets under the standard 



1 Very recently the KamLAND experiment has announced the first results, which exclude all oscillation 
solutions but the 'Large Mixing Angle' solution to the solar neutrino problem 
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gauge symmetries, they can have Majorana masses as well as Yukawa couplings to the 
left-handed leptons (I) and Higgs (<p) doublets:?] 



£ = --MN T N+ (-hN T l V + R.c.) , (1.1) 



2 

where we have taken the N to be the Majorana mass eigenstates. Then by integrating 
out the heavy right-handed neutrinos, we obtain the small Majorana masses for light 
neutrinos, m„ = h 2 (ip) 2 /M. 

A crucial observation in the lagrangian Eq. fll.lj ) is that the Majorana mass term 
MN T N represents nothing but a lepton number (L) violation, actually a (B— L)-violation. 
In fact, the heavy right-handed neutrinos have two distinct decay channels into leptons 
and anti-leptons, and produce B — L asymmetry if the Yukawa couplings h violate the 
CP and if the decay is out of thermal equilibrium. The produced lepton asymmetry is 
partially converted into baryon asymmetry via the aforementioned sphaleron effect, which 
explains the baryon asymmetry in the present universe. This is the original idea of the 
leptogenesis ||. 

We should also note that if we assume a gauged XJ(1)b-l symmetry, the existence 
of three generations of the right-handed neutrinos are automatically required in order to 
cancel the gauge anomaly. As well known, the U(1)#-l symmetry is the unique extra 
U(l) symmetry which can be gauged consistently with the standard model. Furthermore, 
the breaking of this U(1)b-l symmetry naturally provides large Majorana masses to the 
right-handed neutrinos, which leads to the tiny neutrino masses via the seesaw mechanism. 
It should be noticed that the \J(1)b-l symmetry is consistent with the SU(5) GUT, and 
is embedded in larger GUT groups such as SO (10). 



Meantime, the supersymmetry (SUSY) [jrjj has been attracting wide interests in par- 
ticle physics as one of the best candidates for new physics beyond the standard model. 
It protects the huge hierarchy between the electroweak scale and unification (or Planck) 
scale against the quadratically divergent radiative corrections, and the particle contents 
of the minimal SUSY standard model (MSSM) lead to a beautiful unification of the three 
gauge couplings of the standard model at the scale ~ 2 x 10 16 GeV |[15|| , which strongly 
suggests the SUSY GUT [IB]. The MSSM also gives a natural framework to break the 
electroweak symmetry radiatively [IT]]. On the other hand, from the viewpoint of cosmol- 



ogy, SUSY provides an ideal dark matter candidate, the lightest SUSY particle It 
also protects the flatness of the inflaton potential against the radiative corrections, which 
is inevitable for successful inflation. 



2 Here, we have omitted the family- and SU(2)i-indices for simplicity. 
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However, SUSY also causes a cosmological difficulty if the reheating temperature of 



the inflation is too high, that is, the cosmological gravitino problems |]T9|, [s(J. After 
the end of the inflation, the gravitinos are produced by scattering processes of particles 
from the thermal bath, and its abundance is proportional to the reheating temperature. 
Because the gravitino's interaction is suppressed by the gravitational scale, it has a very 
long lifetime unless it is completely stable, and its decay during or after the big-bang 
nucleosynthesis (BBN) epoch (t ~ 1-100 sec) might spoil the success of the BBN ]19|, |2T| . 
On the other hand, if the gravitino is completely stable, its present mass density must 



be below the critical density of the present universe ||20|| . In both cases, the primordial 
gravitino abundance should be low enough, and hence the reheating temperature of the 
inflation is severely constrained from above, depending on the gravitino mass. 

In this thesis, we study in detail several leptogenesis scenarios in the framework of 
the SUSY. There have been proposed, in fact, various leptogenesis scenarios depending 
on the production mechanisms of the right-handed neutrinos: 

• In the simplest and the most conventional leptogenesis mechanism, the right-handed 
neutrino is produced by thermal scatterings [f| |Z|. The delayed decay of the right- 
handed neutrino satisfies the out-of-equilibrium condition a la the original GUT 
baryogenesis. This mechanism requires a relatively high reheating temperature to 
produce the heavy right-handed neutrino, and the aforementioned gravitino problem 
makes it somewhat difficult, depending on the gravitino mass. 

• Another mechanism is given when the right-handed neutrinos are produced non- 
thermally in inflaton decay ||23|| . In this scenario the reheating temperature can be 



lower than the case of thermal production, and the gravitino problem is avoided in 
a wider range of gravitino mass. 

The third mechanism is inherent in the SUSY. The lepton asymmetry is produced 
by the decay of coherent oscillation of the right-handed "s" neutrino, which is the 
super symmetric scalar partner of the right-handed neutrino [f23], |25|| . If the right- 
handed sneutrino's oscillation dominates the energy density of the universe, the 



gravitino problem is drastically ameliorated p6 



In the framework of SUSY, there is yet another, completely different mechanism. The 
lepton asymmetry is produced not by the right-handed neutrino decay, but by a coherent 
oscillation (actually, a rotation) of a flat direction field including the lepton doublet L: 

• The leptogenesis via LH U flat direction [^Sfl , which is based on the Affleck-Dine 
mechanism P?| . 
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1.2 Outline of this thesis 
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The outline of this thesis is as follows. The rest of this chapter is devoted to some reviews. 
We briefly mention the observed baryon asymmetry in Sec. [L3|. A review of the sphaleron 



effect, which is a crucial ingredient of the leptogenesis, is given in Sec. |i.4[ We also briefly 
review the results of the cosmological gravitino problems in Sec. |1.5| . 

In Chapter 0, we discuss leptogenesis scenarios by the decay of the right-handed 
(s)neutrino. First, we study the asymmetric decay of the right-handed neutrino into 
leptons and anti-leptons in Sec. |2.1| . The conventional leptogenesis mechanism where the 
right-handed neutrinos are thermally produced is briefly discussed in Sec. |273. Then we 



perform a comprehensive study of the leptogenesis mechanism in inflaton decay in Sec. 2.3, 



adopting various SUSY inflation models. In Sec. |2]4], we investigate the leptogenesis from 
coherent right-handed sneutrino. 

In the latter half of this thesis, in Chapter |3], we perform a detailed analysis on the 
leptogenesis via LH U flat direction. This scenario may require another overview, which 
will be given in Sec. [3.1| . Here we mention one point, that the most important parameter in 
this scenario which determines the baryon asymmetry is the mass of the lightest neutrino, 
m v \. (Notice that the data from neutrino-oscillation experiments |Lp, [11], [12|] suggest the 
difference of the neutrino mass squared, indicating the masses of the heavier neutrinos, 
m U 2 and m„3.) It is amazing that the observed baryon asymmetry, which was generated 
in the very early universe, is directly related to such a low-energy physics, the neutrino 
mass. 

Chapter £| is devoted to conclusions and discussion. Some of our notations are sum- 
marized in Appendix |A|. Some ingredients of the standard cosmology are reviewed very 
briefly in Appendix. [B]. 
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1.3 Baryon asymmetry 



We will use the following value for the baryon asymmetry in the present universe p8| , |2~9fl : 



Ub _ r^ni = (04 _ 1)xlo -io (L2) 
s s 

where n& and ni are baryon and anti-baryon number density, respectively, and s is the 
entropy density of the universe. This ratio takes a constant value, as long as the baryon 
number is conserved and no entropy production takes place. Notice that ni <C n b and 
hence n# ~ n b in the present universe. (We sometimes call the % just "baryon number 
density," for simplicity.) This value is determined from the big-bang nucleosynthesis 
(BBN). The BBN occurs at the temperature of T ~ (1-0.1) MeV, or equivalently at 
the cosmic time t ~ (1-100) sec, and generates light nuclei, D, 3 He, 4 He, and 7 Li. (For 



reviews, see, for example, Refs. p8| , |29j.) The agreement between the predictions of the 
BBN theory for the abundances of these light nuclei and the primordial abundances of 
them (which are inferred from observational data) is one of the most important successes of 
the standard big-bang cosmology. Actually, the theory of the BBN has basically only one 
free parameter ,0 i.e., the baryon asymmetry, and all the above light-element abundances 



are well explained with ns/s in the range of Eq. ( |1.2|) . 

To demonstrate the smallness of the baryon asymmetry, let us calculate the baryon- 
number and anti-baryon number density at high temperature. For example, consider a 
temperature between the QCD phase transition and electroweak phase transition, say, 
T ~ 10 GeV. At this temperature, all the baryons and anti-baryons are expected to exist 
as quarks and anti-quarks, and they are well in thermal equilibrium. The ratio of the 
baryon-number and anti-baryon-number densities to the entropy density are then given 
by 

^~^~0.02, (1.3) 



where we have used Eq. ( p,19| ). Compared with Eq. ( |1.2|) , this means, (n b —ni) / (n b +ni) 



e>(10" 9 ), i.e., there were (1,000,000,000+1) baryons per 1,000,000,000 anti-baryons at this 
epoch. 



3 We assume the number of neutrino generations to be three. As well known, this fact itself is one of 
the most important implications of the BBN. 
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1.4 Sphaleron 



In this section, we give a brief review of the "sphaleron" effect, which is a crucial ingredient 
of the leptogenesis. If the lepton asymmetry is successfully generated, it is partially 
converted into the baryon asymmetry at equilibrium thanks to the sphaleron- induced 
baryon and lepton number violating process, which then explains the observed baryon 
asymmetry in the present universe. 

At the classical level, the Lagrangian of the standard SU(3)c x SU(2)l x U(l)y gauge 
theory clearly conserves the baryon (B) and lepton (L) numbers. However, as was first 
pointed out by 't Hooft ||, both B and L are violated by quantum effects: 

« = <V£ = Nf^e^TrF^ , (1.4) 

where Nf is the number of fermion generations and #2 arid are the coupling constant 
and the field strength of the SU(2)l gauge group, respectively. It was shown later that 
the SU(2)x x U(l)y gauge theory has a non-contractible loop (path) in the field con- 
figuration space, which connects topologically distinct vacua with different baryon and 
lepton numbers [0O|1 , and it was found that the highest energy configuration along this 
"path" corresponds to a spatially localized and static, but unstable solution pl|, |32fl. This 
solution was named "sphaleron". Since this sphaleron solution is a saddle point of the 
field potential energy, its energy represents the height of the barrier between the vacua 
with different baryon and lepton numbers (see Fig. |1 . 1|) . 

At zero temperature, the rate of baryon (and lepton) number violating process via 
a tunneling between the topologically distinct vacua (for example, proton decay) is ex- 
tremely small f5|, since it is suppressed by the factor of exp(— 167r 2 /g|) ~ 1CT 170 . However, 
as suggested by Kuzmin, Rubakov, and Shaposhnikov || Q , such processes are not sup- 
pressed and can even be efficient at temperatures close to (and above) the electroweak 
phase transition. The transition of the fermion number occurs simultaneously for each 
fermion doublets, so that the total changes in the numbers of baryons and leptons are 
AB = AL = Nf. Therefore, B — L is conserved. (This is clear from the absence of the 
anomaly in B — L in Eq. ( |1.4| ).) 

Let us roughly estimate the temperature T* above which this (B + L)- violating process 
becomes in thermal equilibrium, according to Ref. |J. We assume that the electroweak 
phase transition is the second order, and consider the temperature T < Tc, where Tc is 
the critical temperature above which the vacuum expectation value v (T) of the Higgs field 
vanishes. In such a situation, the transition from one vacuum (say, B = 6 an d L = l ) 
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E 




Figure 1.1: A Schematic behavior of the energy dependence on the configuration of 
the gauge and Higgs fields [A(x),(p(x)] M. The minima correspond to topologically 
distinct vacua with different baryon (B) and lepton (L) numbers. The configuration 
[A sph (x), (p S ph(x)} represents the saddle point of the energy functional, the sphaleron 
solution. 



to the next vacuum (B = bo ± Nf and L = Iq ± Nf) occurs at the rate 



T = C{T) Texp 



M sph (T) ' 
T 



;i.5) 



where dimensionless factor C{T) depends on the ratio v(T)/T and the coupling constants.^ 
M sv \i{T) represents the free energy of the sphaleron configuration (at temperature T), 



which is given by |3T[ 



M sph (T) = 4tt5(T) 



v(T) 
92(T) 



1.6) 



where B(T) depends on the gauge coupling g2(T) and the 4-point coupling constant of 
the Higgs potential A(T) as B = B(X/ varying from 1.5 (\/g% — > 0) to 2.7 (A/p| — > 
oo) [PHI . The rate in Eq. ( |T~5| ) should be compared with the Hubble expansion rate 
# = (n 2 gj90)^ 2 x T 2 /M G . (M G = 2.4 x 10 18 GeV is the reduced Planck scale and g m is 
defined in Appendix [B.2| .) Then, it is found that the sphaleron rate in Eq. ( |1.5| ) indeed 
exceeds the Hubble expansion rate for T > T*, where is given by 



^2(^) 



hi 



/M, 



V T 



-1 



;i.7) 



See comments below. 
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Here, we have neglected the In C(T) and ln(7r 2 g*/90) compared with the large number 
ln(M G /T*). Recalling that v(T) approaches to zero for T — > T c , this occurs below the crit- 
ical temperature, i.e., T* < Tc- Therefore, the sphaleron-induced (_B+L)-violating process 
becomes in equilibrium for T* < T(< Tc). Precisely speaking, a more careful treatment 
including the prefactor C(T), which turns out to be proportional to (v(T)/T) p4 



is 



required. A numerical estimation shows that T* is just below the critical temperature, 
(T c - T m ) < 10 GeV, for 50 GeV < T c < 200 GeV p|. 



The evaluation of the "sphaleron" rate for T > Tc is a complicated problem, since 
the sphaleron configuration no longer exists in the symmetric phase. | Naively thinking, 
there seems no reason for the (B + L)-violating processes to be suppressed, since the po- 
tential barrier vanishes. Actually, theoretical arguments as well as numerical calculations 
suggest [36| that the sphaleron rate per unit time per unit volume is T/V ~ ct|T 4 for 
T > Tc, where a<i = p|/(47r). By using their result, it is found that the sphaleron rate 
exceeds the Hubble expansion rate for T<10 12 GeV. 

Therefore, the sphaleron-induced (B + L)-violating process is in thermal equilibrium 
in the range of 



100 GeV-T* <T<10 12 GeV. (U 



• relation between baryon and lepton asymmetry 

Let us calculate the relation between the baryon and lepton number in the presence of 
the sphaleron process, by means of the analysis of the chemical potentials |37|, [38|. At 
first sight, it seems that the relation is just given by B + L = since the sphaleron effect 
violates B + L, preserving B — L. However, as we will see, a nontrivial relation is derived 
if there exists a non-vanishing B — L asymmetry. 

We first consider the standard model without supersymmetric particles, but with 
Higgs doublets, <fi (i = 1 • • • N v ), and Nf generations of fermions, i,e., left-handed quark 
doublets qj, right-handed up- type and down- type quarks Uj and dj, left-handed lepton 
doublets lj and right-handed charged leptons Ej (j — 1 • • • Nf). Hereafter, we consider the 
symmetric phase T > T c . Hence, the SU(2) L x U(l)y symmetry recovers. We denote the 
asymmetry of the number density of particle i by Arij, which is related to the chemical 



5 Although the saddle-point sphaleron solution does not exist in the symmetric phase, we will call this 
anomalous (B + L)-violating process "sphaleron" also for T > Tc- 
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potential of that particle as follows [see Eq. ([B.20|) 1 



Arii = — n,- "* 



-giT 3 ) for fermion . 



6" VT 



-giT 3 ( ^ ) for boson . 



(1.9) 



3" \T . 

Therefore, the baryon and lepton number asymmetries are given by 

= An/. - An^ 

= ir 2 x (2^-^) , (1.10) 

where we have used the fact that the chemical potentials of particles in the same gauge 
multiplet are the same, which is ensured by the gauge interactions. (Note that the SU(2) L 
is also recovered, since we consider T > T c .) 

In the following, we derive the relations between the chemical potentials /!». First of 
all, we assume that all the Higgs doublets have the same chemical potentials due to the 
mixings between themselves, i.e., = (If there is a Higgs doublet with a conjugate 
quantum number, like in the MSSM, we redefine the chemical potential of that Higgs 
doublet with an additional minus sign.) Then, the interactions via Yukawa couplings 



y e:j ip* lj ej , y u ,jk^qjU k} y dtjk ip* qjd k , (1.11) 



lead to 



- AV + + Hej = j = 1 • • • Nf , 

AV + H + Vu k = j, k = 1 ■ • • Nf , 
-^ + ii qj +ii dh = j,k = l---N f . (1.12) 

Here, we have taken a basis of gauge eigenstates for the quarks. Thus, because of the 
mixing in the Yukawa couplings, the chemical potentials of the quarks become generation 
independent: fA q . = fi q , = and /i^ = /ij (j = 1 • • -Nf). Notice that the relations 
in Eq. ( |1.12|) hold as long as these interactions are in thermal equilibrium, and even the 
electron Yukawa coupling, which is the smallest one, is in equilibrium for T ^ 10 4 GeV. 
Next, the charge neutrality of the universe requires vanishing total U(l)y charge: 

E(An i y i ) = 0, (1.13) 
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where denotes the U(l)y charge of particle i. This leads to 



;i.i4) 



Finally, the sphaleron interaction can be understood as 

< — > vacuum . 

i = all SU(2) L 
doublet fermions 



which leads to 



3N f fi g + = • 



(1.15) 



;i.i6) 



From Eqs. ( |1.12| ), ( |1.14 ) and ( 1.16j) , all the chemical potentials can be written in terms 
of Nf independent ones, say, . Then, from Eq. ( |1.10| ), we obtain 



n B 
n Lj 



2N f + 3N L 



f j 



2N f + 3N 



J2 Pei. - 3 ^ 



i 

6 

-T 2 x 
6 

V y -UN f -9N 

6 2N f + 3N v y^ e] ' 



f k 



1.17) 



Because the conserved quantity is the B — L asymmetry in the present analysis,^ it is 
suitable to rewrite the above relations in terms of n B — n^: 



n B = C (n B - n L ) , 



where C is given by [jT?H 



8N f + AN V 
22N } + IZN V 



1.18) 



;i.i9) 



In the case of the standard model, Nf = 3 and iV^ = 1 yield C = 28/79. 

We can see from Eq. ( [1.18 ) that, if B — L asymmetry is absent, any baryon asymmetry 
vanishes at equilibrium in the presence of the sphaleron effect. (See the arrow (i) in 
Fig. |L2| .) On the other hand, if there occurs a successful leptogenesis, i,e., if nonzero 



6 Actually, not only the total B — L but each (l/Nf)B — Lj (j = 1 • • • Nf) is conserved in the present 
situation. 
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B 




Figure 1.2: The relation between baryon (B) and lepton (L) number. Thin dotted lines 
correspond to constant B — L, along which B and L can move via the sphaleron process. 
At equilibrium, B and L reach the thick dashed line, which represents B = C(B — L). 



lepton asymmetry (and hence B — L asymmetry) is generated in an out-of-equilibrium 
way, it is partially converted into baryon asymmetry |§. (See the arrow (ii) in Fig. |1 . 2| . ) 
The amount of the baryon asymmetry at equilibrium is obtained from Eq. ( |1.18| ) as 

n B 
s 



n B - n L 



eq 



_ c n. 



eq 



(1.20) 

initial 



where we have normalized the number densities by the entropy density, so that they 
become constant against the expansion of the universe. 

In the case of the MSSM, we have an additional Higgs doublet as well as supersym- 
metric partners. Let us calculate the coefficient C in the presence of those particles. 
First, the (gaugino)-(fermion)-(sfermion)* interactions ensure that the chemical poten- 
tial of each sfermionQ is the same as that of the corresponding fermion. (Note that the 
chemical potentials of the gauginos vanish since they are Majorana particles.) Next, the 
super symmetric mass term W = jj,H u H d makes the chemical potentials of the up-type and 



7" 



sfermion" denotes a scalar partner of the fermion. 
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down-type Higgsinos have opposite signs, and hence the relation induced by the sphaleron 
effect, Eq. (|1.16|) , does not change. However, the condition of the vanishing U(l)y charge 
[Eq. Q1.14Q 1 changes. (Note that the contributions from bosons and fermions are different. 
See Eq. (|1.9|) .) Consequently, the coefficient C becomes again C = 28/79, in spite of 
the presence of two Higgs doublets. In the actual thermal history, however, the Higgsino 
and sfermions likely become massive and are essentially absent at the time of electroweak 
phase transition. Hence, it is more appropriate to calculate C without them, but with 
two Higgs doublets. This (N f = 3, N v = 2) leads to C = 8/23. 

So far, we have considered the symmetric phase T > Tc- The actual ratio of the 
baryon to lepton asymmetry in the present universe depends on how the electroweak 
phase transition occurs f3~8fl . More precisely, it depends on the value of v(T)/T just 
before the sphaleron decoupling T ~ T*. Fortunately, however, numerically it does not 
change much from Eq. ( 1.19 ). (The coefficient C changes from 8/23 {y <C T) to 10/31 
(v ^> T) for Nj = 3 and = 2 fl3~8"f , so that the difference is at most a few percent.) 
Thus, we will take C = 8/23 ~ 0.35 throughout this thesis, for simplicity. 

Finally, we should also note the sign of the baryon asymmetry. We know that the sign 
of the present baryon asymmetry is positive, i.e., hb/s > 0.0 Thus, from Eq. (|1.20|) , it 
is found that the leptogenesis must generate a lepton asymmetry with a negative sign, 
n L /s < (i.e., more anti-leptons than leptons should be produced). In principle, the 
sign of the generated lepton asymmetry could be determined if we know the sign of the 
effective CP-violating phase in each leptogenesis mechanism. In the case of leptogenesis 
by the decay of right-handed neutrino (discussed in Chapter 0), it depends on the phases 
of the Yukawa couplings of the right-handed neutrino [see Eq. ( |2.14j )1. On the other hand, 
in the case of leptogenesis via LH U flat direction (discussed in Chapter |3|), the effective 
CP- violating phase is determined by the relative phase between the phase of the SUSY- 
breaking term and the initial phase of the flat direction field 0, which depends on the 
coupling of cf) to the inflaton field. 

In both of those cases, however, it is highly difficult to determine the sign of the 
phases. Thus, we will simply assume that the leptogenesis mechanisms we will discuss 
produce the lepton asymmetry with a correct sign. Keeping in mind the discussion above, 
we will omit the relative sign in Eq. ( |1.2U| ) for simplicity, and use the following relation 
throughout this thesis: 

(1.21) 



Tin fir 

-Z- = 0.35 x — 



initial 



8 This is not a matter of definition, since we can distinguish the matter from antimatter by CP- 
violating processes in the laboratory experiments (e.g., the asymmetry in the decay — > 7r ± e =F ^), which 
is independent of the definition of "matter" we would name from the cosmological baryon asymmetry. 
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1.5 Cosmological gravitino problems 



In this section, we briefly review the results of the cosmological gravitino problems ob- 
tained in the literature, since they give a very important and severe constraint on the 
baryogenesis, i.e., the upper bound on the reheating temperature. For a review, see 
Ref. ra. 



• unstable gravitino 

There are two cases; unstable and stable gravitino. (See Fig. [L3|.) Let us first consider the 
unstable gravitino. Since the couplings of the gravitino with ordinary matter are strongly 
suppressed by the gravitational scale, it has a very long lifetime: 

r 3/2 ^4x 10 5 x (^^) 3 sec, (1.22) 

where m 3 / 2 denotes the gravitino mass, and we have assumed that the gravitino dom- 
inantly decays into a photon and a photinc] 9 ] and omitted the phase space suppression 
of the decay rate, for simplicity. Therefore, it decays after the big-bang nucleosynthesis 
(BBN) epoch (t ~ 1-100 sec), unless gravitino is heavier than ~ 10 TeV |fl~0"| |. Then the 



energetic photon (or some charged particle) produced in the gravitino decay induce elec- 
tromagnetic cascade process. This cascade might destroy the light elements and change 
their abundances, and spoil the success of the BBN |1| . Since the abundance of the grav- 
itinos produced from thermal environment at reheating epoch is roughly proportional to 
the reheating temperature T R , usually there are upper bounds on the T R depending on 
the gravitino mass. Recently, a detailed calculation of the gravitino production rate in 
SUSY QCD at high temperature has been done [fflUf^l A recent analysis of the effect of 
the radiative decays of massive particles on the BBN is found in Ref. W^ - By using their 
results, the upper bounds are given by Tr ^ 10 6 , 10 9 , and 10 12 GeV for m 3 / 2 = 100 GeV, 
1 TeV and 3 TeV, respectively^ 

If the gravitino is even heavier, as in the anomaly mediated SUSY breaking models |43[] , 
it decays during or near the BBN epoch (t ~ 10~ 2 -10 2 sec). In this case, there is another 
constraint which comes from hadronic decay ||2T|. If energetic hadrons are emitted at 
this epoch, they interconvert the neutrons (n) and protons (p) in the background even 



9 If the gravitino mainly decays into a neutrino and a sneutrino, the upper bound on the reheating 
temperature becomes higher. See discussion in remarks below. 

10 Their calculation |ll| shows a slightly smaller abundance of the produced gravitinos than in earlier 
works. 

. ' , \ \ \ ' i • 1 I . * \ - / • \ - \ 1 - . ■ \ 1 \ \ 1 I j- \ / I 1 1 1111 M f \ I • > <. ■ ^ ■ { : , t \ i \ ' ' f '■ 



11 Here, we have taken the gluino mass to be nig = 1 TeV. Precisely speaking, the abundance of the 



gravitinos produced by thermal scatterings depends on the gluino mass as oc {1 + m?(T) / (3my 2 )} [|_[ 
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Figure 1.3: A schematic figure of the thermal history of the universe relevant to the 
cosmological gravitino problems. 



after the freeze-out time of the n/p ratio (t ~ 1 sec), which results in the change of 
the abundances of the light elements. Thus, there are upper bounds on the reheating 
temperature Tr times the branching ratio of the gravitino decay into hadrons B^. By 
using the results of Ref. [O] and a recent analysis of the effects of the hadronic decay 



on the BBN in Ref. |0J, the upper bounds on the reheating temperatures are given by 
(Sfc/0.1) x T R < 10 9 -10 n GeV for m 3/2 ~ (a few - 100) TeV.Q 



• stable gravitino 

Now let us turn to discuss the case of stable gravitino, which is the case if the gravitino 
is the lightest SUSY particle (LSP). (We assume the .R-parity and hence the LSP is 
stable. A gravitino much lighter than the weak scale is realized in a low-energy SUSY 
breaking scenario such as gauge- mediated SUSY breaking models j45|{3) in this case, 

12 We have taken the gluino mass to be m§ = 1 TeV. See footnote [lT[ 

13 This is not the case if the SUSY breaking is mediated by a bulk gauge field in higher dimension 
spacetime [[h|. 
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the upper bound on the reheating temperature Tr comes from the requirement that the 
energy density of the gravitino not overcloseQ the present universe p0 |. From the result in 



Ref . |^T[ , the relic abundance of the gravitinos which are produced by scattering processes 



of particles from the thermal bath after the inflation is given byQ 



Here, rrig is the gluino mass, h is the present Hubble parameter in units of 100 km 
sec -1 Mpc -1 and Q 3 / 2 = Pz/2/Pc ■ (P3/2 and p c are the present energy density of the 
gravitino and the critical energy density of the present universe, respectively.) It is found 
from Eq. (1.23) that the overclosure limit fi 3 / 2 < 1 puts a severe upper bound on the 



reheating temperature Tr, depending on the gravitino mass 1713/2- Here, we have omitted 
the contribution from the decays of squarks and sleptons into gravitinos. This effect makes 
the upper bound on the reheating temperature slightly severer for a smaller gravitino mass 
region, m 3 / 2 ~ (1-100) keV PD |. For a even lighter gravitino m 3 / 2 ~ 1 keV, there is no 



cosmological gravitino problem [i7] , since in this case the gravitino does not overdose the 
energy density of the universe even if it is thermalized. 

Meanwhile, we should also take care of the decays of the next-to-lightest SUSY particle 
(NLSP) PD[ . Since the NLSP can decay only to the LSP gravitino through the suppressed 
interaction, it has a long lifetime. Then its decay during or after the BBN might destroy 
the success of the BBN, just like the case of unstable gravitino. Thus, there are constraints 
on the abundance of the NLSP at the decay time depending on the lifetime of the NLSP. 
(The lifetime of the NLSP depends on the NLSP mass as well as the gravitino mass, and 
its abundance at the time of the decay is determined by its annihilation cross section.) 
This constraint leads to an upper bound on the gravitino mass m 3 / 2 [p0[ , depending 
on the mass and couplings of the NLSP. Then combined with the constraint from the 
overclosure limit of the gravitino explained above, we obtain an upper bound on the 
reheating temperature Tr. Detailed analyses show that reheating temperature can be as 



high as T R ~ 10 9 -10 1U GeV for m 3/2 ~ 10-100 GeV [it . 



14 The word "overdose" might not be appropriate, since the open or flat universe (fl < 1) cannot 
change to a closed (57 > 1) universe. The actual problem is that the Hubble expansion would be too high 
compared with the observation at a temperature of T = To = 2.73K, if the calculated value of f2 exceeds 
unity ||. 

15 Besides the factor in Eq. ( |l.23| ), the abundance of the gravitino depends on the reheating temper- 
ature also through the SU(3)c running coupling, 0:3 (Tr). Precisely speaking, the final abundance is 
proportional to ^(T^) 3 x {a logarithmic correction (known)} for rrig 3> TO3/2 EjJ]. 
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• remarks 



Several comments are in order. First, all of the above arguments have assumed that no 
extra entropy production takes place after gravitinos are produced. If there is a dilution 
of the gravitino by a late-time entropy production, the bounds discussed in this section 
are relaxed. Care has to be taken in this case, however, since baryon asymmetry is also 
diluted if the baryogenesis occurs before that entropy production. 

Next, in the case of unstable gravitino, if the gravitino decays mainly into a neutrino 
and a sneutrino, the upper bound on the reheating temperature becomes weaker since the 
neutrino have only weak interactions. In this case, the high energy neutrinos emitted from 
the gravitino decay scatter off the background neutrinos and produce charged leptons, 
which cause electroweak cascade and produce many photons. The requirement that those 
photons do not alter the abundances of light elements gives an upper bound on the 
reheating temperature, T R < 10 10 -10 12 GeV for m 3/2 ~ 100 GeV-1 TeV |9| . 



In some cases, considerably higher reheating temperatures are allowed. For example, if 
the LSP is the axino (which is a fermionic superpartner of the axion) and the gravitino is 
the NLSP, the reheating temperature can be as high as 10 15 GeV for 7713/2 — 100 GeV |50 |. 



Another interesting case is given when the leptogenesis takes place from the universe 
dominated by the coherent oscillation of the right-handed sneutrino pH] (see Sec. [2.4|). 



See also a model in Sec. |3.3.2| . In the context of gauge-mediated SUSY breaking models, 
an attractive scenario to solve the gravitino problem has been proposed ||51|| , recently. 

Finally, we mention the nonthermal production of the gravitino during the preheating 



epoch after the inflation [52, [53|]. At the preheating epoch, gravitinos can be produced 



either through the scattering of particles which are created by the parametric resonance of 
the oscillating inflaton [52] or directly from the oscillating inflaton [53|. These nonthermal 
productions may increase the primordial abundance of the gravitino and hence might make 
the upper bound on the reheating temperature Tr severer than discussed in this section. 
However, in both cases the produced gravitino abundance depends on the amplitude 
and coupling of the oscillating inflaton. In particular, as for the second mechanism, 
it was shown that the gravitino abundance produced directly from oscillating inflaton 
is sufficiently small as long as the two sectors, the one responsible for supersymmetry 
breaking at true vacuum and the one for the inflation, are distinct and coupled only 
gravitationally ||54j| . (Notice that this is the case for the SUSY inflation models which 
will be discussed in Sec. |2.3| .) We will use, therefore, the constraint from the thermally 
produced gravitinos conservative bound. 
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Chapter 2 

Leptogenesis by the decay of 
right-handed neutrino 

In this chapter, we discuss leptogenesis scenarios by the decays of right-handed neutrino. 
The lepton asymmetry is produced by the CP- violating decay of the right-handed neutrino 
N into leptons L and anti-leptons L. As discussed in Sec. [L.4j , the produced lepton 
asymmetry is partially converted to the baryon asymmetry |J by the sphaleron process || . 

In Sec. |2.1| we discuss the amount of the lepton asymmetry produced in the decays of 
right-handed (s)neutrino. Then we turn to discuss each scenario depending on the produc- 
tion mechanism of the right-handed (s)neutrinos. The original, and the most extensively 
studied mechanism is the thermal production. We briefly discuss it in Sec. |2.2| . Next, 
we investigate the production of right-handed (s) neutrinos in inflaton decay in Sec. |!0| , 
adopting several SUSY inflation models. Here, we also study in detail the inflation dy- 
namics in each inflation model. In Sec. [2.4| , we discuss the leptogenesis from coherent 
oscillation of the right-handed sneutrino. In particular, we mainly discuss the most inter- 
esting case, the leptogenesis from the universe dominated by the coherent oscillation of 
the right-handed sneutrino. 
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2.1 Asymmetric decay of the right-handed neutrino 



Let us start by introducing three generations of heavy right-handed neutrinos to the 
minimal supersymmetric standard model (MSSM), which have a superpotential; 

W = ^MiNiNi + h ia NiL a H u , (2.1) 

where iVj (i = 1,2,3), L a (a = e,fi,r) and H u denote the supermultiplets of the heavy 
right-handed neutrinos, lepton doublets and the Higgs doublet which couples to up-type 
quarks, respectively. (Here and hereafter, we omit the SU(2)i indices for simplicity.) Mj 
are the masses of the right-handed neutrinos. Here, we have taken a basis where the mass 
matrix for Ni is diagonal and real. 

As can be seen in Eq. (|2.1| ), the masses of the right-handed neutrinos violate the 
lepton number. This gives rise to the following two important consequences. First, 
the tiny neutrino masses, which are now strongly suggested by the neutrino-oscillation 
experiments, are explained via the seesaw mechanism |13| . From Eq. (|2.1| ) we obtain the 
mass matrix for the light neutrinos by integrating out the heavy right-handed neutrinos: 

{m v ) a p = hiahip-j^— ■ (2.2) 

i 1 

The second one is the production of the lepton asymmetry by the right-handed neu- 
trino decay || which we discuss now. Because of their Majorana masses, the right-handed 
neutrinos have two distinct decay channels into leptons and anti-leptons. At tree level, 
these two kinds of decay channels have the same decay widths: 

]r Ni = T(N i ^L + h u ) = V (Ni -> L* + hu) 



4 



T ( Ni ^l + Hu ) = T(Ni^J + H*) = ^-T a \h ia \ 2 M t , 

l07T 



\v^, =T{N l ^L + H u ) =T{N l ^l + h u ) = -!- 

Z 1 07T 



(2.3) 



^r~* = r (N* - 1 + h u ) = r (n* — > l* + h*) = ^-j2 a \h ia \ 2 M i} 

where Ni = v c Ri + v c n 1 \ {v c R , /, h u } and {N, L, H u } denote fermionic and scalar compo- 
nents of corresponding supermultiplets {iV, L, H u }, and y and Y* represent antiparticles 
of fermion y and scalar Y, respectively. Here, we have summed the final states over flavor 
(a = e, /i, r) and SU(2)i indices. We can symbolically write the above widths as 

l -T Ni = T(iV J - L + H u ) = T{Ni L + TT U ) = i- J] \h ia \ 2 Mi , (2.4) 
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where N, L and H u (L and H u ) denote fermionic or scalar components of corresponding 
supermultiplets (and their anti-particles). 

If CP is violated in the Yukawa matrix h ia , the interference between decay amplitudes 
of tree and one- loop diagrams results in lepton-number violation || . Hereafter, we con- 
centrate on the decay of lightest right-handed (s)neutrino Ni, since we will consider only 
the Ni decay in the following sections.0 The lepton asymmetry produced in the Ni decay 
is represented by the following parameter e\\ 



ei 



r(JVi -> L + H u ) - Y{N^ ^L + H u ) 



(2.5) 



which means the lepton number asymmetry produced per one right-handed neutrino de- 



cay. Summing up the one-loop vertex and self-energy corrections [fjS], the e\ has the 
following form: 



1 1 v j 

61 ~ ^Wui 



{(Wit) }S 



Ml 



Ml 



(2.6) 



where f v (x) and f s (x) represent the contributions from vertex and self-energy corrections, 
respectively. In the case of the non-supersymmetric standard model with right-handed 
neutrinos, they are given by 



V 

non-SUSY 



X 



I (.r + l)ln(l + ~ 

x 



non-SUSY 1 



x) 



X 



x — 1 



(2.7) 



while in the case of MSSM plus right-handed (s)neutrinos, they are given by 



a 



fx 



SUSY 1 



X 



In 1 + - 

x 



fi 



fuSY 1 



X 



x — 1 
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Hereafter, we use Eq. ( |2.8|) since we assume SUSY. Notice that all of the N\, N\ and N\ in 
Eq. ( |2.3| ) produce the same amount of lepton asymmetry with the same sign [0, |56], 
given by the above e\. Assuming a mass hierarchy M\ <C M 2 ,M 3 in the right-handed 
neutrino sector (i.e., x 3> 1), the above formula is simplified to the following one: 



=2,3 



Mi 



(2.9) 



1 We will assume Mi <C M2, M3. In this case, even if the heavier right-handed neutrinos iv2,3 produce 
lepton asymmetry, it is usually erased before the decays of Ni . 
2 See also Refs. 

3 Here, we have neglected the three body decay of the right-handed sneutrino, e.g., Ni — > L* + Q + u, 
which gives only tiny corrections [Q. We also neglect the effects of the soft-SUSY breaking terms, since 
the mass scale of soft terms (~ 100 GeV-1 TeV) are much smaller than right-handed neutrino mass. 
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Now let us rewrite this e± parameter in terms of the light neutrino mass m Vl so that 
we can relate the lepton asymmetry (and hence the baryon asymmetry in the present 
universe) to the neutrino mass, which is observed by the neutrino-oscillation experiments. 
First of all, by using Im[{(htf) n } 2 } =0, we can write the 6\ as follows: 



Mi 



8tt (htf) 



-Im 



11 



hh^—h*h T 
M 



11 



(2.10) 



where a matrix notation M = diag(M 1; M 2 , M 3 ) is adopted. Next, by using the seesaw 
formula in Eq. ( |2.2| ), it is reduced to 



3 M l lra\h(m* u )h T 



fi 



ii 



87T (H u f (htf) 



(2.11) 



ii 



The neutrino mass matrix (m u ) a/ 3 can be diagonalized by an unitary matrix U a i (Maki- 
Nakagawa-Sakata matrix |59| ) as (m u ) a p = J2i U^UpiiwZ)^ where is a diagonal mass 
matrix = diag(mj,i, m„ 2 , m u3 ) and m v \ < m u2 < m v3 . Then, with rotated Yukawa 
couplings 



(2.12) 



we obtain 



ei 



3 M 1 Im 


h(m v )h T 


11 


Svr (H u ) 2 


M n 





3 M 1 



Im 




2 

+ 


m v2 






^13 


2 

+ 


^12 


2 

+ 


/ill 


2 



(2.13) 



(2.14) 



~~ I tt \ 2 m ^3<5eff , 

8tt (#„) 

where the effective CP-violating phase <5 e fj is defined by 



<5cfT = 



As clearly seen from the above explicit expression, 5 e R is always less than one ||26| , but 
it is in general order one unless the phase of the Yukawa coupling h\ 3 is accidentally 
suppressed or the couplings have a inverted hierarchy, i.e., \h\ 2 \ 2 ^> (^2/^1/1) | ^13 1 2 or 
I ^11 1 2 3> (rri^ 3 / ) I /ii3 1 2 . From Eq. ( [2.13|) , the e\ parameter is given by 



El 



2.0 x 10 



-10 



M 1 



m v3 



10 6 GeV/ V0.05 eV 



5, 



cS 



(2.15) 



This relation is consistent with the one obtained in Ref. []60]]. Here, we have used {H u ) = 
174 GeV x sin/3, where tan/3 = (H u ) / (H d ). (H d is the Higgs field which couples to 
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down-type quarks.) Here and hereafter, we will take sin/3 ~ 1 for simplicity.^ As for the 
heaviest neutrino mass, we take ~ 0.05 eV as a typical value throughout this chapter, 
suggested from the atmospheric neutrino oscillation observed in the Super-Kamiokande 
experiments [|K]]. 

We should stress here that the e\ parameter has an explicit formula given in Eq. ( |2.15|) 
with <5 e ff in Eq. ( |2.14|) (as long as M\ -C M 2 , M 3 ), although in the literature it is sometimes 
treated as a free parameter. In particular, it is proportional to the right-handed neutrino 
mass M 1 for fixed values of m u3 and S e g(< 1). 

Let us mention one last point, the possibility of an enhancement of the asymmetry 
parameter e\. For the self-energy contribution f (x) in Eq. fl2.8|) (and Eq. ( |2.7|) ), we have 
assumed the masses of the right-handed neutrinos are not so degenerate, i.e., \M{ — M\\ ^> 
rVj- However, if the mass difference becomes as small as the decay width, \M~ M\\ ~ r^., 
one expects an enhancement of the self-energy contribution f s (x) [Bl|. Actually, it was 
shown that the asymmetry parameter e\ can reach its maximum value of e± ~ 0(1) for 
|Mj — Mi| ~ TnJ2 [|6l[| . (ei cannot be arbitrarily large. Notice that the lepton asymmetry 
vanishes in the limit where the right-handed neutrinos become exactly mass degenerate, 
since in this case CP- violating phases of the Yukawa couplings can be eliminated by a 
change of basis.) Nevertheless, it requires an extreme degeneracy of right-handed neutrino 
masses |Mj — M\\/M\ ~ h 2 /(8iT), and hence we do not consider this possibility in the 
following discussion. 



This is the case as long as tan/3 > 1. Even for tan/3 ~ 1, the final lepton asymmetry changes 
(increases) by a factor of 2. 
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2.2 Leptogenesis by thermally produced right-handed 
neutrino 



In this section, we briefly discuss the leptogenesis by thermally produced right-handed 
neutrinos. After first suggested by Fukugita and Yanagida 0, this scenario has been ex- 
tensively studied. (See, for example, Refs. |57|, |JT], |63|]. For reviews and references, 
see Refs. J22, ^kJ.) From the viewpoint of the production mechanism of the right-handed 



neutrino, compared with other production mechanisms, this scenario has a very attractive 
point that it requires no extra assumption to create the right-handed neutrinos, besides 
high enough temperature. 

Suppose that the right-handed (s) neutrinos N\ are produced thermally and become 
as abundant as in thermal equilibrium at temperature T > Mi. (Hereafter, we consider 
the lightest right-handed neutrino Ni, since lepton asymmetry produced by the decays of 
the heavier two right-handed neutrinos iVjjta) is likely to be erased before the iVi's decay.) 
Then the ratio of the total number density of the right-handed (s) neutrinos to the entropy 
density is simply given by the following number [see Eq. ( |B.19| )1:^ 

n Nl <_45<(3)W3 + , 1 _ 



s ~ s 2tt 4 g* 1 * N iJ ~ 240 

where = g~ = 2 are the numbers of degrees of freedom for right-handed (s)neutrinos 



and g* = 232.5 for the MSSM plus right-handed neutrino multiplet [see Eq. (|B.15| )1 



Boltzmann equations (a toy model) The evolution of the number density of the 
right-handed neutrino and that of the lepton number density til are described by coupled 
Boltzmann equations. In order to understand qualitative behaviors of the and n^, in 
particular to understand the "out-of-equilibrium condition," let us consider the following 
simple set of Boltzmann equations: 

d n Nl + 3Hn Nl = - 7 (iV a -> lip) - 7 (JVi -> lip*) 



dt 



+ _ Nl ) + 7 (7^* ^ Nl ) ; (2.17) 



^-m + SHni = 7(^1 -»■ lip) - j{lip -»• Nt) 



dt 

d_ 

dt 



- i{lip ^lip*)+ iQtp* -> lip) , (2. 18) 

nj + 3Hn T = 7 (JVi -> !</?*) - 7^* -> ^) 

+ 7'(^-W^)- 7 '(V->Z^), (2.19) 



Numerically the value in Eq. (2.16) is almost the same as l/g* 
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which describe the evolutions of the number densities of right-handed neutrino Aq, lepton 
I, and anti-lepton I. For a while, we will consider a non-SUSY case, for simplicity. The 
terms proportional to the Hubble parameter H describe the effect of the expansion of 
the universe. 7 (A — > B) denotes the reaction density, which means the rate of the 
process A — > B per unit time per unit volume. Thus, the 7-terms in the right-hand sides 
of the equations describe the change of the number densities due to the corresponding 
interactions. We have included only the three kinds of processes (Aq <-> lip, Aq <-» lip* 
and lip «-> lip*) which are necessary in the following discussion. (In actual calculation, 
one must include many other interactions. For other processes, which include the super- 



partners Aq, L and h u , see Ref. fl58| .) As for the primes in j'(l(p <-» lip*), we will give an 
explanation below. 

The reaction densities of the decay and inverse decay are given by 
^N^lip) = -(l + ei )(T Nl )n Nl , 
7ON1-V) = ^(l-e 1 )(T Nl )n Nl , 



7(^-^1) = J(l-6 1 )(r iVl )n^x^, 

Til 



7 (V-iVi) = -(l + eiXT^n^x-i. (2.20) 
/ rq 

Here, (Tni) denotes thermally averaged total decay rate of Aq. For low temperature 
T <C Mi, it is just given by the decay rate at rest frame: (r^) ~ r^, while for high 
temperature T ^> Mi, it is given by (r^) — (Mi/2T) due to the time-dilatation 
effect. 



From Eqs. ( |2.17| )- fl2.20p , the Boltzmann equations of Y Nl = n Nl /s and Y L = n L /s 



(rii — nj)/s are given by 



£y Nl = -(r^^-y^-k^jy^x^, (2.21) 



- 2 7 '(^ -»■ 7^*) + 2 7 '(V -> Zy>) . (2.22) 



Here, there is a subtle point one should take care |>|, As can be seen from Eq. fl2.20| ) 



the inverse decay processes (lip — > iVi and /(/?* — > Aq) produce a net lepton asymmetry 
with the same sign as the decay processes (Aq — > Zyj and Aq — > !<£>*) themselves. This 
can be seen by applying the CPT invariance to the matrix elements of those processes: 
M(N ± -> lip) = M{lip* -> Aq) and .M(Aq -> 7^*) = «M(/y> JVi). Therefore, if there 
would exist only those decay and inverse decay processes, lepton asymmetry would not 
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vanish even in thermal equilibrium. In other words, if we ignore the ^'(lip «-> lip*) terms, 
the right-hand side of the Eq. ( |2.22| ) would not vanish even for Yl = and Y^ 1 = Y^. 

On the other hand, there are lepton number violating two-body scattering processes; 
lip — > lip* and lip* — > lip. By using unitarity, one can show that these processes do not 
produce a net lepton asymmetry if we consider the N\ only as a virtual particle |65| , 
|66] . However, the Boltzmann equations already include Ni as a real particle. Thus one 
should subtract from the above two body scatterings the resonant s-channel contribution 
mediated by Aq (which is understood as an on-shell real particle) to avoid a double 
counting of reactions. After subtracting the resonant contribution, the reaction densities 
of the two-body scatterings in Eq. (|2.22|) leads to |65], |66| 



- 2i{iip -> lip*) + 2Y(V -> lip) 

-2e a (r Nl )n e ^-n L n^\{va'(lip 



lip*)) + -> iip) 



(2.23) 



where the primes mean that the contribution from resonant s-channel Ni exchange has 
been subtracted, (a is the cross section, v is the relative velocity between the initial 
particles and the bracket denotes the thermal average.) Notice that there appears a term 
— 2ei (r^) , which originates in the subtraction of the resonant contribution [f)5|| . After 
all, the Boltzmann equations are reduced to the following forms. 



d 
dt 



Y Nl 



d_ 

It 



Yr 



•(r.Y:)(Vv •- v-) l -(<^) ^ 





(2.24) 



^2.25) 



We see that, for Y^ 1 = Y^ (and Yl = 0), the lepton asymmetry produced by [decay 
and inverse decay] processes is canceled out by the lepton asymmetry from [two body 
scatterings minus resonance] processes, which ensures vanishing lepton asymmetry at 
equilibrium. 

Now let us discuss the evolutions of Y^ 1 and Yl- As can be seen from Eq. ( |2.25| ), 
a deviation of the number density of the Ni from its equilibrium value (Y^ ^ Y^) is 
mandatory in order to produce a net lepton asymmetry Yl ^ from Yl — 0. This can be 
realized when the temperature of the universe T cools down and becomes below the mass 
of the right-handed neutrino M\. To see the out-of-equilibrium condition for N l7 let us 
adopt a dimensionless variable z = M 1 /T instead of the cosmic time t. Then Eq. ( |2.24p 
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becomes 

dY Nl (T Nl ) 



Here, we have omitted the higher order term proportional to €\(Yl / Y^) , for simplicity. 
For T < Mi (z > 1), the equilibrium value of the Aq's abundance decreases due to 
the Boltzmann suppression oc exp(— Mi/T). At this stage, if the dimensionless prefactor 
(r^Vi) /zH is small enough, can no longer catches up the decreasing equilibrium 
value Y^. Namely, deviates from Y^, if (r^) / zH 1 for z^l. Thus, the out-of- 
equilibrium condition is roughly given by 

Tni <1 . (2.27) 



H 



T=M 1 



If this condition is satisfied, Y Nl > Y^ is realized, and a net lepton asymmetry is produced 
[see Eq. ( ggg )]. 



The ratio T^ 1 / H{T = Mi) is related to a mass parameter mi, which is defined as |22 



la 



Mi 

-3 „,r v , / F JVi 



~ 0.8 x l(T d eV x — — . (2.28) 
\ h \t=mJ 

Notice that the above formula looks similar to the neutrino mass given in Eq. (|2.2| ) but 
different from that. In terms of mf, the out-of-equilibrium condition in Eq. (|2.27|) is 
roughly equivalent to m^ ^ 10~ 3 eV. 

It is possible to show an important constraint on this parameter mj [ 57 |: 

ml > m v \ . (2.29) 

To show this, let us rewrite the m{ parameter in terms of the rotated Yukawa couplings 
hik defined in Eq. ( p.l2j ); 

^ = El^l 2 ^- (2-30) 

h 1 



On the other hand, from the seesaw formula Eq. fl2.2|) , we obtain 



m uj 5 jk = -J2^ifiik^-- ( 2 - 31 ) 



Let us define here a matrix as follows: 



^■Ej^-SL. (2.32) 



Mi7n V j 
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Then from Eq. ( [2.31| ) one can show that 

(X?X) 13 = Sij = , (2.33) 

while Eq. (|2.30|) gives rise to 

m{ = ^m uk \Xi k \ 2 

k 

> mm{m uj }J2\X lk \ 2 



> min{m^} 



Ik 



min {m U j} . (2.34) 



In the last equation, we have used Eq. ( p.33|) . Therefore, the mi parameter is bounded 
from below as > m v \ WA - 



Lepton asymmetry Now let us discuss the amount of produced lepton asymmetry 
in the present scenario. If the decaying right-handed neutrinos are as abundant as in 
thermal equilibrium [see Eq. (|2.16|) 1, and if there is no wash-out process of the produced 
lepton asymmetry afterwards, the final lepton asymmetry would be given by a very simple 
formula, 

s s 240 



with the asymmetry parameter e± given in Eq. (|2.15 ). In the actual case, however, a 
suppression factor k should be multiplied: 

The suppression factor k represents two effects. The first one is the wash-out effect of 
the produced lepton asymmetry (and the strength of the deviation from thermal equi- 
librium). If the interactions of right-handed neutrinos N\ are too strong, the produced 
lepton asymmetry would be washed out by the interactions mediated by the N\ itself. 
[See the second term in Eq. ( [2.251 ).] Weakness of the interaction is also required from the 
out-of-equilibrium condition for N\ in Eq. ( |2.27|) . 

The second effect represented by k is the efficiency of the production of the right- 
handed (s)neutrinos N\. Although we have assumed that Ni are produced as abundant 
as in thermal equilibrium ~ for T > Mi, if the Nx's interactions are too weak, the 
thermal scatterings cannot produce enough amount of N\ and hence the number density 
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njVi cannot become as abundant as that in thermal equilibrium, .Q Notice that those 
two effects have double-edged behavior, i.e., the right-handed neutrino N\ should have 
couplings hi a with intermediate strength. 

In order to determine the precise value of the suppression factor k, one has to solve the 
Boltzmann equations including production, decay and inverse decay of the right-handed 
(s)neutrinos, and all the relevant lepton-number violating (and conserving) scatterings 



of lepton fields [32, 22]. A detailed numerical calculation solving the coupled Boltzmann 



equations in the case of SUSY has been done in Ref. p8| . (See also Ref. [Jj2| and references 



therein.) It was shown that both of the production rate of the right-handed (s) neutrinos 
and the rate of the wash-out process of the lepton asymmetry are proportional to the mass 
parameter m{ defined in Eq. ( |2.28| ). It turns out that the suppression factor k becomes 
as large as k ~ 0.05-0.3 for [|58|1 

10~ 5 eV<m^<5 x 10~ 3 eV. (2.37) 

For a smaller value of m~i, the Ni is not produced enough and generated lepton asymmetry 
is suppressed. On the other hand, for a larger value of m~i, enough amount of N\ is 
produced but wash-out of the lepton asymmetry is too strong, and hence final lepton 
asymmetry is reduced. 

It is quite interesting to observe that the favored range of the mass parameter fn[ 



is just below the neutrino mass scale observed by the atmospheric and solar 
(+KamLAND |fl2"| ) neutrino-oscillation experiments, m u ~ O(10 _3 )-O(10 _1 ) eV. We 



should note that although m 1 is not directly related to the neutrino masses m ui , they 
are still indirectly related. (For instance, if we adopt a Froggatt-Nielsen (FN) model [|68[ 



which will be discussed in Sec. [2.3. 2| , mi is estimated as mi ~ 0{m uZ ) ~ (9(0.05 eV). 



Although this is a bit larger than the range in Eq. (|2.37|) , it can be consistent when we 



include 0{1) ambiguities in the FN model. See also Ref. |69|], where a more general case 
was discussed for two generations of neutrinos.) 

After being produced, a part of lepton asymmetry is immediately converted to the 
baryon asymmetry [[| via the sphaleron effect || discussed in Sec. |1.4| . Then the present 
baryon asymmetry is given by 

^ = 0.35x^ 



s 



* 0.3 x 10- ( ^) x ) S m , (2.38) 



0.1/V10 9 GeV/ V0.05eV 



6 One can consider that enough amount of right-handed neutrinos exist from the beginning and discuss 
their decay in thermal background. We do not consider such a case in this section, however, and assume 
that the right-handed neutrinos are produced by thermal scatterings from ~ [p2[ . Notice that as 
long as an inflationary epoch is assumed, right-handed neutrinos should be produced at some stage from 
n Nl ~ 0. 
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where we have used Eqs. ( |2.15| ) and ( |2.36| ). Therefore, the present baryon asymmetry 
ub/s ~ (0.4-1) x 1CT 10 is naturally explained with the mass of the lightest right-handed 
neutrino Mi ~ 10 9 -10 10 GeV, for k ~ 0.05-0.3 and 5 cS ~ 0(1). 

In order to produce enough amount of right-handed neutrinos (i.e., n Nl ~ n e ^), the 
temperature of the universe should be higher than their mass M\. This leads to a lower 
bound on the reheating temperature as Tr^C(10 10 ) GeV. Thus the overproduction 
of gravitinos might cause a difficulty depending on the gravitino mass, as discussed in 
Sec. |Q , If the gravitino is unstable, it must be relatively heavy, m 3 / 2 ~ (a few) TeV. 
When the gravitino is stable, a consistent thermal history can be obtained with gravitino 
mass 10-100 GeV, avoiding the problem of the decays of next-to-lightest SUSY particles 
after the big-bang nucleosynthesis [^8|.[] 

Finally, we comment on the absolute upper bound on the neutrino masses from thermal 
leptogenesis, which has been shown recently in Ref. [Q. The crucial observation in 
Ref. is that the suppression factor k is determined only by three parameters k = K,(rfTi, 
Mi, m), and hence the final baryon asymmetry in thermal leptogenesis depends only on 



four parameters: e±, rax, Mi, and m. [See Eq. (|2.36| ).l Here, m = ym^ + m^ 2 + m 



l u3- 



Furthermore, the maximal value of e± is also determined by Mi and m, ei < e™ ax (Mi, m), 
once the mass squared differences of atmospheric and solar neutrino oscillations are given 
[see Eq. ( |2.14| ) and Eq. ( [2.15D 1 [f26j |70| , [63| . Thus, the maximal baryon asymmetry is 



determined only by the set of three parameters (m^, Mi, m). By solving the Boltzmann 
equations for different points in this parameter space (m^, M 1; m), and using the bound 
Si > m„i, it has been shown in Ref. |)3) that the maximal baryon asymmetry can be 
larger than the empirical value only if 

m<0.30eV. (2.39) 

Therefore, if the baryon asymmetry in the present universe was indeed generated by 
thermal leptogenesis, we have stringent constraints on the absolute neutrino masses, 

m^<0.18eV. (2.40) 



'See also Ref. 
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2.3 Leptogenesis in inflaton decay 



In this section,^ we discuss the leptogenesis scenario where the right-handed neutrino 
is produced non-thermally in inflaton decays [p3]] . We will find that this scenario is 
fully consistent with existing various SUSY inflation models such as hybrid, new, and 
topological inflation models.0 

The crucial difference between this scenario and the case of the thermally produced 
right-handed neutrino discussed in the previous section is that, the heavy right-handed 
neutrino N can be produced with relatively low reheating temperatures Tr of the inflation. 
We find that the required baryon asymmetry can be obtained even for Tr ~ 10 6 GeV in 
some of the inflation models, and hence there is no cosmological gravitino problem in the 
interesting wide region of the gravitino mass — 10 MeV-10 TeV.[^] (See Sec. |1.5| .) 

On the other hand, the amount of the produced lepton asymmetry (and hence baryon 
asymmetry) in the present scenario crucially depends on the physics of the inflation, such 
as the mass of the inflaton m x and the reheating temperature Tr. Therefore, detailed 
analyses on the inflation models are necessary. 



In Sec. p.3.1| , we calculate the amount of the resultant lepton asymmetry in inflaton 
decay. Then we introduce a Froggatt-Nielsen model in Sec. [2.3.2| to estimate the mass 
(and decay rate) of the right-handed neutrino. The subsequent subsections are devoted to 
each SUSY inflation model. The leptogenesis in a hybrid inflationary universe is discussed 
in Sec. |2.3.3| , where we consider two different types of SUSY hybrid inflation models. In 
Sec. |2.3.4| we discuss the leptogenesis in a SUSY new inflation model. The case of a SUSY 
topological inflation is considered in Sec. p. 3. 5 . Finally, we will briefly comment on the 
production of right-handed neutrinos at preheating |T4| in Sec. [2.3.6| . 



2.3.1 Lepton asymmetry 

Let us first estimate the produced lepton asymmetry in the present scenario. The result 
obtained here is a generic one, which can be applied to all the inflation models discussed 
in subsequent subsections. 

After the end of inflation, the inflaton decays into light particles and the energy of 
the inflaton is transferred into the thermal bath. Then it is very plausible that the right- 
handed neutrino is also produced in inflaton decay, if its decay channel is kinematically 



8 This section is based on the works in a collaboration with T. Asaka, M. Kawasaki, and 
T. Yanagida @. 

9 Leptogenesis in a "natural chaotic inflation model" |72| was also investigated in Ref. [7^ ]. 
10 Here, we consider both of the unstable and stable gravitino. 
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allowed. We assume that the inflaton decays into two right-handed neutrinos, which leads 
to the following constraint: 



m x > 2M 1 , (2.41) 

where m x is the inflaton mass. Here and hereafter, we consider only the N\ decay, 
provided that the mass Mi is much smaller than the others (Mi <C M 2 , M 3 ). As discussed 
in Sec. |2J], the decay of N\ into leptons and anti-leptons produces a lepton asymmetry. 
We will consider the case where the right-handed neutrino is heavy enough compared 
with the reheating temperature, i.e., Mi ^ Tr. In this case the produced right-handed 
neutrino N\ is always out of thermal equilibrium and it behaves like frozen-out, relativistic 
particle with energy — m x /2. The ratio of the number density of the right-handed 
(s)neutrinos to the entropy density s is then given by 

nNi Prad n x n Nl 
— - ~ - x — x — - 



s s p x n x 

= -T R x — x2B r = \B r -^, (2.42) 
4 m x 2 m x 

where p ra d is the energy density of the radiation just after the reheating process completes, 
and n x and p x = rn x n x ~ p ra d are the number and energy density of the inflaton just 
before its decay, respectively. B r = B r (x — * NiN\) denotes the branching ratio of the 
inflaton decay into Ni channel. 

As we will see, when we adopt the model which will be introduced in the next sub- 
section, the Ni decays immediately after produced by the inflaton decays. The lepton 
asymmetry is then given by 

n L n Nl 

— - ei 

s s 



3 x 10- 10 B r -±±- ^ x 5eff , (2.43) 




where we have used the asymmetry parameter e\ given in Eq. ( |2.15| ). Notice that there 
is no wash-out effect of the produced lepton asymmetry,0 since N\ is out of thermal 
equilibrium. 

Let us mention one point here. Although we will assume an explicit model in the 
next subsection to make the subsequent discussion concrete, the formula of the lepton 
asymmetry in Eq. ( p. 43 ) is a generic one, unless Mi ^$Tr or N\ has a very long lifetime 



and dominates the energy density of the universe before it decays. If Mi J$ Tr, the N\ 



11 Lepton-number violating 2-body scatterings mediated by TV; are out of thermal equilibrium as long 
as <(10 keV) 2 (10 9 GeV/T fl ) 
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is produced more or less by thermal scatterings, and the resultant lepton asymmetry 
is expected to be reduced to the case discussed in Sec. p72| . As for the case where Ni 
dominates the energy density of the universe before its decay, we give a brief discussion 
in Appendix |2.A. 



As discussed in Sec. |1.4| , the sphaleron process converts this lepton asymmetry into 
baryon asymmetry as %/s ~ 0.35 x ul/s. In order to explain the observed baryon 
asymmetry ur/s ~ (0.4-1) x 10~ 10 , we should have lepton asymmetry 



— ~ (1 -3) x 10~ 10 . (2.44) 



From Eqs. ( |2.41 ) and fl2.43| ), we can see that the reheating temperature of the inflation Tr 
is bounded from below as Tr ^ 10 6 GeV, since otherwise the produced lepton asymmetry 
is too small as n L /s < 10~ 10 . 



2.3.2 Froggatt-Nielsen model for neutrino 

As we have just shown, the amount of generated lepton asymmetry depends on the mass 
of the inflaton m x and the reheating temperature Tr (and the branching ratio B r ), as 
well as the mass of the right-handed neutrino Mi. Before discussing the inflaton sector, 
we here introduce the Froggatt-Nielsen (FN) model in order to settle the mass (and 



decay rate) of the right-handed neutrino and to make the discussion in the subsequent 
subsections concrete. 

The FN model is one of the most attractive framework for explaining the observed 
hierarchies in quark and lepton mass matrices. This model is based on an U(1)fn sym- 
metry that is broken by the vacuum-expectation value of S, (S) ^ 0. Here H is a gauge 
singlet field carrying the FN charge Q= = — 1. Then, all Yukawa couplings are realized as 
nonrenormalizable interactions including S, and are given by the following form; 

= y %] e q * +q i ^jH u{d) , (2.45) 

where y^j are 0(1) coupling constants, Qi are the FN charges of various chiral superfields 
ipi and e = (S) /Mq. Here, we have assumed that Higgs doublets H u and Hd (which couple 
to up-type and down-type quarks, respectively) have zero FN charges. The observed mass 
hierarchies for quarks and charged leptons are well explained by taking suitable FN charges 
for them. For instance, we assign FN charges (a + 1, a, a) for lepton doublets Lj, while 
giving charges (2, 1, 0) to the right-handed charged leptons with e ~ 0.05-0.1 |75|. |77|. 
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■0t 


L\ L2 L% 


E\ E2 E3 


Ni N 2 N 3 






a + 1 a a 


2 1 


deb 


-1 



Table 2.1: The Froggatt-Nielsen charges of various supermultiplets. a 
assume d > c > b > 0, i.e., M x < M 2 , M 3 . 



or 1. We 



The charges are listed in Table. |2.1| . We will take a = or 1 according to Ref. fTSfl . The 
Charges of the quarks are also determined in the framework of the SU(5) grand unified 



theory WE 



Let us apply the above mechanism to the neutrino sector as well. The mass matrix 
for the heavy right-handed neutrinos iVj is given by 



M, 



Rij 



t Jl-Q, M( 



• 



(2.46) 



where M represents some right-handed neutrino mass scale and are coupling constants 



of order unity like jjij. Charges for right-handed neutrinos are also found in Table. \27L 
We take d > c > b > 0, i.e., M x <C M 2 , M 3 . Hereafter, we will take a basis where the 
mass matrix for the charged leptons is diagonal.^ Then, the neutrino Dirac mass matrix 
mu and the right-handed neutrino mass matrix Mr are given by the following forms; 



M. 



R 




( 



-a+1 







where {rno)ia is defined as W 
then given by 



Uii Vi2 m 

2/21 V22 y™. 

2/31 V32 V33 

Cn ^12 £13 
C12 ^22 C23 

\ £13 £23 £33 / 

e, /i, r). The neutrino mass matrix is 



V 

/ e d 


V 








e c 











(2.47) 



(m D ) ia NiL a (a 



m 



-2a 



(HuY 



( e 



Mn 



-2a 



\H u y 



M 



1 

V 1 

0(e 2 ) 0(e) O(e) 

0{e) 0(1) 0(1) 

0(e) 0(1) 0(1) 



-1 




(2.48) 



12 Hcre, one might wonder if the mixing matrix from the charged lepton sector would change the 
discussion, since the mass matrix for the charged leptons has off-diagonal elements in the above FN 
mechanism. However, the correction from this effect yields higher order terms in e, and hence we can 
safely neglect it. 
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As shown in Ref. [[T!J, this mass matrix can naturally lead to a large v^—v^ mixing angle, 
which is suggested from the atmospheric neutrino oscillation [1C]. It is remarkable that the 
FN charges of the right-handed neutrinos are completely canceled out in the neutrino mass 
matrix in Eq. ( |2.48| ) and hence the hierarchy in the neutrino mass matrix is determined 
only by the charges of the lepton doublets, (a + 1, a, a). 

In this model, the mass scale of the right-handed neutrino Mq is estimated as 
M ~ ~ 2a{Hu)2 



1-2 n„\- <•,., ,. - i (2.49) 



{ 0.6 x 10 15 GeV for a = 
\ (1 - 6) X 10 12 GeV for a = 1 . 

Here, we have used m u3 ~ 0.05 eV and e ~ 0.05 - 0.1. Then the mass of the lightest 
right-handed neutrino Mi is given by 

Ma ~ M e 2d 

j (1 - 6) x 10 12 GeV for a + d = 1 
\ (0.3 - 6) x 10 10 GeV for a + d = 2 

In the following analysis, we only consider the case a + d = 1 or a + d = 2. Thus, our 
assumption of out-of-equilibrium condition M 1 ^Tr is justified as far as Tr^S 10 9 GeV. 
On the other hand, the total decay width of the Ni, T^, is given by 

T Nl ~ -L £ 2 ^)M 1 

47T 



io r<„\- I'.,,- ., ./ - ■-> (2.50) 



J (0.3 - 5) x 10 9 GeV for a + d = 1 , 

j (0.02 - 5) x 10 5 GeV for a + d = 2. ^ ' 

In deriving the formula of the lepton asymmetry in Eq. ( [2 .43] ) , we have assumed that the 
Ni decays immediately after produced in inflaton decay, which corresponds to ^> T x . 
(T x is the decay rate of the inflaton x-) m terms of the reheating temperature Tr ~ 
(7r 2 (yf*/90) _1 / 4 y / MGr x , this requirement corresponds to 

T fi «2x 10 10 GeV f in F"0 ^ ■ (2.52) 



V10 3 GeV, 

Therefore, our assumption is again guaranteed for Tr ^ 10 9 GeV. 

Finally, let us estimate the asymmetric parameter e\ in the present FN model. From 
Eq. (p|), it is given by || 



1 



fi 



8vr (e a + d f 

3 c 2(a+d) 

8tt 

3 m v \M\ 
8^ (H u ) 2 



f a+d a+b^ z b 



■ 6 



£ 2b 



(2.53) 
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We see that the above estimation agrees with the general result for e\ in Eq. ( |2.13| ), with 
5 e fj ~ 1. Thus, we can use the formula of the lepton asymmetry in Eq. ( 2.43Q consistently. 



2.3.3 Hybrid inflation 

In this subsection we perform a detailed analysis on hybrid inflation models and examine 
whether they can produce sufficient lepton asymmetry to account for the baryon asym- 
metry in the present universe, avoiding the overproduction of the gravitinos. 

2. 3. 3. a hybrid inflation with a B — L symmetry 

Before discussing hybrid inflation models, let us first show a particle-physics model for 
the heavy right-handed neutrinos iVj. A simple extension of the SUSY standard model 
is given by considering a gauged \](1)b-l symmetry, in which, as mentioned in Sec. [j~T| , 
right-handed neutrinos iVj are necessary to cancel B — L gauge anomaly. We introduce 
standard-model gauge-singlet supermultiplets \l/ and \l> carrying B — L charges +2 and 
—2, respectively, and suppose that the B — L symmetry is spontaneously broken by the 
condensations (^} = at high energies. = (ty^ is required from the D-term 

flatness condition of the U(I)b-l-) Then the heavy neutrinos N iy which carry B — L 
charge —1, acquire Majorana masses through the following superpotential^ 

W = y i 9N i N i . (2.54) 

Here, we have assumed that \1/ and \1/ have zero FN charges. Then the right-handed 
neutrino scale M in Eq. ( |2.49| ) is given by 

M = (*) = (¥) . (2.55) 

Thus, the right-handed neutrino mass scale Mq derived in the FN model in Eq. ( |2.49| ) is 
explained in term of the B — L breaking scale about (9(10 15 ) GeV or 0(1O 12 ) GeV. 

A superpotential causing the B — L breaking is given by 

W = $ (A^-/x 2 ) , (2.56) 

where $ is a gauge-singlet supermultiplet, A a coupling constant and fi a dimensionful 
mass parameter. (We will take a basis where /i 2 and A are real and positive by using 



Leptogenesis in the hybrid inflation discussed in Ref. ]7q | assumes a nonrenormalizable superpotential 
W = (g' i /2)N i N i ^^, where the B - L charge of * is taken to be +1. 
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the phase rotations of $ and Notice that this superpotential possesses a U(l) R- 

symmetry where the $ and have U(1)r charges 2 and 0, respectively. The potential 
for scalar components of the supermultiplets \1/ and \I/ is given by, in supergravity, 



V 



,K/M* 



d 2 K 
dZ~dZ* 



D Z WD Z *W* -3 



1 



\W\ 



+ terms 



;2.57) 



where 



D Z W 



dW 1 &Ff 
+ — — W, 



Ml dZ, 



(2.56 



and Zj denote scalar components of supermultiplets. We assume the /^-invariant Kahler 
potential for $, \l/ and \I/ 



|<|)| 2 _|_ |i[/| 2 _)_ |v[/| 2 _|_ 



4M£ 



l$l 



(2.59) 



where the ellipsis denotes higher-order terms which we neglect in the present analysis. 
From the above scalar potential, we have the following SUSY-invariant vacuum (hereafter, 
we denote the scalar components of the supermultiplets $ and \l/ by the same symbols as 
the corresponding supermultiplets): 



(*> = o, 



(2.60) 



where we have chosen to be real and positive by using the U(l) b-l rotation. 

It is quite interesting to observe that the superpotential (|2.56|) is nothing but the one 
proposed in Refs. |79, [80[ for a SUSY hybrid inflation model. In this context, the real 
part of $ is identified with the inflaton field pj\f2. The scalar potential is minimized at 
— ^ — when (p is larger than the following critical value: 



(2.61) 



and hybrid inflation occurs for cpc < <p ^ M G and k = — k% > 0. Including one-loop 
corrections [79[|, the potential for the inflaton <p is given by, for p > pc-. 



V 



<P 



k 

1+ 2\M G 



16 V / \M G 



+ 



1287T 2 



2 ^ ln (w) + (^ 2 -^) Zln ( 1 



2A 2 



V 2 , 



(2.62) 
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where A denotes the renormalization scale. Here, we have included the higher order terms 
which come from the Kahler potential ( |2.59| ) up to fj, 4 ((p/Ma) 4 term, since the initial value 
of the inflaton field is close to 

Let us now discuss the inflation dynamics. The slow-roll conditions for inflation are 
given by p8 

2 

< 1, 




y < 1, (2-63) 

where the prime denotes the derivative with the inflaton field <p. These conditions are 
satisfied when (p > tpc, A < 1 and < k < 1. Therefore, while the inflaton p> rolls down 
along the potential in Eq. (|2.62|) from <pi (<p c < V 9 / ~ Ms) to ipc, the vacuum energy /i 4 



of the potential dominates the energy of the universe and hence the hybrid inflation takes 
place |81] . 



After the inflation ends, the vacuum energy is transferred into the energies of the 
coherent oscillations of the following two fields: the inflaton p> and a scalar field £ = 
(\? + \&)/ y/2. Notice that the inflaton <p forms a massive supermultiplet together with the 
£ field in the vacuum in Eq. (|2.60|) , whose masses are given by 



= ms = v2A/i 

= V2 X ($) . (2.64) 

The radiations of the universe are produced by the decays of the tp and/or £ field. 
In order to estimate the reheating temperature Tr we should know the total decay rates 
of these scalar fields. Through the interactions in the superpotentials in Eqs. (|2.54|) 
and (|2.56|) , the inflaton <p decays into scalar components Ni of the Ni supermultiplet, if 
kinematically allowed, with the rate 



_L_ (MA 2 ( 4M 1 2 \ 1/2 
64n m *{m) [ ml) 



(2.65) 



while the £ field decays into scalar (JVi) and fermionic (Ni) components of the iVi with 
the rates 

r ^ m) = ±. ±(m 2 ( 1 _^y\ ( , 66) 

16tt m s V (*) / V m% J 



14 However, we have neglected the /i 4 (<p / Mg) 4 terms coming from the higher order interactions in the 
Kahler potential (e.g., K — k'\$>\ 6 /M G ), for simplicity. 
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Here, we have used Mi = £1 and m^, = m s = a/2 A ( 1 i r ). 

Since the £ field has a non-zero vacuum expectation value, it can decay also through 
nonrenormalizable interactions in the Kahler potential, 

^ = E^IW, (2.68) 

where ipt denote supermultiplets of the SUSY standard-model particles including the 
right-handed neutrinos Ni, and c, are coupling constants of order unity. Then the partial 
decay rate through these interactions is estimated as 



2 



r(E - tfik) - j^C I ^ 1 m|, (2.69) 

where C = yjc^ is a parameter of order unity. In the following analysis we take C ~ 1, 
for simplicity. Comparing this rate fl2.69|) with the rates in Eqs. (|2.66|) and ( 2.671) , we see 



that the total decay rate of the £ field is determined by that in Eq. (|2.67|) , i.e., 



r(E -> JViJVj) , (2.70) 
since T(S -> AqiVi) > T(S -> JV^) > T(E -> ^) for 2 Mi < m s < M X (M G / W) 2 - 



We assume 2Mi < = as discussed in Sec. [2.3. 1| . Then the inflaton (p and the E 
field decay at the almost same time because of ~ r^. Thus, the reheating temperature 
Tr is estimated by the total decay width of the E field as 



7T 



2 \ -1/4 
.2 



1R "{90 7T ) V F ^g. (2.71) 

Here it should be noted that the branching ratio of the decay rate of the E field and that 
of the inflaton tp into two N x (Ni) is almost unity, and hence we have B r ~ 1 in Eq. fl2.43|) . 

The above hybrid inflation must explain the following two observations:^ 

(i) the e-hold number N e of the present horizon, and 

(ii) the density fluctuations observed by the cosmic background explorer (COBE) satel- 
lite 



15 



The spectrum index in the hybrid inflation is predicted to be scale- invariant, n s ~ 1. 
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First, while the inflaton ip rolls down along the potential ( |2.62| ) from ipN e to ip G , the scale 
factor of the universe increases by e e . This e-fold number N e is given by 

In order to explain the present horizon scale, the e-fold number should be 

^= 67+ H^r)' (273) 

where H in{ = \/V / (\/3M G ) ~ p 2 /(\/3M G ) denotes the Hubble parameter during the 
inflation. Next, the amplitude of the primordial density fluctuations Sp/p predicted by 
the hybrid inflation, 

sj, i v"Hf N .) (2 74) 

should be normalized by the data on anisotropies of the cosmic microwave background 



radiation observed by the COBE satellite f82 |, which gives the following constraint: 



M G \V'(cp Ne )\ 



5.3 xlO" 4 . (2.75) 



The parameters of the inflation are determined by the constraints in Eqs. ( |2.73|) 
and (|2.75|) . We have performed numerical calculations using the full scalar potential 



in Eq. ( |2.62| ). The scale p of the hybrid inflation is determined for given couplings A 



and k, and it is shown in Fig. p.lj .p 5 ! Here and hereafter, we exclude the region where 
y?7v e ~ Mq, because in that region the higher order terms of (y/Mg)" become large and 
our effective treatment of the inflaton potential in Eq. (|2.62|) would be invalid. 

In Fig. |2.2] and Fig. |2.3| , we also show the B — L breaking scale and the inflaton 
mass m^p. It is interesting that, as shown in Fig. |2.2| , the scale of the B — L breaking is 



predicted as ~ (1~5) x 10 15 GeV in a wide parameter region, 10~ 6 ^$ A ^$ 10~ 2 and 
k 10~ 3 , which is consistent with Mq ~ 10 15 GeV (i.e., a = 0) derived from the observed 
neutrino mass [see Eq. (|2.49| )1. On the other hand, the lower value of the B — L breaking 
scale of Mo ~ 10 12 GeV (a = 1) cannot be obtained in the present hybrid inflation model. 

The reheating temperature Tr depends on the mass of the right-handed neutrino, 
Mi, since the decay rate of the £ (and ip) depends on Mi. We take Mi ~ 10 12 GeV 
(a + d = 1) and Mi ~ 3 x 10 9 GeV (a + d = 2) for representation. The obtained reheating 
temperatures Tr are shown in Fig. fZJ and Fig. 



16 Although the scale /x slightly depends on the mass of the right-handed neutrino Mi, we find no sizable 
difference in the scale fi between the cases of a + d = 1 and a + d = 2, which are discussed below. This 



is because Mi enters in the calculation only through the small correction of ln(T^) in Eq. ( 2.73| ) 
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Figure 2.1: The contour lines of the scale /i in the hybrid inflation model. The contour 
lines are shown by the dotted lines and corresponding values of \x are represented in unit 
of GeV. The upper bound on k from the requirement <Pn £ < Mq is shown by the thick 
solid line. 

Here, we require the reheating temperatures Tr to be lower than 10 8 GeV to avoid 
overproduction of the gravitinos in a relatively wide range of gravitino mass (see Sec. |1.5| ). 
It is found that, for the region of the inflaton mass = ^> 2 Mi, the reheating 
temperature Tr ^ 10 8 GeV is obtained only for the case Mi ~3x 10 9 GeV (a + d = 2). 
In the case of Mi ~ 10 12 GeV (a + d = 1), the desired low reheating temperature is 
obtained for the region m v < 2Mi because the decay into Ni is kinematically forbidden 
and the decay rate is determined by the suppressed decay width in Eq. ( |2.69| ). However, 
such cases are not interesting since the Ni are not produced in the ip and E decays and 
leptogenesis does not take place.Q 

Now let us examine whether the leptogenesis works well or not in the above hybrid 
inflation model. Since the heavy Majorana neutrinos Ni are produced in the decays of 
the inflaton ip and the E field, the mass of the inflaton should satisfy = > 2Mi. 



As derived in Sec. |2.3.1| , the ratio of the produced lepton number to the entropy n L /s is 
given by Eq. ( [2.43 ). Here, we stress again that the branching ratio is automatically given 
by B r ~ 1 in the present hybrid inflation. We show in Fig. |2.6| and Fig. |2.7| the Ul/s for 
the cases a + d — 1 and a + d — 2, respectively. 



17 See, however, also Sec. 2.3.6 
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Figure 2.2: The contour lines of the fi?) in the hybrid inflation model. The contour lines 
are shown by the dotted lines and corresponding values of are represented in unit of 
10 15 GeV. 
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Figure 2.3: The contour lines of the inflaton mass in the hybrid inflation model. The 
contour lines are shown by the dotted lines and corresponding values of m v are represented 
in unit of GeV. 
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Figure 2.4: The contour lines of the reheating temperature Tr in the hybrid inflation 
model for the case Mi ~ 10 12 GeV (a + d = 1). The contour lines are shown by the 
dotted lines and corresponding values of Tr are represented in unit of GeV. The upper 
bound on k from the requirement ifN e < Mq is shown by the thick solid line. The k 
yielding m v = ms = 2M\ is also shown by the thick solid line. 
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Figure 2.5: The same as Fig. |]4|, but for the case M 1 ~ 3 x 10 9 GeV (a + d — 2). 
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Figure 2.6: The contour lines of the lepton asymmetry n^/s in the hybrid inflation model 
for the case Mi ~ 10 12 GeV (a+d = 1). The contour lines are shown by the thin solid lines 
and corresponding values of n^/s are represented. We also show the contour lines of the 
reheating temperature by the dotted lines and corresponding values of Tr are represented 
in unit of GeV. The upper bound on k from the requirement <fN e < Mq and the lower 
bound on k from = ms > 2M\ are both shown by the thick solid lines. 
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Figure 2.7: The same as Fig. |1| but for the case M 1 ~ 3 x 10 9 GeV (a + d = 2). 
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First, we consider the case of Mi ~ 10 12 GeV (a + d = 1). We find from Fig. |2l] that 



the lepton asymmetry enough to explain the present baryon asymmetry can be generated 
in a wide parameter region. However, we have a too high reheating temperature of 
Tr ~ 10 9 -10 12 GeV as mentioned before. Therefore, only the small region of ~ 2 Mi, 
where the reheating temperature becomes as small as Tr ~ 10 6 -10 8 GeV due to the phase 
space suppression is free from the cosmological gravitino problem for a wide range of 
gravitino mass. We can obtain in this very narrow region the required lepton asymmetry 
of Ul/s ~ (1-3) x 10~ 10 to account for the present baryon asymmetry. 

Next, we consider the case of Mi ~ 3 x 10 9 GeV (a + d = 2). It is found from Fig. 2.7 



that the required lepton asymmetry of ni/s ~ 10~ 10 as well as the low enough reheating 
temperature of Tr ~ 10 7 -10 8 GeV are naturally offered in the region of k ^5 10 -3 and 
A ~ 10~ 6 -10~ 3 . Therefore, we conclude that the hybrid inflation with Mi ~3x 10 9 GeV 
can produce a sufficient baryon asymmetry, giving a reheating temperature low enough 
to solve the cosmological gravitino problem. However, even lower reheating temperature 
of Tr ~ 10 6 GeV, which is required to avoid the cosmological difficulty for a very wide 
range of gravitino mass m 3 / 2 — 10 MeV-10 TeV, is achieved only in the narrow parameter 
region of = mj ~ 2Mi, where Tr is reduced due to the phase suppression as in the 
previous case (a + d = 1). 

2.3.3.b hybrid inflation without a B — L symmetry 

We have, so far, identified the U(l) gauge symmetry in the hybrid inflation model with 
the B — L symmetry. We now consider the case where the U(l) symmetry is not related 
to the B — L symmetry and even completely decoupled from the SUSY standard-model 
sector. The role of the U(l) gauge symmetry here is only to eliminate an unwanted flat 
direction in the superpotential in Eq. ( |2.56| ). 

We reanalyze the leptogenesis in hybrid inflation in the absence of the superpotential 
in Eq. ( p. 54 ). In this case the E field decays through the nonrenormalizable interactions 



in Eq. (|2.68 ). On the other hand, the decay of the inflaton (p is much suppressed due to 
the absence of the interaction in Eq. ( |2.54j ). Thus we introduce a new interaction in the 
Kahler potential as 

K = -£-$*H u H d + H.c. . (2.76) 

Through this interaction the inflaton tp can decay faster than the S field for the coupling 
h^C ((\I / ) /Mg) 2 — 10~ 6 , and the reheating temperature Tr is given by the decay of the 
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In this region, we should include the decay rates in Eqs. (2.66) and (2.69) in estimating Xj 
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E field. (Notice that the reheating temperature is determined by the slower decay rate.) 
Since the decay rate of the £ field in Eq. (|2.69| ) is very small compared with Eq. ( |2.67|) , 
the reheating temperature Tr becomes much lower than in the previous model. We show 



the obtained reheating temperature Tr in Fig. |2Tg . 

The inflation dynamics is almost the same as in the previous hybrid inflation model. 
The differences in the results for parameters, //, ($?) and m v , only come from the the 



small correction of ln(Tft) in Eq. ( 2.73 ), and we find no sizable difference in these values 
between the present and the previous models. 

Now let us turn to estimate the lepton asymmetry in this hybrid inflation model. First, 
we find that too small lepton asymmetry is obtained for Mi ~ 3 x 10 9 GeV (a + d — 2). 
This is because the Tr and Mi are too low to produce enough lepton asymmetry [see 



Eq. ( p. 43 )1- Hence, we concentrate on the case of Mi ~ 10 12 GeV (a + d = 1) in the 



following analysis. 

We show the obtained lepton asymmetry til/s in Fig. |2.9| . Here, we have assumed B r = 
1. Notice that the E field decays not only into the heavy neutrinos iV\ but also into the 
SUSY standard-model particles through the nonrenormalizable interactions in Eq. ( |2.68|) 



and hence B r ~ 1 is not automatic in this model. Therefore, the lepton asymmetry 



shown in Fig. |2J)| is understood as a maximal value. It is found that the required lepton 
asymmetry to account for the empirical baryon asymmetry can be generated in a wide 
parameter region k ^ 10~ 2 and A ~ 10~ 3 - 10~ 2 with the reheating temperature of Tr ~ 
10 6 -10 8 GeV. It is remarkable that we can obtain rii/s ~ 10 10 and Tr ~ 10 6 GeV 
simultaneously for A ~ 10~ 3 and /c^SlO -3 . The overproduction of gravitinos can be 
avoided in the full gravitino mass region of 777.3/2 — 10 MeV-10 TeV with such a low 
reheating temperature Tr ~ 10 6 GeV.0 

2.3.4 New inflation 

In the previous subsection, we have seen that the SUSY hybrid inflation can successfully 
produce the lepton asymmetry to explain the baryon number in the present universe, 
even with low reheating temperatures of Tr ~ 10 6 -10 8 GeV. In was found, however, 
that relatively small couplings of k ^ 10~ 2 and A ^5 10~ 2 are necessary to realize such low 
reheating temperatures. Although this is not so problematic, it is important to consider 
also other inflation models which naturally realize a low reheating temperature. The new 
inflation is a well-known candidate for such an inflation. In this subsection, therefore, 



19 Although gravitinos are also produced in the reheating process of the inflaton tp decay, they are 
diluted by the subsequent entropy production of the E decays and become negligible. 
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Figure 2.8: The contour lines of the reheating temperature Tr in the hybrid inflation 
model without the B — L symmetry. The contour lines are shown by the dotted lines and 
corresponding values of Tr are represented in unit of GeV. 
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Figure 2.9: The contour lines of the lepton asymmetry n^/s in the hybrid inflation model 
without the B — L symmetry for the case M x ~ 10 12 GeV (a+d—1). The contour lines are 
shown by the thin solid lines and corresponding values of n^/s are represented. Here we 
have assumed the branching ratio B r (E — > N\N\) is B r — 1 for the estimation of n^js. 
We also show the contour lines of the reheating temperature by the dotted lines and 
corresponding values of T R in units of GeV. The upper bound on k from the requirement 
(PN e < M G and the lower bound on k from m v = m s > 2M 1 are both shown by the thick 
solid lines. Aa 




we investigate the leptogenesis in the new inflation. 

In order to make the discussion concrete, we propose a SUSY new inflation model 
which has the following superpotential W and Kahler potential K: 

W = X{v 2 -g^), (2.77) 

G 



K = |$p + |X| 2 + -^|$| 4 + -^|$| 2 |X| 2 + -^|X| 4 + ---, (2.78) 



where $ and X denote supermultiplets, v is the energy scale of the inflation, g, K\, K2 
and k 3 are constants of order unity, and the ellipsis denotes higher order terms. The 
superpotential in Eq. fl2.77| ) is naturally obtained, for example, by imposing U(1)r X Z„ 
symmetry. Hereafter, we will take g and v 2 real and positive by using the phase rotations 
of X and $. 

The SUSY vacuum of the scalar potential is obtained from Eqs. fl2.57[ ), ( |2.77| ) and 
( |2.78| ), which is given by 

(2 \ V" 

Here and hereafter, we use the same symbols $ and X for the scalar components of 
corresponding supermultiplets. The scalar potential for |$|, \X\ -C Mq is given by 



, , ; , |$| \ 2 2 $ n + $* rt 2 |$ 



2n 



• (2-80) 

Here, we have neglected irrelevant higher order terms. As for the X field, it is found 
that if k 3 < —1/3, X receives a positive mass squared which is larger than H? n{ , where 
Hmi = v 2 1 (V3Mg) is the Hubble parameter during the new inflation. Thus X settles 
down at X = for k% < —1/3. Hereafter, we assume K3 < —1/3 and set X = 0. 

For g > and k = K2 — 1 > 0, we can identify the inflaton field <p with the real part 
of the scalar field $. The inflaton potential near the origin (ip ~ 0) is obtained as,[^] 

4 k 4 f (f \ 2 g 2 ip n 



V(<p)~v*- ~v* -r^ -tt^t^^W- (2-81) 

vv ^ ; 2 \MqJ 2 n / 2 ~ l M2T 2 y ' 

The slow-roll conditions for the inflation in Eq. (|2.63|) are satisfied when 

<£;<!, (2.82) 



l/(n-2) 

< = v^Mg 



gn(n - 1) M ( 



( - 2 



(2.83) 



2U 



The new inflation model proposed in Ref. mm offers a scalar potential for the inflaton field similar to 



Eq. ( |2.8l[) . Thus, our discussion in this subsection can also be applied to that model. 
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and the new inflation takes place when the inflaton (p rolls down along the potential in 
Eq. ( E3ID from (p ~ to <p f .^\ 



Now let us discuss the decay of the inflaton. We assume that the inflaton decays 
through nonrenormalizable interactions in the Kahler potential such as 

K = Y J c'm 2 \^\\ (2-84) 

i 

where ipi denote supermultiplets for SUSY standard-model particles including the right- 
handed neutrinos, and are coupling constants of order unity. (Here, a Kahler potential 
K = (1 / Mg)&* H u Hd is forbidden by suitable charge assignment of the Z n symmetry 
for the Higgs supermultiplets H u and Ha-) With the interactions in the above Kahler 
potential, the decay rate is estimated as 



2 



T ^^ c '\wJ ml ' (2 ' 85) 



where C = J2i c'- 1 and we take C — 1 in the following discussion. Here, the inflaton mass 
rritp in the true vacuum Eq. fl2.79| ) is given by 



llnM o(^r) • (2-86) 



/ v \ 2-2/n 

m v ~ ng^'"M G ^— ■ 

Then the reheating temperature Tr is given by Tr ~ (ir 2 g* /90) -1 / 4 JT '^Mq with the decay 
rate in Eq. ( ggg ). 

The scale v of the new inflation is determined from (i) the e-fold number N e of the 
present horizon, (ii) the amplitude and (iii) the spectrum index n s of the primordial 
density fluctuations Sp/p. First, the number of e-foldings N e should satisfy Eq. ( |2.73|) to 
explain the present horizon scale. In the present new inflation model, N e is given by 

N e = I d^Twr 1 — 



^ r M G \V'(<p)\ 



d<P m- + L d V — . (2.87) 

VNe kv A — J <p n 9 2 f_ 



where 



k v 



2 



l/(n-2) 



«=^ 2M °{n g Ml) • (2 ' 88) 



21 For a successful new inflation, the initial value of the inflaton field <p should be taken near the local 
maximum of the potential (p ~ j85fl . Although this seems unnatural, there have been found some 
dynamical mechanisms p6[ to solve this initial value problem. 
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Figure 2.10: A Schematic behavior of the new inflation potential for ip -C Mq. 



and ifN e is the value of the p field when the present universe crossed the horizon. Here 
we have assumed that (pN e < (p (see Fig. [2.10|) . This condition corresponds to 

k>l/[N e (n-2)}. (2.89) 



From Eqs. ( |2.83| ) , ( |2.87|) , and ( |2.88j ) we obtain the field value of p^e'- 



V2M G — 



k v 



2 \ V(«-2) 



ng Mf 



exp 



G , 



-k N e + 



nk — 1 



n - 2)k(l - k) 



^2.90) 



This (fN e should satisfy the COBE normalization in Eq. ( |2.75| ): 

5.3 x 10" 4 . 



kM G ip Ne 



(2.91) 



On the other hand, the COBE observations show the spectrum index n s as n s 



1.0 ± 0.2 [82], while in the present new inflation model it is given by [84] 



n, ~ 1 — 2k . 



[2.92) 



Therefore, we take 0.01 ^$ k ^0.1 in the following analysis assuming the coupling k not 



be extremely small. (Notice that this justifies the assumption <p Ne < <p [see Eq. (|2.89| )1.) 

From the above relations, we calculate the inflation scale v for given n, g and k, which 
are shown in Fig. 2TTT . The vacuum expectation value ($), the inflaton mass and the 
reheating temperature Tr are found in Fig. 2.12 , Fig. 2.13 , and Fig. 2. 14] , respectively. 
From Fig. |2.14| , we find that the new inflation model with the power indices n = 4, 5 and 
6 naturally offer the low reheating temperature Tr ^ 10 s GeV. For the cases n — 7, 8, we 
obtain T R < 10 8 GeV for the region k ~ 0.05-0.1 and 0.07-0.1, respectively. 
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Figure 2.11: The scale v of the new inflation for g — 1. We take the index n as n = 4, 
5, 6, 7 and 8 from the bottom to the top. 
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Figure 2.12: The vacuum expectation value ($) of the new inflation for g — 1. We take 
the index n as n = 4, 5, 6, 7 and 8 from the bottom to the top. 
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Figure 2.13: The inflation mass m$ of the new inflation for g — 1. We take the index n 
as n = 4, 5, 6, 7 and 8 from the bottom to the top. 




k 



Figure 2.14: The reheating temperature Tr of the new inflation for g — 1. We take the 
index n as n = 4, 5, 6, 7 and 8 from the bottom to the top. 
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Now we turn to discuss the leptogenesis in this new inflation model. The ratio of 
the lepton- number density to the entropy density s is again given by Eq. (|2.43|) . 
Notice that the new inflation model with n = 4 gives such a low reheating temperature 
as Tr ~ 1 GeV-10 GeV (see Fig. j2.14| ) that the required amount of lepton asymmetry 
cannot be generated, and hence we discard this case. 

The obtained lepton asymmetry for the cases n = 5, 6, 7 and 8 are shown in Fig. [2.15| 
and Fig. [Tl| by taking M 1 ~ 3 x 10 11 GeV (a + d = 1) and M l ~ 10 9 GeV (a + d = 2), 
respectively. Here we have taken relatively smaller values of Mi [see Eq. (|2.50|) ] to obtain 
a wider allowed region, taking account of 0(1) ambiguities in the FN model. 

Interesting results are given for M\ ~ 3 x 10 11 GeV. (Here, we find that the condition 
> 2Mi excludes the regions k < 1.2 x 10~ 2 and k > 4.8 x 10~ 2 for the case n = 5.) 
From Fig. 2.14 and Fig. 2.15 , we see that the sufficient lepton asymmetry can be generated 
for n = 5, 6, 7 and 8 with the low enough reheating temperature of Tr ~ 10 6 -10 8 GeV. 
In particular, for n = 5, the required lepton asymmetry ni/s ~ — 10~ 10 is obtained 
for B r ~ 1 with low reheating temperature Tr ~ 10 6 GeV. Thus, we can avoid the 
overproduction of gravitinos for a wide range of gravitino mass 1713/2 — 10 MeV-10 TeV 
in this case. 



On the other hand, when M\ 
leptogenesis does not work well. 



10 9 GeV (a + d = 2), we find from Fig. |JJ that the 
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Figure 2.15: The lepton asymmetry produced via the decays of Ni in the new inflation. 
The index n is taken as n = 5,6,7 and 8 from the bottom to the top. We take Mi ~ 
3 x 10 11 GeV (a + d= 1). For n = 5, the regions k < 1.2 x 10~ 2 and k > 4.8 x 1(T 2 are 
excluded since m v < 2Mi. 
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Figure 2.16: The same as Fig. fTT5] but for M x ~ 10 9 GeV (o + d = 2). 
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2.3.5 Topological inflation 



In this subsection, we discuss the leptogenesis in a topological inflation. If the vacuum 
expectation value of the inflaton exceeds the gravitational scale, the inflaton potential 
of the new inflation in Eq. ( |2.81| ) becomes nothing but the inflaton potential for the 
topological inflation pTf . 



We adopt the SUSY topological inflation model proposed in Ref. |8q] . The superpo 
tential and Kahler potential in that model are given by 

1 G 



W = v*X(l-g w ), (2.93) 



K = \X\ + l$l H — l$l H — \X\ l$l H —\X\ -\ (2 94) 



where v is the energy scale of the inflation and g, k±, k 2 and k 3 are coupling constants of 
order unity. Hereafter, we take g and v to be real and positive by using the redefinition of 
the field X and $. These potentials possess U(l)i? x Z 2 symmetry; the U(l)# charge of 
$ and X are and 2, respectively, and $ (X) is odd (even) under the Z 2 . This discrete 
symmetry is crucial to allow a topological defect (in this case, a domain wall) and hence 
to realize the topological inflation [J87[ . 

From Eqs. ( ^.57[ ), ( |2.93p and ( |2.94| ) we find a SUSY- invariant vacuum: 



(X) = 0, ($) = i = M G . (2.95) 

where the scalar components of the supermultiplets are denoted by the same letters as 
the corresponding supermultiplets. The potential for the region \X\, |$| Mq is written 
approximately as 



V ~ v 



4 



1-9 



$2 2 



Ml 



+ (2.96) 



Hereafter, we set X = assuming K3 < — 1/3 as we have done in the new inflation. 

For g > and k% < 1, the inflaton field (p is identified with the real part of $ and the 
potential around the origin is given by 



V{<p) ~ 1 - , (2.97) 




where A; = 2(7 + /t 2 — 1 . 

A topological inflation takes place if the vacuum expectation value of $ is of order of 
the gravitational scale Mq P7| . It was shown in Ref. that, in order for a topological 
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defect (in our case, a domain wall) to innate, its vacuum expectation value should be 
larger than the following critical value. 



($) > ~ 1.7M G . (2.98) 

The slow- roll conditions in Eq. (|2.63|) are satisfied when < k < 1 and <p?S<pf, where pj 
is expected to be of order of the gravitational scale. Hereafter, we take <pf = Mq. 

As for the decay of the inflaton ip, we assume that p decays through the interactions 
in the Kahler potential as shown in Eq. ( |2.84j ), and hence the reheating temperature Tr 
of the topological inflation is again given by T R ~ (iT 2 g*/9Q)~ 1 / A ^Jr ip MG with the decay 
rate in Eq. fl2.85| ). Here, the inflaton mass in the true vacuum in Eq. ( |2.95D is estimated 

as 

^v 2 2v 2 , . 

m ^ 2 ^w G = W)- (2 ' 99) 

As in the other inflation models, the scale of this topological inflation is determined by 
(i) the e-fold number N e , (ii) the density fluctuations 5p/p and (iii) the spectrum index 
n s . The e-fold number in the present model is estimated as 

N e ~ / d^-— — ~ -In [*L , 2.100 

J Vf M G W k \<p N J 

which should satisfy the present Hubble horizon scale in Eq. fl2.73| ), while the COBE 
normalization in Eq. ( |2.75| ) requires 



v 2 



kM G p Ne 



5.3 xlO" 4 . (2.101) 



The spectrum index n s is again given by n s ~ 1 — 2k, and hence we take 0.01 ^5 k^z 0.1 
as in the new inflation model. 

From the above relations, we can calculate the inflation scale v and the mass m^, and 



hence the reheating temperature Tr. We show the results in Fig. |2.17| . Here, we have 
taken ($) = — 1.7M G . A crucial deference between this model and the previous 

new inflation model is the reheating temperature. In this topological model, as shown in 
Fig. |2.17| , relatively high reheating temperatures are obtained compared with those in the 



new inflation. This is because the large value of the vacuum expectation value ($). (See 
Eq. (pp.) 

Finally, we show the obtained lepton asymmetry ul/s in Fig. |2.18| . Here, we have 
taken M 1 ~ 3 x 10 11 GeV (a + d = 1) and M l ~ 10 9 GeV (a + d = 2) as in the 
new inflation model discussed in the previous subsection. We can see that for the both 
cases enough lepton asymmetry n^/s > 10~ 10 is obtained for B r ~ 1. On the other 
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k 



Figure 2.17: The inflation mass and the reheating temperature Tr of the topological 
inflation model. We take (<&) = {<&) c — \.1Mq. The solid and dot-dashed lines represent 
m v and T R , respectively. 



<^)>=<(|)> 




k 



Figure 2.18: The lepton asymmetries, ul/s, in the topological inflation model for the 
cases of M x ~ 3 x 10 11 GeV (a + d = 1, the solid line) and M x ~ 10 9 GeV (a + d = 2, the 
dot-dashed line). 
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hand, the requirement of a low enough reheating temperature (Tr ^ 10 8 GeV) is satisfied 
only in a small range of the coupling 0.092 ^ k ^ 0.1. (The lowest reheating temperature 
Tr = 6.3 x 10 7 GeV is obtained for k = 0.1.) For k ~ 0.1, the spectrum index n s deviates 
from the scale invariant one as n« ~ 0.8, and it will be testable in current and future 



satellite experiments on anisotropics of the cosmic microwave background radiation p0 |. 

One final comment: in the case of Mi ~ 3 x 10 11 GeV, we may accommodate a 
late-time entropy production of order 10 2 -10 3 since the produced lepton asymmetry is as 
large as n^/s ~ 10 -7 as shown in Fig. |2.18| . If this is the case, the energy density of the 



gravitinos is diluted by the factor 10 -10 , which allows the region of Tr ~ 10 GeV in 



Fig. 2.17 



2.3.6 Production of right-handed neutrinos at preheating 

Here, we briefly mention the production of right-handed neutrinos at preheating epoch, 
which was investigated in Ref. J74] . 



So far, we have assumed that the reheating process after the inflation occurs in such a 
way that energy density of the zero mode of the inflaton decays perturbatively into light 
particles. However, the heavy right-handed neutrinos might be predominantly produced 
nonperturbatively, through the parametric resonance PTl B3]. 



In Ref. 0], production of heavy right-handed neutrinos via preheating of fermion |92j 
was discussed, adopting a chaotic inflation with the inflaton potential V((p) = {\/2)m ip ip 2 . 
im^p is the inflaton mass.) The point is that, if there is a Yukawa coupling between the 
inflaton (p and heavy fermion X 

C = g(pXX, (2.102) 

then the effective mass mxe& of the fermion X depends on time t through the oscillating 
inflaton <£>(£) during the inflaton-oscillating period: 

m XeS (t) =m x + 9<p(t). (2.103) 

Therefore, the effective mass mx e s{t) oscillates with time, and can even be zero if the 
amplitude of the oscillating inflaton is large enough, which is very crucial to produce very 
massive fermions heavier than the inflaton itself. 

The production of the X fermion crucially depends on the parameter 

8S 9*2)1. (2.104) 
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where ip(0) denotes the initial amplitude of the inflaton oscillation. An efficient production 
of very massive fermions requires a large value of q ||74|| , and hence large enough initial 
amplitude <p(0). In the case of chaotic inflation, ip(0) — Mq and m v ~ 10 13 GeV can lead 
to q ~ g x £>(10 10 ) > 1. 

The maximum mass of the produced fermion is given by J74 



(™*Lax- 0.25^(0). (2.105) 

If mx ~ ( m x) max , the effective mass mx e s(t) can vanish more than once [see Eq. ( |2.103| )1P| 
and the production of X through parametric resonance occurs efficiently, whereas if oth- 
erwise the production is much less effective. The production continues during vanishing 
effective mass can be achieved. Namely, the production stops when the amplitude of the 
oscillating inflaton tpo becomes below the following critical value: 

Wc* = ~ g m x • (2-106) 

The energy density of the inflaton at this time is thus given by p v = (l/2)m^ ((po) ait - On 
the other hand, the energy density of the produced fermion X at this frozen-out time is 
given by [73 1 

p x ~-^m>;. (2.107) 

After that, i.e., for ip Q < (v 5 o) cr it> the ratio of the energy density of the X to the inflaton 
(p, px/Pip, is fixed. From the above equations, this ratio is simply given by 

— (2-108) 

Now let us identify the X field to the right-handed neutrino N±. From above equation, 
the resultant lepton asymmetry after the reheating process of the inflation completes is 
given by 

^ (2.109) 

Thus, the baryon asymmetry %/s ~ 0.35 x n^j s is given by 
n B . , n 2 w 1 n-10 ( \ . . ( m u3 



22 Their numerical calculation shows ]74| that the minimum of the inflaton amplitude <^q ^ during the 
first oscillation is given by (p^ ' ~ —0.25(^(0). 
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where we have used Eq. ( |2.15| ). 

Before closing this section, let us see what happens if we apply the above mechanism 



to the first one of the SUSY hybrid inflation models discussed in Sec. |2.3.3.a| , where the 
inflaton has a Yukawa coupling to the right-handed neutrino. In this case, the parameter 
q is estimated as 

q ~ 

fnt- 

1 / O T\/T \ 2 

(2.111) 

Therefore, in the region we have studied, where the perturbative decay is kinematically 
allowed (m^ > 2Mi), preheating takes place only through a narrow resonance, and hence 
only small fraction of the inflaton oscillating energy can be converted to heavy right- 
handed neutrinos via preheating. For the case of Mi ~ 10 12 GeV (a + d — 1), in the 
small coupling region k < 10" 4 and A < 1(T 4 where ^(l/2)Mi [see Fig. |2.3|| , preheating 
mechanism might work effectively, although it requires another detailed analysis. 
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2.4 Leptogenesis from TV-dominant universe 

In the previous sections, we have discussed the leptogenesis mechanisms by the decay 
of right-handed (s) neutrino, which is produced as a particle. However, there is a new 
possibility once we introduce the SUSY, that is, the condensation of the scalar component 
of the right-handed neutrino, iV p5|.0 The lepton asymmetry is produced in the same way 
as the previous mechanisms, by the asymmetric decay of the right-handed "s" neutrino into 
leptons and anti-leptons. However, the thermal history of the universe in this scenario is 
quite different from the others. In this section^ we investigate this scenario in detail. In 
particular, we mainly discuss the case in which the coherent oscillation of the right-handed 
sneutrino iV dominates the energy density of the universe before it decays. Remarkably, 
the amount of produced baryon asymmetry is determined almost only by the decay rate of 
the right-handed neutrino, whatever happened before the coherent oscillation dominates 
the universe. (Thus, the baryon asymmetry is independent of the details of the inflation !) 

Furthermore, as a big bonus, the abundance of the gravitinos which are created in the 
reheating phase after inflation are diluted by the entropy production due to the decay of 
the coherent right-handed sneutrino, so that the cosmological gravitino problems can be 
avoided even when the reheating temperature Tr of the inflation is higher than 10 11 GeV, 
in a wide range of the gravitino mass 7713/2 ~ 10 MeV-10 TeV. In particular, this dilution 
of the gravitinos has great advantages in the gauge-mediated SUSY breaking (GMSB) 



models [45 



We will also comment on the cosmic density perturbation caused by the quantum 
fluctuation of the right-handed sneutrino, which may be detected in the future. 

2.4.1 Lepton asymmetry from iV-dominant universe 

The main character in the present scenario is the right-handed sneutrino, N. During 
inflation, it can acquire a large amplitude [^5], |93[ if the Hubble expansion rate of the 



inflation ifi n f is larger than the mass of the N. (The mechanism which makes a scalar 
field develop a large initial amplitude during inflation will be explained in Sec. |3.2.2| , 



in the context of the Affleck-Dime condensation.) Let us assume that there exists (at 
least) one right-handed neutrino with a mass lighter than H in {, and that it develops a 
large expectation value during the inflation. Hereafter, we focus on the lightest right- 
handed sneutrino Ni for simplicity. (Possible contributions from the heavier right-handed 



23 See also Refs. |2J, [)3|, where the right-handed sneutrino is considered as an inflaton of the chaotic 
inflation. 

24 This section is based on the work in a collaboration with H. Murayama and T. Yanagida f2q|. 
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sneutrinos N 2 and N 3 will be discussed in Sec. [2.4.4. b| .) It is assumed here that the 
potential for the right-handed neutrino is given simply by the mass term 

V = Ml\N x \ 2 . (2.112) 

We do not settle the initial amplitude of the Ni at this stage, which will be discussed 
later. Because of the large amplitude of the Aq, the lepton L and Higgs H u doublets 
receive large masses, and the expectation values of their scalar components (^L^ and (H u ) 
vanish. 

The evolution of the Aq is as follows. After the end of the inflation, the Hubble 
parameter H decreases with cosmic time t as H oc t -1 , and Aq begins to oscillate around 
the origin when H becomes smaller than the mass of the right-handed sneutrino Mi. 
(H ~ Mi.) Then, the coherent oscillation eventually decays when the Hubble expansion 
rate H becomes comparable to the decay rate T~ of Aq. (H ~ T^.) As explained in 
Sec. |2.1| , the asymmetric decay of Aq into leptons and anti-leptons produces a net lepton 
asymmetry. The generated lepton number density is given by 

n L = ei M 1 |iV u | 2 , (2.113) 
where \Ni d \ is the amplitude of the oscillation when it decays, and e\ is given by Eq. ( |2.15|) . 



The fate of the generated lepton asymmetry depends on whether or not the coherent 
oscillation of Aq dominates the energy density of the universe before it decays |[25|| . Here, 
we focus on the leptogenesis scenario from the universe dominated by Aq. (We will 
mention the case where Aq does not dominate the universe in Sec. |2.4.4.d .) As we shall 



show soon, once the Aq dominant universe is realized, the present baryon asymmetry 
is determined only by the properties of the right-handed neutrino, whatever happened 
before the Aq dominates the universe. We first derive the amount of the generated lepton 
asymmetry just assuming that the Aq dominates the universe, and after that we will 
discuss the necessary conditions of the present scenario. 

Once Aq dominates the universe before it decays, the universe is reheated again at 
H ~ by the decay of Aq. The energy density of the resulting radiation, with a 
temperature T Nl , is given by the following relation; 

^2* = Ml\N ld f 

= 3M£l%, (2.114) 



while the entropy density is given by [see Eq. ([B.17Q ] 

2vr 2 
S= l5 



9*T 3 Nl . (2.115) 
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From the Eqs. ( |2.15| ), ( p . 1 1 3| ) , ( [2.114| ) and ( |2.115| ), the ratio of the lepton number density 



to the entropy density is given by the following very simple form; 

ul_ 3 Tivi 
s ~ 4 6l Mi 

Here, it is assumed that the decay of the N\ occurs in an out-of-equilibrium way, namely, 
TjVi < M\, so that the produced lepton- number asymmetry not be washed out by lepton- 
number violating interactions mediated by N±. As explained in Sec. |1.4| , this lepton 
asymmetry is partially converted into the baryon asymmetry through the sphaleron effect: 

— ~ 0.35 x — 

Therefore, as stressed previously, the baryon asymmetry in the present universe is indeed 
determined only by the decay temperature of the right-handed sneutrino (and the 
effective CP violating phase S e s). Thus it is independent of unknown parameters of the 
inflation such as the reheating temperature Tr. Assuming the effective CP violating phase 
^efr (< 1) to be not too small, the empirical baryon asymmetry %/s - (0.4-1) x 10~ 10 is 
obtained by taking 



■Ni 



10 6 - 10 7 GeV. (2.118) 



2.4.2 Necessary conditions 

So far, we have required the following two conditions; 

(i) Ni dominates the universe before it decays. 

(ii) Ni decays in an out-of-equilibrium way. 

Let us first consider the easier one, the second condition (ii). By taking the in 
Eq. Q2.118D , the condition of the out-of-equilibrium decay is simply given byQ 

Mi > T Nl ~ 10 6 - 10 7 GeV . (2.119) 



25 The condition Mi > Tjf ± is formally equivalent to 



if we take the usual relation 3MqH 2 — (ir 2 /30)g*T 4 . However, this interpretation is not appropriate in 
the present case, since the universe is not radiation-dominant for H >Tn 1 . In any case, the condition 
Eq. ( gjUj| is necessary to avoid the thermal production of N% after the decay of the coherent Ni. 
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Notice that this condition gives rise to a constraint on the Yukawa couplings hi a , since 
the temperature is determined by the decay rate of the Ni, and hence is related to 
the mass and couplings of Ni. Actually, from Eq. fl2.114|) and = (l/47r) J2a \ hi a \ 2 Mi, 
the following relation is derived: 

, rp . 1/2 ,rp . 1/2 

>-6 / \ MATi\ 



^M^xiO-^j . (2,20) 

Hence, it is found from Eq. ( ETTT51 ) that Yukawa couplings h\ a of N\ should be smaller than 
O(l0~ 5 ). (Possible explanations for the smallness of h la will be discussed in Sec. 2.4.4. a| .) 



Now let us turn to discuss the first condition (i). In order to discuss whether or not 
Ni dominates the universe, it is necessary to consider the history of the universe before 
it decays. Here, we assume that the potential of the Ni is "flat" up to the Planck scale, 
namely, the potential is just given by the mass term Mf\Ni\ 2 up to the Planck scale. 
(This may not be the case when the masses of the right-handed neutrinos are induced by 
a breaking of an additional gauge symmetry. We will discuss such a case in Sec. |2.4.4.d| .) 

Assuming the flatness of the iVi's potential up to the Planck scale (i.e., only the mass 
term), the initial amplitude of the oscillation is naturally given by \Nn\ ~ Mq, since 
above the Planck scale the scalar potential is expected to be exponentially lifted by the 
supergravity effects [see Eq. ( 2.571) ].^ Then the energy density of Ni when it starts the 



coherent oscillation is given by pjvi = Mf\Nn\ 2 ~ M\M G . 

The rest of the total energy density of the universe at H = Mi is dominated by (I) 
the radiation or (II) the oscillating inflaton \, depending on the Mi and the decay rate 



of the inflaton T x . (See Fig. |2~T9 



If T x > Mi [case (I)], the inflaton decay has already completed before H = Mi, and 
the energy density of the radiation at H = Mi is given by p md = (p to tai — PjvJ — 2Mf M G . 
In this case, the oscillating N\ dominates the universe soon after it starts the oscillation 
and hence before its decay. 

On the other hand, if Mi > T x [case (II)], the reheating process of the inflation has 
not completed yet at H = M%, and the inflaton \ is still oscillating around its minimum, 
whose energy density is given by p x ~ 2M\M 2 G . The ratio of the energy density of N\ to 
that of the inflaton, Pn 1 /Px — 1/2, takes a constant value until either of these oscillations 
decays. Because the energy density of the radiation p ra d resulting from the inflaton decay 
is diluted faster than p^ 1) the oscillating N\ can dominate the universe if its decay rate 
Y Nl is slow enough compared with that of the inflaton T x ; T x ^> r^-EO 



26 Even though it is possible that Ni has a larger initial amplitude \Nn\ > Mq (see, e.g., Ref |24], |93|), 
it depends on the scalar potential beyond the Planck scale, so that we do not discuss this possibility here. 



27 More precisely, pjvi > Prad occurs for < (l/2\/3)r x in case (II). 
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Figure 2.19: Schematic behaviors of the energy densities. R denotes the scale factor of 
the expanding universe. It is found that the coherent oscillation Ni can dominate the 
energy density of the universe if the condition r x 3> r^i (Tr TjvJ is satisfied. 

Therefore, the condition for Ni to dominate the universe is just given by T x ^> T^. 
In terms of the reheating temperature Tr, it is 

T R > T Nl ~ 10 6 - 10 7 GeV , (2.121) 

which is easily satisfied in various SUSY inflation models. (Recall that the SUSY inflation 
models in Sec. [T3] naturally lead to Tr 3> 10 6 -10 7 GeV.) Thus, the present leptogenesis 
scenario from N% dominated early universe is almost automatic as long as the right-handed 



neutrino has suitable mass and couplings given in Eqs. (2.119) and (2.12C). 



2.4.3 Dilution of the gravitinos 

A very attractive feature in the present "A?i-dominant universe scenario" is the dilution 
of the gravitino due to the entropy production caused by the Ni decay. We discuss this 
point in this subsection. 

Let us first calculate the dilution factor of the entropy production caused by the N±s 
decay. It depends on the two cases discussed in the previous subsection, M 1 < T x [case (I)] 
and Mi > T x [case (II)]. We first consider the latter case, M x > T x [see Fig. |TT| (II)]. At 
the time of inflaton decay (H = T x ), the ratio of the energy density of Ni to the entropy 
density is given by 



PNi 

s 



before JVi decay Prad 



PNl ~ — Tr , (2.122) 
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where we have used p^/Prad — 1/2 and p ra ,d/s — (3/4)Tr at H — Y x . This ratio takes a 
constant value until the iV\ decays, since p^ t and s are diluted at the same rate R~ 3 as 
the universe expands. Then at the time of the Ni's decay, the energy density of the N\ is 
converted to the radiation. The ratio of the energy density to the entropy density in the 
new radiation is given by 



Prad 



= jT Nl . (2.123) 

after Ni decay '« 



By equating the p^ (just before the Ni's decay) in Eq. fl2.122|) with the p ra d (just after 



the Nis decay) in Eq. ( |2.123| ), we obtain the dilution factor of the entropy production: 

A = ^^ = -^. (2.124) 

■^before 2 

Next, we consider the case T x > Mi [case (I)]. In this case, the ratio p^/s takes a 
constant value after the iVi starts its oscillation, i.e., for H ^ M\. This ratio is evaluated 
in the same way as Eq. fl2.122|) : 



s 



_ = \t Mi (2.125) 

before TVi decay o 



where Tm 1 is the temperature at H = Mi in the case (I), which is given by 

/ /If \ !/ 2 

r * 3 7 x 10 " (lo^v) GeV ' (2126) 



Then the dilution factor at H = is obtained in the same way as in Eq. ( [2.124D , which 
results in A = T Ml /(2T Nl ). 

To summarize, the dilution factor A is given by 

Tr 

— — (for T R <T Ml ) 
A^<( ^ , (2.127) 



2T, 



(for T R >T Ml ] 



Ni 



where we have used the fact that T x < Mi (T x > Mi) corresponds to T R < T Ml {T R > 
T Ml )- 

As discussed in Sec. [T7|, the abundance of the gravitino produced at the reheating 
epoch after the inflation is proportional to the reheating temperature of the inflation 
T R . Therefore, there are usually severe upper bounds on T R depending on the gravitino 
mass. However, in the present scenario, the number density of the gravitino is diluted 
thanks to the entropy production caused by the AVs decay, which was discussed just 
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above. Namely, the constraint coming from the gravitino problems should apply, not to 
the reheating temperature T R , but to an effective temperature given by 



2T, 



1 j . . . 

TR e ff = —T R ~ { n rr ( T R 




i / 

7 



2 x 10 7 GeV (for T R < T Ml ) 

2 x 10 6 - 2 x 10 7 GeV x ( (for T R > T Ml ) ,(2 ' 128) 

\Tm 1 J 

which is much below the original reheating temperature T R . Actually, for T R < Tm x , 
the gravitino abundance which was originally produced at the reheating after primordial 
inflation is diluted and becomes comparable to the abundance which is produced at the 
"reheating" after AVs decay. 

In the case of unstable gravitino, Eq. ( |2.128j ) means that the cosmological gravitino 
problem coming from the big-bang nucleosynthesis bounds can be avoided in a wide range 
of the gravitino mass m 3 / 2 ~ 100 GeV-10 TeV, even if the reheating temperature T R of 
the inflation is much higher than 10 11 GeV. (Recall that the Tm 1 is given by Eq. ( |2.126| ).) 
We should stress that the fact that such high reheating temperature is allowed makes it 
very easy to construct SUSY inflation models (see Sec. |2.3| ). 

In the case of stable gravitino as well, we see from Eq. (|2.128| ) that the overclosure 
problem can be avoided in a wide range of the gravitino mass. Actually, in our scenario, 
the relic abundance of the gravitino is obtained by dividing the original abundance in 
Eq. ( |L23D by the dilution factor A in Eq. ( F~l27p : 

1 



a 3/2 h 2 



with Ni decay ^ 



— ^3/2 h 



without Ni decay 

m a Yf m 3/2 \~V TrcS ^ 



Therefore, the overclosure problem can be avoided almost independently of the reheating 
temperature T R , and a reheating temperature even higher than 10 11 GeV is possible for 
m 3/2 ~ 10 MeV. Moreover, it is found from this equation that the present energy density 
of the gravitino is independent of the reheating temperature T R , in a very wide range of 
Tjv x <C T R < Tjvfi . Thus, we can predict the gravitino mass by requiring that the gravitino 
is the dominant component of the dark matter^ 

on ~\ A \T ( m 9 V /^matter/* 2 \~ { T RcS \ 

ms/2 K 20MeVx(— -BIS" [ww) 



2S 



Notice that with such a small gravitino mass as in Eq. (2.130), the problem of the decay of the next 



to-lightest SUSY particle into LSP gravitino after the big-bang nucleosynthesis is also avoided 



See Sec. 1.5 
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(40 - 400) MeV x 



rrih 



'a 



h 2 



for T Nl < T R < Tj 

/» ~ 0.7 ra. 



Mi- 



.1 TeVy V 0.15 
Here, we have taken the present matter density f2 matter 



(2.130) 
0.3 and 



Such a light (stable) gravitino as in Eq. ( |2.130|) naturally exists in gauge-mediated 

The GMSB mechanism has been regarded 



SUSY breaking (GMSB) models |S 



very attractive candidate for the SUSY breaking, since it suppresses quite naturally the 
flavor changing processes, which are inherent problems in the SUSY standard model. It 
is remarkable that the gravitino with mass 7713/2 — 40-400 MeV becomes naturally the 
dominant component of the dark matter in the present scenario, independently of the 
reheating temperature T R for T Nl <C T R < T Ml . This prediction comes from the fact that 
the present energy density of the gravitino is determined by the effective temperature 
TRett = 2T/V! (for Tr < TmJ, while the decay temperature of the right-handed neutrino 
Tjvj is fixed by the baryon asymmetry in the present universe [see Eq. (|2.118|) 1. 



2.4.4 Remarks 

We give several comments on the present scenario. 
2. 4. 4. a neutrino masses m ui 

The first one is about the neutrino mass m v . Let us first rewrite the neutrino masses in 
terms of the three contributions from Ni, N% and iV^: 



ni 



„from Ni 1 „,from No i „,from N3 



where 



rn 



from N; 



(Huf 
Mi 



(2.131) 



(2.132) 



and the denote the Yukawa couplings of the neutrinos in the basis where mass matrices 
of heavy and light neutrinos are diagonal, which are defined in Eq. ( |2.12| ). Then we see 
that the contribution from Ni satisfies the following relation:^] 



rn 



.from Nj 



\hijhij 



lkl 



Mi 

2 l \ H uY 

Mi 



29 r 



The last expression in Eq. (2.133) is nothing but the mi parameter in Eq. (2.2S 
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,2<#u> o 1 n-4 „v /^iVi 



= E l^| 2A T7^ - § X !0~ 4 eV ^ , (2.133) 

Here, we have used the relation in Eq. (|2.120| ). Therefore, it is understood that the mass 
scale of the neutrinos suggested from the atmospheric |T(J and solar ill] (+KamLAND [fT2]| ) 



neutrino oscillations, ~ (9(10 3 )-(9(10 x ) eV, should be induced from the heavier 
right-handed neutrinos, N 2 and iV 3 . 

The relative hierarchy between the mass and couplings of Ni and those of the N 2 
and iV 3 might be naturally explained by a broken flavor symmetry. For example, a 
broken discrete Zq symmetry with a breaking parameter e ~ 1/17 and charges 

Q(L e , L M , L T ) = (a + 1, a, a) and Q(Ni, N 2 , N 3 ) — (3 + d, c, 6) gives rise to the following 
superpotential; 

where ^ and /ij Q are 0(1) couplings. The above charge assignments for lepton doublets 
naturally lead to the realistic neutrino mass matrix including the maximal mixing for 
the atmospheric neutrino oscillation f76|| . (See Sec. [2.3. 2| .) The overall mass scale of the 
right-handed neutrino Mq is determined by m v3 ~ e 2a (H u ) 2 /Mq. By taking a + d — 2, 



this model gives M 1 ~ e 2d M ~ 7 x 10 9 GeV, ^JJ2 a \hia\ 2 ~ £ 5 ~ 7 x 1(T 7 , and hence 
TV, ~ 1 x 10 7 GeV.0 

Recently, it has been shown [96| that the small mass and the couplings of the N%, 



required in the present scenario, naturally arise if the right-handed neutrino is a bulk 
field in higher dimensional theory. 

2.4.4.b contributions from N 2 and N 3 

In the previous subsections, we have only considered the lightest right-handed sneutrino, 
Ni. The heavier right-handed sneutrino iV2( 3 ) can also develop a large amplitude during 
the inflation (if M 2 ( 3 ) < H in f) and it may produce lepton asymmetry in a similar way 
to the Ni. However, it is found that the decay temperatures of the N 2 and N 3 cannot 
satisfy the out-of-equilibrium condition T 2 (3) < M\, since N 2 and N 3 must explain the 
mass scales of the neutrino oscillations. (In other words, if Tjv a(3) would also satisfy the 
condition T/v 2{3) < Mi < M 2 ( 3 ) , the mass scale for the neutrino oscillations could no longer 
be explained by the contributions from N 2 and N 3 [see Eq. ( |2.133| )].) Therefore, even if 
the A^ 2 ( 3 )'s decay produces additional lepton asymmetry, it is washed out and hence it 
cannot contribute to the resultant total lepton asymmetry. 



30 The model here looks similar to the one in Sec. 3.3.1, but different from that. 
31 In terms of mi parameter, this corresponds to m~i ~ £ 6 m^ ~ O(10 -9 ) eV. 
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2.4.4.C thermal effects 



We should also comment on the effects of the thermal plasma, which might cause an early 
oscillation of the right-handed sneutrino Ni before H = Mi. (Notice that there is a dilute 
plasma with a temperature T ~ [T^MqH) 1 ^ even before the reheating process of the 
inflation completes ||28|| .) There are basically two possible thermal effects |J7|, |9*5| . (Those 
thermal effects will be also explained in detail in Sec. |3.2.3| .) 

First, when the temperature T is higher than the effective mass for L and H u , 

\ 2 \Ni\, the Ni receives an additional thermal mass SM 2 } 



i.e., T > m, 



efi 



J2a \hic 



'th 



(1/4) J2 a \ hia\ 2 T 2 from the Yukawa coupling to L and H u pTfl . Thus, the N\ field would 



start an early oscillation if the additional thermal mass becomes larger than the Hubble 
expansion rate before H = Mi. However, even if iV\ receives the thermal mass, the ratio 
of the thermal mass to the Hubble expansion rate is given by 



H 2 



0.07 x — ^ 

V Mi 

!0T Nl 



0.03 x 



Mi 



Mi 

IT 

Mi 

IT 



3/2 



T 



ft 



T 



Mi, 



for 
for 



Tr < T, 



Mi 



(2.134) 



Tr >T 



Mi 



where we have used the relation given in Eq. 00201) . Therefore, we can safely neglect 
the above thermal effect as long as Mi is a bit larger than T^. 

There is also another thermal effect |98| . If the temperature is lower than the effective 
mass for L and H„ 



l u , i.e., T < m eff = yJ2a |^ia| 2 |^i|> the evolution of the running gauge 
and/or Yukawa coupling constants f(n) which couple to them are modified below the scale 
[i = m eS . Thus, these running coupling constants depend on |iVi|, and there appears an 
additional thermal potential for Ni, 



SV{Ni) = aT 4 log 



J^2 



(2.135) 



where a is a constant of order 0(f ). However, again, it turns out that the effective 
thermal mass for Ni is less than the Hubble expansion rate; 



H 2 



aT 4 



H 2 \Ni\ 2 



0.2 x a 



0.05 x a 




for T R < T, 



Mi j 



(2.136) 



for Tr > TV 



and hence this thermal effect is also irrelevant to the present scenario. 
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2.4.4.d more general cases 



So far, we have assumed that the initial amplitude of the iVi's oscillation is \Nn\ ~ Mq- 
This can be realized when the right-handed neutrino has only the mass term up to the 
Planck scale. We have also assumed that the Nx dominates the energy density of the 
universe before its decay. Here, we consider more general cases for completeness. 

Let us assume that the masses of the right-handed neutrinos are dynamically induced 
by a spontaneously broken gauge symmetry, for example, by the U(l) b-l symmetry, 
which is an attractive candidate for such a gauge symmetry as discussed in Sec. |L~T[ We 
denote the chiral superfields whose vacuum expectation values break the XJ(1)b-l by \I/ 
and ^. (We need two fields with opposite charges in order to cancel XJ(1)b-l gauge 
anomalies.) Due to the D-term and the F-term coming from the superpotential which 
gives the right-handed neutrino masses, the scalar potential of the right-handed sneutrino 
Nx is lifted above the XI (1) b-l breaking scale (\&) [P5[ . Therefore, the initial amplitude 
of the iVi's oscillation at H ~ Mi is given by \Nxi\ ~ 

The breaking scale of the XJ(1)b-l gauge symmetry is model dependent. If it is broken 
at the Planck scale, (^f) ~ Mq, the discussion so far does not change at all.0 On the 
other hand, if (&) is below the Planck scale, the initial amplitude of the Nx's oscillation 
is reduced and its energy density when it starts the oscillation is given by 

p Nl = Ml\N u \ 2 < p total = ZMlM 2 G for H = M x . (2.137) 

Compared with the previous case (see Fig. |2.19| ), we can see that the initial energy density 
is much suppressed. 

In Fig. |2.20| , We show the behavior of the energy densities of the inflaton x, the right- 
handed sneutrino Nx and the radiation, for generic cases.0 Figure (la) and (lb) represent 
the case of T x > Mi (Tr > TmJ, (Ha) and (lib) represent the case of Mi > T x > 
(T Ml > T R > T Nl ) and (III) represents that of T Nl > T x (T Nl > T R ). Among them, 
(la) and (Ha) show the case of Nx dominant universe, corresponding to (I) and (II) 
in Fig. p,19| , respectively. Figures (lb), (lib), and (III) represent the cases in which the 
coherent oscillation of the right-handed sneutrino N does not dominate the energy density 
of the universe. 



32 In this case, we need small couplings in order to explain the intermediate right-handed neutrino mass 
scale. For example, a superpotential W — (l/2)gi 1 ^NiNi with (VP) ~ Mq and g$ ~ 10~ 4 gives the mass 
JVf 3 ^ 10 14 GeV to the heaviest right-handed neutrino. Such a small Yukawa coupling could well be a 
consequence of broken flavor symmetries. 

33 For the reduced initial amplitude, the thermal effect from the aT 4 log(|7Vi | 2 ) potential becomes larger 
than the case of lA^I ~ Ma- However, it is still irrelevant as long as the coupling a ~ C(/ 4 ) is small 
enough. [See Eq. ( p!l36| ).] 
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Figure 2.20^ Schematic behaviors of the energy densities of the inflaton x, right-handed 
sneutrino N±, and radiation for general cases. 
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Figure 2.21: The condition of the Ni dominant universe is represented in (T x , |iV"ij|) plane. 
The behaviors of the energy densities for the regions of (la), (lb), (Ha), (lib) and (III) 



are shown in the corresponding figures in Fig. 2.20 



Let us first discuss the condition of the Ni dominant universe, which was given by 
[Eq. ( |2.121| )1 for |iVi| ~ Mq. Because of the smaller initial amplitude, the condition now 
becomes 



~ /r \ 1/4 

N u \ > V3M G ^j for T x >M l 



fla) 



\N U \ > VSM G 



T 



Ni 



1/4 



for T Nl <T X <M 1 



(Ila) 



(2.138) 
(2.139) 



These conditions are summarized in Fig. |2.21| . Notice that the Ni cannot dominate the 
universe for r x < T Nl [the case (III)]. 

Next we estimate the baryon asymmetry in the above general cases. It can be calcu- 
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lated in the same manner as before, and is given by, for each case, 



(Ia)(IIa) — = 0.35x- ei 
s 4 

(lb) — = 0.35 x - ei 

s 4 

(lib) (III) — = 0.35x|ei 







Mi 




( Wu\ ^ 






Mi 


( Wu\ \ 


2 


\VsMgJ 


Mi 



0.5 x 10~ 10 
0.5 x lO" 10 
0.5 x 10~ 10 



-in ( 



VlO 6 GeV 



( \N U \ 
VlO 6 GeV 7 I V3M ( 



Tn 



( \N U \ 
10 6 GeV 7 I V3M ( 



(2.140) 



where 77 = (m^s/0.05 eV)5 e g. Notice that the baryon asymmetry which was given in 
Eq. ( |2.117| ) does not change and independent of the initial amplitude \N U \ as long as the 
Ni dominant universe is realized [(la) and (Ila)]. 

Finally, we comment on the gravitino problem. The effective temperature of the 
cosmological gravitino problem, which was given by [Eq. ( [2.128 )1, is now given by 

2 



^3/2 h 2 OC T; 



T, 



T, 



RcS 



T, 



/V3Mg\ 

1 V \N U \ ) \ T Mi 
( V3M G \ 2 

1 I \N U \ ) 



ft 



(la), 
(Ila) , 

(lb) (lib) (III) 



(2.141) 



2.4.4.e cosmic microwave background (CMB) anisotropy from N\ 

Let us mention the cosmic density perturbation coming from the fluctuation of the initial 



amplitude of the right-handed sneutrino, \N U \ } according to the recent study in Ref. |99 

In the present scenario, the CMB radiation we observe today has two origins; the 
inflaton \ an d the right-handed sneutrino iV\. Furthermore, if the Hubble expansion 
rate of the inflation is larger than the mass of the right-handed sneutrino Mi, as we 
have assumed, the initial amplitude of the N\ has a non-vanishing fluctuation during the 
inflation. Therefore, the CMB anisotropy also has two origins, the primordial fluctua- 
tions of x an d that of N\. If the N\ completely dominates the energy density of the 



universe, the perturbation becomes adiabatic |100|. On the other hand, if the radiation 



coming from N\ is negligible compared to that from x, purely baryonic isocurvature per- 
turbation is produced from iVj. Particularly interesting is between these two cases (i.e., 
the region near the border line between the two regions (la)- (Ila) and (lb)- (lib)- (III) in 
Fig. |2.21|), in which correlated mixture of the adiabatic and isocurvature fluctuation is 



generated, that can be observed in the current and future satellite experiments on the 
CMB anisotropies 
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2. A TV-dominant universe from inflaton decay 



In this appendix, we briefly discuss the case where the N± which are produced in inflaton 



decay dominates the universe before it decays. (This case was also discussed in Ref. [[74).) 

The baryon asymmetry produced by the N% decay in this case is the same as that of 
the iV-dominant universe (Sec. |2.4| ), and is given by Eq. (|2.117|) : 



n B _ . n r . . i r.-lO ( ^JVl A ( m u3 



which requires the decay temperature of the Ni to be — 10 6 -10 7 GeV for 6 c s ~ 0.1- 
1.Q 

On the other hand, the condition for the Ni to dominate the energy density of the 
universe in this case is derived as follows. The point is that the ratio n^/s, which is 
given by Eq. ( |2.42| ), takes a constant value until the time of Nis decay. Thus, the ratio 
of the energy density of the N% to that of the radiation at Aq's decay is given by 

PNl _ Mi UN Lx S 



Prad $ Prad 



M x -B 



2 m x 3Tat x 

= 2B r ^^L. (2.143) 

m x T Nl 

Thus, the dominant condition p^/prad > 1 gives a lower bound on the Mi as follows: 

TTl 

M x > ^7^T Wl , (2.144) 

which is stronger than the out-of-equilibrium condition (Mi > T^) for m x > 2B r T R . 
Then using the relation in Eq. ( |2.120| ), we obtain a constraint on the Yukawa couplings: 

£ < 5 x 10-B, (^-) ( 10 ' 3 GeV ) x (**-) . (2.145) 
V ~ \10 9 GeVj\ m x ) \10 6 GeV J V ; 

If this is satisfied, Ni produced from inflaton decay can dominate the energy density of 
the universe, and the produced baryon asymmetry is given by Eq. ([2.142 ) . 



:!4 



This point was not discussed in Ref. |74 
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Chapter 3 

Leptogenesis via LH U flat direction 



In the previous chapter, we have investigated various leptogenesis mechanisms in which 
lepton asymmetry is produced by the CP- violating, asymmetric decay of the heavy right- 
handed (s)neutrino into leptons and anti-leptons. In this chapter|] we will discuss a 
completely different leptogenesis scenario. The main character in this chapter is not a 
heavy right-handed (s)neutrino, but a flat direction scalar field including the charged 
lepton doublet L. Actually, as we will see, most of the necessary ingredients for this 
scenario exist in the minimal supersymmetric standard model (MSSM), and basically it 
does not rely on any new physics beyond supersymmetry, the existence of small Majorana 
neutrino masses and inflation. 



1 This chapter is based on the works with T. Asaka, M. Fujii and T. Yanagida JlQl], |M IIM Il03| 
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3.1 Overview 



Affleck and Dine |27j proposed a completely new mechanism to produce cosmological 



baryon asymmetry, which is inherent in the SUSY framework. In the SUSY theories, there 
are many flat directions in the scalar potential which contain scalar quark (squark) and 
scalar lepton (slepton) fields, carrying baryon and/or lepton charges. During the inflation, 
such a scalar field (which will be denoted by 0) can acquire a large vacuum expectation 
value along the flat direction, since the potential is flat. After the inflationary stage ends, 
the flat direction field starts a coherent oscillation. At this stage, if there exists baryon 
(or lepton) number violating operator in the scalar potential, field can have a nonzero 

phase rotational motion: — [arg(0)] 7^ 0, which means, a nonzero 0-number 

= i (0*0 - 0*0) ^0, (3.1) 

is produced. (The dot denotes a derivative with time t.) Since the field carries baryon 
(or lepton) number, this represents nothing but baryon (or lepton) number production. 

Among various flat directions, the LH U flat direction proposed by Murayama and 



Yanagida [25| is especially attractive and had been investigated in detail later as well [104 



105| , |101| , p5j |102| , |103|1 . The most interesting point is that it is closely related to the 



neutrino physics. Actually, as we will see, the most important parameter which determines 
the present baryon asymmetry is the mass of the lightest neutrino, m u i- Furthermore, 
the LH U flat direction is very special among various flat directions, since it is free from a 
serious problem of the Q-ball formation.^ 

Meantime, it has been pointed out [|104| , |97| , |98|| that the dynamics of a flat direction 
field is affected by thermal effects. Although production of baryon and/or lepton number 
usually takes place before the reheating process of the inflation completes, there is a dilute 
plasma even in this epoch. This dilute plasma can modify the scalar potential [07, 



affect the dynamics of 0, and change (actually, suppress) the resultant baryon asymmetry. 

In Sec. \5.'2[ therefore, we will perform a detailed analysis on the dynamics of the 
LH U flat direction field 0, including all the relevant thermal effects, in order to calculate 
the amount of produced lepton (and hence baryon) asymmetry. We see that evolution 
of the field is indeed affected by the thermal effects, and that the produced lepton 
asymmetry is suppressed by those thermal effects. On the other hand, we will find a 
very interesting result: the amount of produced baryon asymmetry is almost independent 
of the reheating temperature Tr, and is determined mainly by the mass of the lightest 



2 We do not discuss the Q-ball problem in this thesis. For details and references, see, for example, 



Refs. [106, 102, 107 1 
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neutrino vn v \. We show that the empirical baryon asymmetry hb/s ~ (0.4-1) x 10 -10 
predicts an ultralight neutrino m v \ ~ (0.3-1) x 10~ 9 eV for a wide range of reheating 
temperature Tr ^ 10 5 GeV. In Sec. [O] we show some explicit models to explain such a 
small mass of the lightest neutrino. 



Then we study in Sec. |3.4| what happens in this leptogenesis scenario via LH U flat 
direction when we introduce a gauged U(l) b-l symmetry. We will see that if the initial 
amplitude of the flat direction field exceeds the B — L breaking scale, we can avoid the 
suppression of the lepton-asymmetry production due to the thermal effects. 



Finally, in Sec. [3.5| , we show that if we assume a very small mass of the lightest neutrino 
as suggested by the present leptogenesis mechanism, the rate of the neutrino-less double- 
beta decay can be predicted in a high accuracy. 

A comment: we do not discuss other lepton-number carrying flat directions (see, e.g., 



early works in Refs. [ |108| , pTH , other general directions in Ref. ||104| , and a recent one in 
Ref. [|109|1 ) nor the case of multi-scalar manifolds |110fl . Notice that the LH U direction is 
special not only because it is a leptonic direction, but also because it is directly related 
to observable physics: the neutrino mass. 
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3.2 Lepton asymmetry from LH U flat direction 

In this section, we will discuss the dynamics of the LH U flat direction field in detail and 
calculate the amount of generated lepton asymmetry. 



3.2.1 The LH U flat direction 



Let us start by writing down the effective dimension-five operator in the superpotential 

777 . ) 

w = ^tfV = ou— ' ( 3 - 2 ) 

2 {H u ) 2M eS>i 
which induces neutrino masses m V i after the neutral component of the Higgs field H u 
obtains its vacuum expectation value (H u ) = 174 GeV x sin (3. Here, we have taken a 
basis where the neutrino mass matrix is diagonal. As we will see, the effective mass scales 
M eff j play crucial roles in the present scenario. Notice that they are directly related to 
the neutrino masses m V i by the following equation: 

Meff { = = 3 x 10 13 GeV ( —) . (3.3) 

m ui V m vi ) 

Here and hereafter, we take sin/3 ~ 1, for simplicity.^ 



We adopt the following scalar field fa as the Affleck-Dine flat direction 2j| for lepto- 



genesis, which was originally proposed in Ref. [25 



a -M)' i3A) 

where the factor y/2 is introduced to ensure the canonical kinetic term for fa. Here, we 
have chosen the neutral components of the Lj and H u for illustration, which is always 
possible by using the SU(2) L gauge transformation. We also used an U(l)y gauge rotation 
to make their phases to be the same. 

There are three flat directions (i = 1,2,3). (See Appendix |3.A] for more details.) 



Hereafter, we will suppress the family index i for simplicity. As we will see later, it turns 
out that the relevant flat direction fa for the most effective leptogenesis corresponds to 
the first family i — 1. Thus, our flat direction field is parameterized by one complex scalar 
field, ee fa. We will use the same symbol to denote the supermultiplet corresponding 
to the scalar field 0, hereafter. Its superpotential is then given by 

W = —^—<p\ (3.5) 



3 This is the case for tan/? > 1. Even for tan/3 ~ 1, the final lepton (and hence baryon) asymmetry 
changes (decreases) by a factor of 2. 



78 



while the scalar potential is given by 

Va = m > ? + mt„ M 4 + Hc -) + ik; w °- (3 - 6) 

Here, m$ and m 3 / 2 a m are SUSY-breaking mass parameters,^ and the last term directly 
comes from the superpotential Eq. (|3.5|). Here and hereafter, we assume that the SUSY 
breaking is mediated by the gravity. Therefore, it is expected that ~ m 3 / 2 ~ 1 TeV 
and \a m \ ~ 1. 

3.2.2 Initial amplitude 

The potential for is .D-flat, and it is also F-flat in the limit of m v — > 0. Thus, it is 
expected that the field can develop a large amplitude during the inflaton. It depends 
on the coupling between and the inflaton \ whether or not the field can have a large 
initial amplitude. Suppose that there are general nonminimal couplings of the field to 
the inflaton \ in the Kahler potential: 



K = 0V + X f X 



+ 



( c <t> 



+ H.c.) + ^X f X0V + " " " , (3.7) 



where and are complex and real coupling constants, respectively, and the ellipsis de- 
notes higher order terms which are irrelevant to the following discussion. During inflation, 
the energy density of the universe is dominated by the inflaton x- Thus, the total energy 
density of the universe is given by 3H^ u{ Mq ~ | | 2 , where H in f is the Hubble parameter 
during the inflation and F x is the F-component of the inflaton supermultiplet \- Because 
of this nonzero energy density of x, there appears additional SUSY-breaking effects during 
the inflation [IU1|. Using the scalar potential in the supergravity in Eq. ( |2.57| ), we obtain 



5V = 3 (1 - + M 2 ) f£f |0| 2 - (c,0 4 + H.c.) + ■ ■ • , (3.8) 

where we redefined the phase of to include the phase of F x . A crucial observation here 
is that the field gets a large mass term of the order of the Hubble parameter. (We 
will call it a Hubble mass term.) In general, the coefficient of this Hubble mass term, 
3(1 — bcf, + |c^| 2 ), is expected to be order unity. If it is positive, therefore, the field is 
driven exponentially towards the origin during the inflation, and the present mechanism 



4 Here, we have included the possible contribution from the /x-term, W = fiHuHd, which gives a mass 

2 

<t>- 



term (l/2)/i 2 |</>| 2 , in the mass m 2 
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cannot work. (Notice that this would happen even if the Kahler potential is minimal, i.e., 
Ccf, = = 0.) Hence, we will assume that it is negative, and hereafter take 



5V = -cuHl^ + ^g- (a H 4 + H.c. j , (3.9) 
with ch — | a# | — 1. 

The negative Hubble mass term causes an instability of at the origin. From Eqs. (|3.6|) 
and (|3.9|), the minima of the potential are given by 



|0| ~ ^JM eS H ini , (3.10) 

ia\ - arg(a H ) + (2n + l)vr 
arg(0) ~ , n = 0---3, (3-11) 



where we have used m^, m 3/ / 2 |a m | <C H in {. Since curvatures of the potential around the 
minimum along both the radius and phase directions are positive and of the order of the 
Hubble parameter H in f, the flat direction runs to one of the four minima from any given 
initial value and is settled down it. 

Here, we comment on another possibility that the flat direction field does not receive a 
mass term of the order of the Hubble parameter, which occurs with a Kahler potential of 



the no-scale type supergravity ||111|| . The evolution of the field in this case was discussed 
in Ref . [|105|| , and it turns out that the initial amplitude is likely to be given by Eq. (|3.10|) 
in this case as well. 



3.2.3 Evolution of the flat direction 

Before discussing the detailed dynamics of the field, we first roughly describe the evo- 
lution of and note the relevant epoch for the leptogenesis. During the inflation, the 



field takes a large value ~ ^M eS H in f as discussed in the previous subsection. After the 
end of inflation, the value of gradually decreases as the Hubble parameter H decreases 
and then, at some time, the starts its coherent oscillation. As we will see, the net lepton 
number is fixed when the flat direction field starts the coherent oscillation. 

Hereafter, we will assume that the net lepton number is fixed before the reheating 
process of the inflation completes, namely, during the epoch when the energy density of 
the universe is dominated by the oscillating inflaton \. (This assumption will be justified 
after we calculate the amount of the generated lepton asymmetry in Sec. |3.2.4j .) 
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3. 2. 3. a potential for 



Let us first show the total effective potential for the flat direction field relevant to the 
leptogenesis: 

Vtoui = \ml-H 2 + £ c k f 2 k T 2 

V fk\<t>\<T 



+a g a s (T) 2 T 4 In 




J>2 



(3.12) 

We explain each term in turn. 

First of all, the F-term potential and the SUSY-breaking terms give rise to the scalar 
potential at vacuum Vo, which is given by Eq. (|3.6| ). 

Next, because the energy density of the universe is still dominated by the inflaton x 
during the inflaton-oscillation epoch, the field receives additional SUSY-breaking terms 
caused by the finite energy density of x- 

5V H = -c H H 2 \<f)\ 2 + -§r- (a H 4 + H.c.) . (3.13) 



This can be derived in the same way as Eq. (|3.9|) . Recall that we have taken ch — \o>h\ — 1- 
We will call the second term Hubble A-term. 

Finally, the rests in Eq. (|3.12| ) correspond to the thermal effects which we discuss 
now. The crucial point here is that the decay of the inflaton x occurs during the coherent 
oscillation of x, while it completes much later than the beginning of the oscillation. Thus, 
even in the x oscillation period, there is a dilute plasma consisting of the decay products 
of the inflaton x, although most of the energy density of the universe is carried by the 
coherent oscillation of x- The temperature of this dilute plasma is given by p8| 

T ~ {TlM G H) llA . (3.14) 

Notice that in this dilute plasma the temperature decreases as T oc if 1 / 4 , rather than as 
T oc H 1 / 2 like in the usual radiation dominated universe. 



This dilute plasma has crucial effects on the dynamics of the field ||104| , P8 
There are basically two possible thermal effects. First, the fields ip^ which couple to (f) are 
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produced by the inflaton decay and/or by thermal scatterings if their effective masses are 
less than the temperature, 

m eStk = f k \4>\ < T , (3.15) 

and hence they must be included in the plasma. Here, fk denote the Yukawa or gauge 
coupling constants between the flat direction and ipk (we take fk real and positive). 
Therefore, the flat direction receives a thermal mass of order ~ /feT, if the condition 
Eq. ( |3.15|) is satisfied. The induced mass term is given by [0 

< h = E c fc / fc 2 T 2 |0| 2 , (3.16) 

fkW<T 

where Ck are real positive constants of order unity (see below). The summation means 
that only the fields in thermal plasma [i.e., ipk whose couplings fk satisfy the condition 
Eq. ( |3.15|) 1 can induce the thermal mass for 0. 

We include all the couplings relevant for the flat direction H u = L\ = 0/a/2 in the 
SUSY standard model, i.e., the gauge couplings for the SU(2)l and U(l)y gauge groups, 
and Yukawa couplings for up- type quarks and charged leptons. The couplings fk, which 
are redefined so that the effective masses for the fields -0fc become = fk\4>\, are 

listed in Table. |3.1| with the coefficients Ck of the thermal mass for 0. Among the gauge 
supermultiplets, the flat direction does not couple to one U(l) group which remains 
unbroken by the condensation of 0. Thus, there are one Z-like and two PU-like massive 
gauge multiplets as in the SUSY standard model. For the couplings of the charged leptons, 
we should be careful. The flat direction receives its thermal mass from only one linear 
combination of charged leptons. The effective Yukawa coupling yn is determined by the 
Yukawa couplings for the charged leptons and the mixing matrix among gauge eigenstates 
of the neutrinos. (Note that the relevant is the flattest direction which corresponds to the 
first family v\ of mass eigenstates.) Taking large mixing angles both for atmospheric and 
solar neutrino oscillations suggested by the recent data [[TO], [□], |T2"|| , we find yn — 0(0.1- 



1) x y T . We will see in Sec. |3.2.4| that the Yukawa coupling for the up quark induces 



a relevant thermal effect, while other couplings turn out to be too large to satisfy the 
condition Eq. (|3.15|) in most of the parameter space. 



Next, there is another thermal effect which was pointed out in Ref. [98]. Along the 
LH U flat direction, the SU(3)c gauge symmetry remains unbroken and hence gluons and 
gluinos are massless. Furthermore, the down type (s) quarks also remain massless since 
they have no direct coupling to the field. These massless fields produce the free energy 
which depends on the SU(3)c gauge coupling constant gs, and we obtain an effective 
potential 5V oc gs(T) 2 T A at two loop level.[] At first sight, there seems no dependence on 



J For the coefficient of this free energy, see, for example, Ref. [112 



82 













fk (coupling for 0) 


^91 + 91 


92 

V2 


^= ( a = up, charm, top ) 


Vli 
V2 


Ck 


1 


1 


3 


1 


4 


2 


4 


4 



Table 3.1: The couplings of the flat direction field to other fields and the coefficients Ck 
of the thermal mass of induced by these fields. g% and #2 denote the gauge couplings 
for the U(l)y and SU(2)i gauge groups, respectively, and y Uia (a = up, charm, top) and 
Uli are the Yukawa couplings for up-type quarks and charged lepton, respectively. 



the field in this free energy. However, we can see this is not the case by means of the 
following arguments. 

The evolution of the running coupling g$ is given by 

Ar^ = 4- 2 (- 3 C 2+ j:c m y (3,7) 

where C2 = 3 and C{Rj) = 1/2 for fundamental representations Ri. Notice that the 
evolution changes due to the decoupling when the scale \i passes through a mass of a 
field rrii. Since the up type (s) quarks get large masses from the couplings to the field 
fu,a — Uu,a/V2 (a = up, charm, top), they change the trajectory of the running at 
H = = / M) o|0|j as shown in Fig. pTT]. Hence, if f u ,a\4>\ > T, the running coupling 
g s(T) depends on the value of the field: 

9S (T) = 9s (r) U + ^£ Y. C(fl„)ln(%^), (3.18) 

6A7T fu,a\4>\>T V 1 ) 

where Mjj is the ultraviolet scale where gs is fixed. Then, there appears an additional 
potential through the modification of the gauge coupling constant in Eq. (|3.18|) 



|2 N 



5V^ = a g a s (TfT i r- £ C(R u , a ) In ( J^L j , ( 3 .19) 

where a g is a constant which is a bit larger than unityQ and cts(T) = 5 f s(T)^ =0 / (4vt). (We 
will omit the subscript = for simplicity.) Hereafter, we take the factor (2/3) J2C(R) 
to be unity since it does not change the result much.[] 



By using the free energy at two- loop order given in Ref. |112| , a g is given by a g = 1.125 in the case 
of the LH U fiat direction. The author thanks Masaaki Fujii for helpful discussion. 

7 At least the top Yukawa coupling / U) top always satisfies / Uj top|0| > T before the oscillation of <p. 
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Figure 3.1: A schematic behavior of the SU(3)c gauge coupling g s . The dashed line 
represents the running coupling when the (p held does not have a vacuum-expectation 
value. 

3.2.3.b evolution of 6 



Now we eventually obtain the total effective potential Eq. ( |3.12| ): 

Kotai =V Q + 6V h + 8Vl h + 5Vt . 
The evolution of <ft is described by the equation of motion with this V tota \ as 

distal 







(3.20) 



(3.21) 



where the dot denotes a derivative with time. During the inflation, there is no thermal 
plasma and hence no thermal potentials. The vacuum expectation value of the field at 
this stage is given by Eqs. (|3.10|) and ( |3.11| ). Let us call the valley to which the <fi field 
rolls down n = n Q . 

After the inflation ends, the inflaton x starts to oscillate and its decay produces a 
dilute plasma. However, the potential is still dominated by Hubble-induced terms [in 
Eq. (|3.13|) 1 and |0| 6 term [in Eq. (|3.6| )1 at the first stage of the oscillation. Thus, the flat 
direction field <fi is trapped in the following minimum: 

(3.22) 



arg(0) 



M eS H, 
- arg(q g ) + (2n + l)n 
4 



(3.23) 



Thus, the resultant baryon asymmetry changes by only a factor of V3 at most. [See Eqs. ( p. 34 ) and 
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This is because the curvatures around the minimum along both radius and phase directions 
are of the order of H also in this period, and hence the field always catches up the 



instantaneous minimum 104 



Then, as the Hubble parameter decreases, the negative Hubble-induced mass term is 
eventually exceeded by one of the following three terms in the potential:^ 

+ E ^fl T2 +Ws(T)^. (3.24) 
h\4>\<T l<?l 



H 



2 < 



As we show now, it is this time when the oscillation of starts. Let us denote the Hubble 
parameter at this time by H osc . The evolution of the after H ~ H osc depends on which 
term in Eq. ( ^.12[ ) dominates the effective potential. There are basically three cases; the 
potential is dominated by (i) m 2 |0| 2 term, (ii) T 2 |0| 2 term, or (iii) T 4 ln(|0| 2 ) term. 

(i) First, if the potential is dominated by the m 2 |0| 2 term, the equation of motion 
Eq. fl3.21| ) is given by 



+ 3#0 + m 2 d 



0. 



(3.25) 

0) and the amplitude 



It is clear in this case that the field oscillates around the origin (</ 
of the oscillation damps as \(f)\ oc H oc t . 

(ii) Second, when the potential is dominated by the thermal mass term c^flT 2 |0| 2 
the equation of motion becomes 



+ 3#0 + c k f£T z <f> ~ . 
This equation can be solved analytically and the solution is given by 

m = 



(3.26) 









[ Z 2/3 J 2/3W 


+ 02 


[ 2 2/3 J ^/3(^)] 



4 /2 \ 1 / A 

I (\clftM G Tl t 3 



(3.27) 



where (f>± and 02 are constants and J v is the Bessel function. We see that the oscillates 
around the origin in this case as well. The time scale of the oscillation is ~ (f^MaT^) 
and the amplitude of the oscillation is damped as |0| ~ t~^'^ ~ 

(iii) The third case is given when the a g a 2 s T A ln(|0| 2 /T 2 ) term dominates the potential. 
If we neglect the time dependence of T 4 , the damping rate of the oscillation amplitude due 
to such a flat potential, V ~ ln(|0| 2 ), can be estimated by using the virial theorem, and it 



8 Here, we have compared (\/4>)(dV/d4)* ) instead of the potential terms themselves, since the evolution 
of 4> is governed by the equation of motion Eq. (3.21). 
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is given by |0| cx H 2 oc t 2 |TTH]. In the actual case, however, the potential itself gradually 



decreases with time as T 4 oc t 1 . We have numerically checked that the amplitude damps 
as |0| oc H a oc t~ a with a ~ 1.5. 

Thus, in any case, the field starts a coherent oscillation. Notice that, in all the above 
cases, the damping rate is faster than the rate before the beginning of the 0's oscillation 

(|0| oc tf 1 / 2 oc r 1 / 2 ). 

Now let us estimate the cosmic time t osc when the field starts its oscillation, or 
equivalently the Hubble parameter at that time, H osc = (2/3)t~^,. As we will see in 
Sec. |3.2.4| , this oscillation time plays an important role to determine the final lepton 



asymmetry. H osc is determined by the term which exceeds the Hubble mass term at first. 
(Therefore, again, there are three cases.) 

(i) If the soft mass term dominates the potential at first, H osc is just given by H osc = 

(ii) The case when the thermal mass term dominates at first is a bit complicated. Let 
us first consider a simple case in which there is only one field ipi with a coupling f\ . From 
Eqs. ( |3.15|) and ( |3.24|) , it is found that the following two conditions must be satisfied in 



order to cause the 0's oscillation by the thermal mass term 

1 M G T 2 
ft M 2 S 

H < (clffMeT 2 )"- , (3.29) 



H < 71^71^ ( 3 - 28 ) 

1/3 



where we have used the relations \<f>\ ~ VM eS H and T = (T|M G if ^Z 4 . The meanings 
of these conditions are as follows. If the coupling f\ is small enough, the ?p\ can easily 
enter in the thermal plasma {fi\<f)\ < T), but the induced thermal mass for is smaller 
than the Hubble-induced mass H at the beginning (cif 2 T 2 < H 2 ). However, since the 
temperature decreases more slowly than the Hubble parameter H (T oc if 1 / 4 ) as the 
universe expands, the thermal mass eventually exceeds H (c\f 2 T 2 > H 2 ). Then the 0's 
oscillation starts at the time when the condition Eq. ( |3.29p is satisfied. On the other 
hand, if f\ is very large, the would-be thermal mass can be large enough to exceed the 
Hubble mass [c\]\T 2 > H 2 ). However, the thermal mass does not appear until the \<f>\ 
becomes small enough to satisfy fi\<f)\ < T. In this case the starts to oscillate at the 
time when the condition Eq. ( |3.28| ) is satisfied. Therefore, we find the Hubble parameter 
H nsr to be 



H osc = min 



lM 0T 2 R ( r 2 f A MT 2\^ 



(3.30) 



It is easy to apply the above discussion to the case of more than one couplings. If there 
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is another coupling fa which can satisfy the both conditions Eqs. ( |3.28| ) and ( (3.29|) earlier 
than fx, the flat direction starts its oscillation earlier. Therefore, the Hubble parameter 
is given byf] 



H osc = max (Hi) 



where 



Hi 



mm 



1 M G T\ 
ft M ( 



cff 



(clftMaT, 



1/3 



(3.31) 



(3.32) 



(iii) The third case is again given when the a g a 2 T A ln(\<f)\ 2 /T 2 ) term dominates the 
potential at first. In this case, it is found from Eq. ( |3.24| ) that the oscillation time of the 
field is given byQ 



H n 



a s T, 



R 



a g M G 



M, 



cff 



1/2 



After all, the Hubble parameter at the time when the 
lation is given by 



(3.33) 

field starts its coherent oscil- 



H osc ~ max 



m^, Hi, a s T R 



OgMc 



1/2 



(3.34) 



where Hi is given by Eq. ( |3.32j) . We stress here that the thermal effects always make the 
oscillation time earlier. (Notice that if there is no thermal effect, the oscillation time is 
just given by H osc = m^.) As we shall show soon, this "early oscillation" drastically affect 
the amount of the produced lepton asymmetry. 



3.2.4 Lepton asymmetry 

Now we are at the point to calculate the lepton asymmetry produced by the field. 
Since the field carries lepton charge, its number density is related to the lepton number 
density n L as 



TIt, 

2 



^(0*0-0*0) . (3.35) 



From Eqs. ( |3.12| ), ( |3.21| ) and ( |3.35| ), the evolution of til is described by the following 
equation: 



h L + 3Hn L = ra^/ 2 Im ( am( fA + -JL—\ m (a H 4 ) . (3.36) 



9 We assume here hierarchical couplings and neglect effects of the summation. 

10 The running coupling as(T) here is obtained by solving the equation H 2 — a g ots(T)T 4 /|0| 2 itera- 
tively. 
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The production of the lepton asymmetry occurs as follows. Suppose that the original 
A-term would vanish, i.e., 7713/2 a m — 0. If this is the case, the flat direction is always 
trapped in one of the valleys induced by the Hubble A-term [see Eq. ( [3.23j )1 during the 
both periods t < t osc and t > t osc , and the direction of the valley does not change with 
time. Therefore, there is no force which causes the motion of along the phase direction 
and no lepton-number asymmetry is produced (i.e., the right-hand side of Eq. ( |3.36| ) 
vanishes.) However, we have the original A-term, and the phase of is kicked by the 
relative phase difference between a m and a# in Eq. ( |3.12|) . The phase of changes during 
its rolling towards the origin, because the Hubble parameter H decreases and the direction 
of the true valleys changes with time. Therefore, the original A-term, which corresponds 
to the first term in Eq. ( |3.36| ), plays a role of the source of the lepton asymmetry. 

One might wonder if the Hubble A-term, which corresponds to the second term in 
Eq. ( 3.36Q , gives a larger contribution to the lepton asymmetry since H 3> 7723/2. How- 
ever, since the almost traces one of the valleys determined mainly by the Hubble A-term, 
lm(a#0 4 ) in Eq. ( |3.36| ) is highly suppressed compared with lm(a m 4 ). In fact, we have 
found numerically that the contribution of the second term in Eq. ( |3.36| ) is always com- 
parable or less than that from the first term. Thus, we neglect the lepton asymmetry 
produced from the Hubble A-term in our analytic calculation, for simplicity. By integrat- 
ing Eq. (|3.36| ), we obtain the resultant lepton number at the time t, 

'R 3 n L ](t) ~ ^-j'dt'RH-m^m^) , (3.37) 

where R denotes the scale factor of the expanding universe, which scales as R 3 oc H~ 2 oc t 2 
in the universe dominated by the oscillation energy of the inflaton x- We can see that 
the total lepton number increases with time as R^n^ oc t until the oscillation of starts 
(t < t osc ), since 4 oc H 2 and hence i? 3 4 ~ const in this stage. 

On the other hand, after the starts its oscillation, the production of lepton number 
is strongly suppressed. This is because Im (a m 4 ) changes its sign rapidly due to the 
oscillation of 0, and also because the amplitude of 0's oscillation is damped as i? 3 4 ~ 
t~ n with n > 1. [See discussion around Eqs. ( ^.25[ )-( p.27| ).] Therefore, the net lepton 
asymmetry is fixed when the oscillation of starts. The generated lepton number at this 
epoch (t = t osc ) is given approximately by 



^l(^osc) = ^7^-Im [a m 4 (t 



OSC ; 



tc 



7: (rn 3/2 \a m \J M cS H osc 5 ph , (3.38) 



3 

where £ p h ~ sin(4arg0 + arga m ) represents an effective CP violating phase. Here, we 
have used 4>(t osc ) ~ \J M eS H osc and t osc = {2/?>)H^ c . Recall that we have assumed that the 
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production of net lepton asymmetry occurs before the reheating process of the inflation 
completes, i.e., H osc > T x where T x is the decay rate of the inflaton \- Thus, the lepton 
number when the reheating process completes (t — tu,H~ T x ) is given by 

n L (tn) = 



H osc 



2 



n L (t osc )(ir^) ■ (3.39) 

Then, the lepton-to-entropy ratio is estimated as 

n L M eS T R f m 3/2 \a r . 



'in 



s 12M% \ H Q 



V , (3.40) 



when the reheating process of inflation completes. Here, we have used 3MqT x = p ra d(tR) = 
(3/4)T R x s(t R ). 

As explained in Sec. pL .44 this lepton asymmetry is partially converted into the baryon 
asymmetry through the sphaleron effect. Thus, after all, the present baryon asymmetry 
is given by 

— ~ 0.35 x — 

s s 

n nor. M e gT R / m 3/2Wm\ \ x ro A1\ 

= a029 x ^r j 6 *- (3 - 41) 



From Eq. (^4ll) , we see that the baryon asymmetry is a monotonically increasing 



function of the scale M eS . [Notice that the H osc depends on M eS as shown in Eq. ( |3.34j ).l 
Since M e fj is directly related to the neutrino mass m u as M c g = (H u ) 2 /m„, this means 
that the baryon asymmetry becomes larger as the neutrino mass m u decreases. Therefore, 
as mentioned at the beginning of this chapter, the most effective flat direction corresponds 
to the lightest neutrino u 1 , i.e., the first family field 4>/V2 = L\ = H u . 

The factor H~l in Eq. ( p.41|) represents the suppression of the lepton asymmetry 



due to the thermal effects. When the 0's oscillation is caused by the thermal effects, 
i.e., if osc > m^, the production of lepton asymmetry stops earlier than the case without 
thermal effects, and hence the resultant baryon asymmetry becomes suppressed. 

We should also note that the amount of the produced baryon asymmetry is propor- 
tional to the gravitino mass, m 3 / 2 . This is because the force which rotates the field along 
the phase direction comes from the SUSY-breaking A-term (m3/2/8M e ff)(a m 4 + H.c). 
Thus, the resultant baryon asymmetry would be very much suppressed if the gravitino 
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Figure 3.2: The contour plot of the baryon asymmetries n^/s in the m v \-T R plane. The 
lines represent the contour plots for %/s = 10~ 9 , 10~ 10 , 0.4 x 10~ 10 , 10~ n , and 10 -12 
from the left to the right. 

is light, like in gauge-mediated SUSY breaking models lj45H.PI This is the reason why we 
have assumed gravity- mediated SUSY-breaking and 1123/2 — 1 TeV. 

Fig. [3.2| shows the contour plot of the produced baryon asymmetry in the m v \-Tn 
plane, given by the analytic formula Eq. ( |3.41|) . Here, we have taken = 777.3/2 = 1 TeV, 
\a m \ = 1 and 5 p h = 1, and we have used the relation m v i = (H u ) 2 /M e s. (Although the 
big-bang nucleosynthesis constraint on the gravitino abundance requires Tr ^ 10 9 GeV 
for m 3 /2 = 1 TeV, we have plotted the contour in Fig. |3]^ up to T R < 10 11 GeV, keeping 
in mind that even Tr ~ 10 11 GeV can be allowed for a slightly heavier gravitino, say, 



m 



3/2 



3 TeV. See Sec. 1.5 



A remarkable observation here is that the present baryon asymmetry %/s is deter- 
mined almost independently of the reheating temperature for a wide range ofTR> 10 5 GeV. 
In particular, for a relatively high reheating temperature T R ^ 10 8 GeV, the baryon asym- 



n If we introduce a gauged U(l)s-i, however, the LH U flat direction can produce enough lcpton 
asymmetry even with a light gravitino. See Sec. 3.4. 
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metry derived from the Eqs. ( |3.34| ) and ( |3.41 ) is given by the following simple from 



- - 10- U ^ x f-^U ( ^4¥) . (3.42) 



10- 8 eVJ V 1 TeV 
The reason why it is independent of the reheating temperature Tr is that the oscillation 
time H osc is determined by the thermal potential T 4 ln(|0| 2 ) in the higher temperature 
regime and becomes proportional to the reheating temperature T R . [See Eq. fl3.34 ).] Thus, 



the Tr dependence is canceled out in Eq. (|3.41| ) Even in the lower reheating temperature 
region 10 5 ^ Tr ^ 10 8 GeV, where H osc is determined by the thermal-mass term poten- 
tial T 2 \(p\ 2 , Tr dependence is still mild, i.e., ur/s oc 1# . The "reheating-temperature 
independence" discussed here is a very attractive feature of the present mechanism since 
the produced baryon asymmetry crucially depends on Tr in many other baryogenesis 
scenarios. 



For the analytic calculation in Fig. |3.2| , we have taken 5 p h = 1. It is expected that 



<5 P h ~ ^(1); sa Y) <5ph — 0.1-1, unless there is an unnatural cancellation between arg(a r , 



and arg(a#). Thus, it is found from Fig. |3.2| that the present baryon asymmetry in our 



universe hr/s ~ (0.4-1) x 10 10 suggests an ultralight neutrino of a mass, 

m vX ~ (0.1 -3) x 10~ 9 eV, (3.43) 
for T R > 10 5 GeV and 5 ph ~ 0.1-1. 

• remarks 

Before closing this section, we give several comments. First, let us justify our assumption 
that the production of the lepton asymmetry stops before the reheating process of the 
inflation completes, namely, H osc > T x . We can see this condition is satisfied in all 
the parameter region in Fig. |3]2] by using the fact that H osc is bounded from below, as 
shown in Eq. ( |3.34j ). For example, it is enough to show asTR(a g Mc / M e sY^ 2 > T x . By 
using T x = (7r 2 g*/90) 1//2 T!/MG, it is found that this condition is indeed satisfied for 
T R < 10 14 GeV x (m^/10- 9 eV) 1 / 2 . 

Next, let us investigate the behavior of Fig. [O] in more detail. Fig. |3.3| shows the 
dependence of the baryon asymmetry on the Yukawa coupling of the up quark. Notice that 
the up-quark Yukawa coupling is given by y UtUp = ^ up / (H u ) = m up /(174 GeV x sin/3). 
In the previous figure |3]^, we have taken m up = 3 MeV (and tan/? = 10). In Fig. |3~3| , we 
have taken m up = 1 MeV and m up = 5 MeV for comparison. It is found that the final 



12 



Precisely speaking, H osc depends on Tr through the running coupling as{T) even in this case. 



However, as can be seen in Fig. 3.2, this dependence is negligibly mild. 
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Figure 3.3: The contour plot of the baryon asymmetries ris/s in the m v -\—Tn plane. The 
lines represent the contour plots for %/s = 10~ 9 , 10~ 10 , 10 -11 , and 10~ 12 from the left 
to the right. The solid lines are the same as those in Fig. p.2| (m up = 3 MeV). The left 
(right) dashed lines correspond to m up = 5(1) MeV. 



baryon asymmetry changes at most one order of magnitude. However, the dependence on 
the reheating temperature is still mild, and hence our conclusion does not change much. 



Meantime, we can also see from Fig. |3]3| that the up-quark Yukawa coupling is the only 
coupling which can give a dominant contribution to the early oscillation of </>, since below 
this "up-quark Yukawa" region the oscillation is caused by the soft mass term (where 
ub/s oc Tnm~i), while in the higher Tr region the oscillation is caused by the thermal 
logarithmic term T 4 ln(|^| 2 ) (where Ub/s oc m v \ ). This is because other couplings fk 
are too large and the fields with those couplings cannot satisfy the condition fk\<f>\ < T. 

We have also solved the equation of motion ( |3.21| ) numerically with the full scalar 
potential Eq. fl3.12| ). The result is shown in Fig. |3.4| , where the final baryon asymmetry 
n B /s ~ 0.35 x n L /s is plotted. (n L is given by Eq. ( |3.35| ), n L = (z/2)(</>*</> — 4>*4>).) Here, 
we have taken = m 3 / 2 |a m | = 1 TeV, arg(a m ) = 7r/3 and arg(a^) = 0. As for the 
initial condition, we have taken H in f = 10 12 GeV, and put the <fi field at the minimum of 
the potential during the inflation. From Fig. ET1 we confirm that the analytic estimation 
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discussed above reproduces very well the result obtained by the numerical calculation. 

Finally, we comment on the possibility to avoid the early oscillation caused by the 
thermal effects discussed in this section, which has been recently pointed out in Ref. ||114 |. 



The point is that, if the scale of the inflation is very low, the Hubble parameter during 
the inflation if in f can be below the H osc which is given by Eq. ( |3.34j ). If H m t < H osc , our 
(implicit) assumption that the Hubble parameter of the universe gradually decreases from 
Hi n f to H osc breaks down. What happens in this case is that the flat direction field <fi starts 
its oscillation just after the inflation ends, i.e., at H ~ H in f. Thus lepton asymmetry is 
also fixed at this time. The final baryon asymmetry is then given by replacing H osc in 
Eq. ( ^4TD with H m{ : 



s 



, , 29x « (3.44) 

H in{ <Hosc M G \ H mi ) 



It was shown [|1 14f that H m f ^ H osc is indeed satisfied for very low ifi n f and high enough 



reheating temperature T R , e.g, H mi ~ 10 5 GeV and T R ^ 10 7 -10 8 GeV. For example, in 
an extreme case where H m { ~ and Tr ~ (^g^/QO)" 1 ^^ 'M G H mi (sudden reheating), 
baryon asymmetry is given by [|1 14 



2- K 0.3 x 10- ( " 2 W V . (3.45) 



m v \ J V i?i n f / VI TeV 

Thus the empirical baryon asymmetry can be explained even with m v \ ~ 10~ 4 eV in this 
extreme case. 
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Figure 3.4: The plots of the baryon asymmetries %/s in the m v \-Tji plane, which are 
obtained by numerically solving the equation of motion Eq. (|3.21|) . The regions with 
points show the result of the numerical calculation. They represent 10 -9 > rig/s > 10~ 10 , 
1(T 10 > n B /s > 10 -11 and 10 -11 > n B /s > 10" 12 from the left to the right. The solid 
lines are the same as those in Fig. |3.2| , which represent the contour plots for %/s = 10 -9 , 
1(T 10 , 10 -11 , and 10" 12 from the left to the right, given by the analytic formula Eq. ( |3.41|) . 
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3.3 Some models 



In Sec. p.2.4| we have shown that the baryon asymmetry in the present universe predicts 



the mass for the lightest neutrino in a narrow region, m v \ ~ (0.1-3) x 10 9 eV. Together 



with the neutrino masses required to solve the solar |TT[ and atmospheric |10| neutrino 
anomalies, this suggests a very large mass hierarchy between the lightest and the heavier 
two neutrinos. 

Although the present leptogenesis mechanism via LH U direction does not rely on the 



origin of the effective operator Eq. ( |3.2|) , it is natural to consider the seesaw mechanism [13 
as the origin of that operator. In this case, the effective mass scale M e g is related to the 
masses of the right-handed neutrinos M;. Then, one might wonder if the discussion so 
far is reliable, since successful leptogenesis requires the effective mass scale M e g to be far 
beyond the Planck scale M G , i.e., M eff = (H u ) 2 /m ul ~ 3 x 10 22 GeV x (m^/10" 9 eV) -1 . 
Here, we first comment on this point. In the framework of the seesaw mechanism, the 
neutrino masses are given by Eq. ((2^ 



m ui = - V hkih ki ^ Hu ) , (3.46) 

k M k 

where we have rotated the Yukawa coupling so that the neutrino mass matrix becomes 
diagonal [see Sec. ^J] and Eq. ( |2.12j )1. Then, from Eqs. (|3.3j ) and (|3.46|) we obtain 

-i 



M eff ee M effil = - I £ • (3.47) 

Therefore, the large mass scale M e g (or equivalently small neutrino mass m v \) can be 
obtained by taking the Yukawa couplings hu to be small, keeping the masses of the 
right-handed neutrinos Mi to be below the Planck scale. 

Nonetheless, the reader might still consider that the mass scale M e g is too large. In 
other words, the mass of the lightest neutrino m v \ ~ (0.1-3) x 10~ 9 eV might seem to 
be too small, compared with the masses of the two heavier neutrinos, m v ~ (9(10~ 3 )- 
(9(10 _1 ) eV suggested from the recent neutrino oscillation experiments |K| [11], [L^]. In 
this section, we propose some models in which such an ultralight neutrino (or large mass 
scale M e ff) can be naturally implemented. 

3.3.1 Model for the ultralight neutrino with a broken discrete 
symmetry 



Here, we show an explicit model based on a Froggatt-Nielsen (FN) mechanism pS 



in 



which the required large mass hierarchy is naturally obtained. We adopt a discrete Z 6 as 
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the FN symmetry instead of a continuous U(1)fn- We see that the discrete symmetry is 
crucial to produce the required large mass hierarchy in the neutrino sector. To see this, 
let us first recall the case of the continuous U(1)fn, which is described in Sec. [2.3.2| . In 
that case, the mass matrix of the neutrino in Eq. ( |2.48| ) gives rise to a mild hierarchy 



m vl : m u2 : m v3 ~ e 2 : 1 : 1 = 0(1O~ 2 ) : 0(1) : 0(1) 



(3.48) 



and hence it cannot explain the large hierarchy required in the present scenario. To 
change the above point, we suppose that the broken FN symmetry is not a U(1)fn but a 
discrete symmetry Z n with n = 2d. Then, the mass matrix for the right-handed neutrino 
Mr changes into the following form; 



£12 


e c+d 


63 


e b+d 


\ 


£22 


e 2c 


£23 


£ b+c 




£23 


£ b+c 


£33 


e 2b 


J 



(3.49) 



M R = M ^e c+d 

\^l3B b+d 

Here, we have assumed d(= n/2) > c > b > 0. Notice that the Majorana mass for N\ 
is no longer suppressed by the power of e, which is a basic point to yield an extremely 
small neutrino mass m v \. Although the mass matrix Mr is modified, the structure of the 
neutrino mass matrix looks similar to the original one: 

m 1 



m. 



m 



°w E m ° 



-2a 



(HuY 



Mn 



,2a 



Mn 



( e 

1 

V 

/ 0{e 2 
0(e] 

V 0(e] 




r 1 











(i 


1 




) 










0(e) 0(e 
0(1) 0(1) 
0(1) 0(1) 



(3.50) 



However, in spite of the large components of the mass matrix, we see that one of the 
mass eigenvalue of this mass matrix strongly suppressed as ~ £ 2 ( 1 + d ) . (This suppression 
is also understood directly by taking the determinant of the above mass matrix.) Thus, 
the required large mass hierarchy between m v \ and 77^2,3 can be naturally obtained by 
taking d = 3 (Z%). 

Notice that the neutrino mass matrix Eq. ( [3.50] ) can still lead to a large v jX -v T mixing 
angle. The mass hierarchy of the charged leptons also remains unchanged.0 Hence, this 
model preserves all the desirable features in the continuous U(1)fn model, explaining the 
mass of the ultralight neutrino m v \ as well. 



13 This is because the charges of the Yukawa couplings of the charged leptons are smaller than n = 6, 
Q(^ i )+Q(i i )<(a + 3)<6. 
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Figure 3.5: The plot for r = 8m 2 ol / <5m 2 tm and m v \ in a Froggatt-Nielsen model with a 
discrete Z 6 symmetry. 



To demonstrate our point, we randomly generate 0(1) couplings ijij and £y. Namely 
we calculate the mass matrix for neutrinos, taking the magnitudes of the couplings ylj 
and £y to be in a range 0.5 - 1.5 and their phases to be - 2n. We also take e = 0.05 - 0.1 
randomly.0 From the obtained neutrino mass matrix, we calculated the mass eigenstates 



m h 
8ml 



and the mixing matrix £7^,0 and then required the parameters 8m: 



2 

sol 



in 



i/3 



m 



v2i 



"atm ~~ m v2 " l vll 0111 il/ atm 

consistent with the data of neutrino-oscillation experiments |T 
large-angle MSW solution to the solar neutrino deficits, since it is the only solution allowed 



m. 



sin 2 26 a 



4|^ 3 | 2 (1 



Ufa] 2 ), and tan 2 # 



sol 



11 



U e2 /U el \ 2 to be 
121. Here, we take 



by the recently announced KamLAND results |fL2| , and adopt the following constraints 



0.9 < sin 2 2# atm < 1.0 



r = 8m 2 ol /8m 



2 

atm 



0.008 -0.15 



We also required 



tan 2 9 sol = 0.2 - 0.9 



\U e3 \ <0.15, 



(3.51) 



(3.52) 



to satisfy CHOOZ limit [|116|| . Fig. |3.5| shows the obtained mass of the lightest neutrino, 
m v \. We can see that an ultralight neutrino of mass m v \ ~ (0.1-3) x 10~ 9 eV is naturally 
obtained. 



14 A similar calculation was done in Ref. 
15 The matrix U a i is defined in Sec. 2.1 



where they adopted the continuous U(1)fn model. 
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3.3.2 Large effective masses for right-handed neutrinos 



In Sec. |3.3.1| , we have shown a model in which an ultralight neutrino with mass m v \ ~ (0.1- 
3) x 10~ 9 can be naturally implemented. Here, we discuss another possibility. The point 
is that the large effective mass scale M e g- is required only during the leptogenesis, not 
necessarily at the true vacuum. Therefore, if the effective masses of the right-handed 
neutrinos in the early universe are dynamical values and they are different from the 
values in the true vacuum, we might naturally obtain a large effective mass scale M e g 
during the leptogenesis. 

Suppose that the right-handed neutrino masses Mj are given by a vacuum expectation 
value of some field X with a superpotential 

W = ^Z i XN i N i , (3.53) 

which gives Mi = £j (X). If (X) takes a large value during the leptogenesis, therefore, Mj 
can be much larger than the values in the true vacuum. To demonstrate this point, let 
us consider that the field X is responsible to the Peccei-Quinn symmetry breaking ||117 . 



Then, as we will see, there automatically exists a flat direction containing the X field. 
The X field can have a large value as X ~ Mq during the inflation whereas it has a 
vacuum expectation value of (X) ~ F a in the true vacuum. Here, the Peccei-Quinn 
breaking scale F a is constrained by laboratory experiments, astrophysics, and cosmology 
as F a ~ 10 10 -10 12 GeV flll8|. 



We consider the following superpotential for the Peccei-Quinn symmetry breaking 
sector: 

W = XY(XX-F^), (3.54) 

where A is a coupling constant, and Y, X, and X are supermultiplets which are singlets 
under the standard-model gauge groups, which have 0, +1, and —1 Peccei-Quinn charges, 
respectively. From Eq. (|3.54f ) we see that there exists a flat direction XX = F%. We 
parameterize this direction by the scalar field a, which is called "saxion". This saxion 
a receives a soft SUSY breaking mass of m a ~ 777,3/2, and the X and X have vacuum 
expectation values of the order of F a ^ 

Along this flat direction the X field can have a large value during the inflation, if the 
additional SUSY breaking effects during the inflation induce a negative mass squared for 
X (see Sec. |3.2.2[). The initial value of X is then expected to be X ~ Mq, since the 



16 We assume here that the soft SUSY breaking masses for X and X are almost the same, i.e., mx — m-^ 
(~ m 3/2 ). 



98 



supergravity effects prevent the X field from running over M G [see Eq. Q2.57| )l, while the 
initial value of the X field is given by X ~ F^/ Ma- 
in order to fix the mass scale M e g at the true vacuum, let us take the mass of the 
lightest neutrino to be m v \ = 1 x eV for a representative value, with which the 
hierarchy of the neutrino masses becomes mild. Then the scale M e g is given by 

H 3xlO-Gev(^V<fV (3.55) 
\ m v i J \ F a J 

Here, we have used M e g oc (X) and M e g = (H u ) 2 jm v \ at the true vacuum (X) ~ F a . 
Although this scale M e g at the true vacuum is too low to produce the enough baryon 
asymmetry (see Fig. |3.2|), the large initial value of X, Xq ~ M G gives rise to 



— 25 -ram©' 

For such a high scale M e g, the evolution of the <p field is completely free from the thermal 
effects and the early oscillation can be avoided, even if the reheating temperature Tr is 
as high as T R ~ 10 8 GeV.0 

The baryon asymmetry is then obtained by using Eqs. ( fj.41 ) and ( |3.56j ) as 



where we have used H osc ~ ~ m 3 / 2 |am|- This shows that too large amount of the 
lepton asymmetry is produced. However, this asymmetry is sufficiently diluted, since 
there exists substantial entropy production by the decay of the saxion a. 

The saxion begins the coherent oscillation with the initial amplitude Xo ~ Mq at 
H ~ m CT ~ 7713/2, just after the production of the lepton asymmetry ends. The oscillation 
energy at that time is given by p a ~ When the reheating process of the inflation 

completes, the ratio of p a to the energy density p ra d of the radiation of the universe is 
estimated as 

— - 1 f^) 2 - (3-58) 

Prad 6 \M G ) 1 ^ 

Notice that p a decreases at the rate R~ 3 as the universe expands, while p ra d decreases as 
R~ A . Therefore, the oscillation energy of the saxion dominates the energy of the universe 
soon after the reheating process completes. Here, it is useful to take the ratio of p a to the 
entropy density s, 

-^-T R {^f]\ (3.59) 
s 8 H \M G J ' v ; 



17 



This can be seen by calculating the oscillation time H osc in Eq. (3.34) with M c g ~ 10 GeV 



99 



which takes a constant value until the saxion decays. 

This energy density of the saxion is transferred into the thermal bath when it decays. 
Here, it should be noted that the saxion might decay dominantly into two axions. If it is 
the case, the extra energy of the axion at the BBN epoch raises the Hubble expansion of 
the universe, which leads to overproduction of 4 He. To avoid this difficulty, the branching 
ratio of the saxion decay into two axion should be smaller than about 0.1. Here, we simply 
assume this is the case.[^] Then, the saxion decays dominantly into two gluons with the 
decay rate 

2 S 

Through this decay the universe is "reheated" again, and its reheating temperature T a is 
estimated as 

T ° = 10 GeV Irrev J {-^—)- (3 ' 61) 

Notice that the saxion decay takes place far before the beginning of the Big-Bang Nucle- 
osynthesis, and hence is cosmologically harmless. 

The saxion decay increases the entropy of the universe by the rate 

A^after Pa ^after 

"^before "^before Pa 



\M G J 3T ff 



. 6 , T R ^/lTeV^ 3/2 / F a \fX 



Because of this entropy production by the saxion decay, the primordial baryon asymmetry 
Eq. (|3.57 ) is also diluted by the rate A. Then the present baryon asymmetry is given by 



Notice that the present baryon asymmetry is independent of the reheating temperature 
T R , while it depends on the Peccei-Quinn scale F a . We see that the desired baryon 
asymmetry is obtained even with the lightest neutrino mass of m v \ ~ 10~ 4 eV for F a ~ 
10 10 GeV. 

Finally, we comment on the cosmological consequences of this model. 



18 This is realized when m \ ~ in—. 

x x 
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(i) The entropy production by the saxion in Eq. ( |3.62| ) ensures that we are free from 
the cosmological gravitino problems, since the number density of gravitinos produced at 
the reheating process is diluted by the rate A. 

(ii) This model does not suffer from the cosmological problem of the "axino", which is 
a fermionic superpartner of the axion. The axinos are produced in the reheating process 



by the thermal scatterings and may lead to a cosmological difficulty [ |119|| . However, in 
this model, the interaction of axino at the reheating epoch is suppressed by Xq ~ Mq, 
not by F a , and hence the production of axino is less effective. Furthermore, the entropy 
production by the saxion decay dilutes the axino abundance. Thus, the axino becomes 
completely cosmologically harmless in this model. 

(iii) The reheating temperature of the saxion decay Eq. ( |3.61|) is high enough to 
thermalize the lightest SUSY particles (LSPs) if the mass of the LSP is less than 200 GeV, 
and hence the stable LSP can be a dominant component of the dark matter. 
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3.4 Effects of the gauged U(l) b-l symmetry 



In Sec. P?2] , we found a very interesting aspect of the leptogenesis via LH U flat direction, 
the "reheating temperature independence of the cosmo logical baryon asymmetry." The 
resultant baryon asymmetry is almost determined by the mass of the lightest neutrino 
(and the gravitino mass), and hence we can predict the mass of the lightest neutrino m v \ 
as m v \ ~ 10~ 9 eV from the observed baryon asymmetry. 

In this section, we discuss the effects of gauging the XI (1) b-l symmetry on the lep- 
togenesis via the LH U flat direction. As mentioned in Sec. the absence of U(1)b-l 
gauge anomaly automatically requires three families of right-handed neutrinos, and the 
breaking of B — L symmetry naturally gives rise to the heavy Majorana masses of these 
right-handed neutrinos, which explains the tiny neutrino masses via the seesaw mecha- 
nism. 

Below the B — L breaking scale, the behavior of the scalar potential of the <fr held does 



not change much from the one discussed in Sec. |3.2|. Therefore, if the B — L breaking 



scale v is larger than the initial value of field obtained in Sec. 3.2.2, 



v>^/M eS H in{ , (3.64) 

the amplitude of the field is always smaller than the B — L breaking scale, and the 
dynamics of the field is basically the same as the one discussed in Sec. |3.2| . Hence, the 
resultant baryon asymmetry reduces to the one obtained in Sec. 3.2.4 [see Eq. ( |3.41| ) and 
Fig. |3.4||. We will call this case "F-term stopping case." 



On the other hand, if v ^5 \/M e gH- m f } the field develops its initial value as large as 
the B — L breaking scale. At this scale, the potential of the flat direction <ft can be lifted 
by the effect of the U(l) b-l -D-term ||106|| . (We denote this case by "D-term stopping 



case.") As we shall show in this section, the dynamics of the (f) field is drastically changed 
in this D-term stopping case. Actually, the obtained baryon asymmetry is much larger 
than the case without a X]{1)b-l symmetry, and hence than the F-term stopping case. 

It is interesting that the resultant baryon asymmetry in the D-term stopping case 
linearly depends on the reheating temperatures of inflation, but it is completely indepen- 
dent of the gravitino mass. This is totally an opposite situation to the result obtained 
in Sec. |3.2|, in which the final baryon asymmetry is almost independent of the reheating 



temperatures of inflation and linearly depends on the gravitino mass [see Eq. ( |3.41| ) and 



Fig. |3.4[|. This "gravitino- mass independence of the cosmological baryon asymmetry" 



provides us a great advantage in gauge-mediated SUSY breaking scenarios K5|. 
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3.4.1 The model 



In order to discuss the effect of the XJ(1)b-l gauge symmetry, we should investigate 
the B — L breaking sector, which gives the masses of the right-handed neutrinos. To 
demonstrate our point, we use the following superpotential: 



W 



hNLH u + ^SNN + XX (SS - v 2 ) + fiH u H d , (3.65) 



where, h, £, A are coupling constants and we assume A = 0(1). X, S and 5* are singlets 
under the MSSM gauge groups and they carry the U(l) b-l charges of 0, —2 and +2, 
respectively, v is the the breaking scale of the U(1)b_l symmetry. Here and hereafter, 
we take the couplings h, £, A and /i, v to be real, by field redefinitions. We will omit 
the family indices, since the "effectively flat" direction relevant for the present baryon 
asymmetry is again given by the flattest flat direction, L — L\. 

The shape of the superpotential in Eq. (|3.65|) is the same as that of the SUSY hybrid 



inflation discussed in Sec. [2.3.3. a| .M Actually, If the X field in Eq. ( ^.65| ) develops a large 
expectation value and the energy density of the universe is dominated by the vacuum 
energy V ~ \ 2 v , hybrid inflation takes place. Here, however, we assume that the B — L 
breaking sector in Eq. Q3.65p is not directly related to the inflation, and discuss the effects 
of this U(1)b-l on the dynamics of the LH U flat direction. 

As we have shown in Sec. |3.2.2 , the Kahler potential must have non-minimal couplings 



of the field to the inflaton x, since otherwise the <fi field gets a large positive mass term 
of the order of the Hubble parameter and it is driven exponentially towards the origin 
during inflation [|104||. We assume that there are also non- minimal couplings of other fields 



to the inflaton in the Kahler potential: 

K = ]Tyt r + x t x 



Y 



+ £(^^ + H.c.)+E^A^ + ..., (3.66) 

where cy and by are complex and real couplings of order unity, respectively, and Y denotes 
L, H u , Hd, X, S, S, and N. Then, the full scalar potential relevant to the flat direction 



19 



The superpotential in Eq. ( |3.65| ) is also similar to the one in the model in Sec. 3.3.2. However, the 



model adopted here is different from that, since the U(l) symmetry here is gauged, while in Sec. 13.3.2 
the Peccei-Quinn U(l) symmetry is global. 
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field L = H U = (j>/y/2 is given by Q 
1 



V 



-g 2 (-2|S| 2 + 2|S| 2 + |Af--|0| 2 



-H 2 + iSN 



+ h 2 \N<j)\' 



+ 



-£N 2 + XX S 
2^ 



+ X 2 



SS-v 2 \ +X 2 \XS\ 2 

+ 3# 2 ^T(l-& y )|y| 2 ^ -V3h(^2c y YW y + R.c^ +V sb , (3.67) 

where 6' y = by — \cy\ 2 and Wy = dW/dY, and we redefined the phase of cy. p = 0(1) is 
the gauge coupling constant of the U(l) b-l- V$b represents SUSY breaking terms in the 
true vacuum. All fields denote the scalar components of the corresponding superfields. 
Here, we omit the potential coming from thermal effects, which will be discussed in 
Sec. 



The scalar potential in Eq. ( |3.67| ) is so complicated that it seems very difficult to solve 
the dynamics of the relevant fields. Fortunately, however, we only need to know the shape 
around the minimum of this potential, if the curvatures around the bottom of the potential 
are as large as the Hubble parameter during inflation. This is because all the scalar fields, 
which have masses as large as the Hubble parameter, settle down at the bottom of the 
potential during inflation and trace this potential minimum in the subsequent evolution. 
Therefore, we first find out the relations to minimize the potential in Eq. ( |3.67| ). Though 
this is also a very hard task because of the complexity of the potential, one can find 
out the minimum of the potential at least when \<p\ and H^zv are satisfied. From 



Eq. ( |3.67| ), we find the approximate minimum of F-terms for 

•• 2 h<f) 2 



< 



S 



X ~ 0, N ~ 



(3.68) 



S' ' 2£S 

Notice that the curvatures associated with these three relations are of the order of the 
B — L breaking scale v, and hence these relations hold as far as H^zv. By using these 



relations, we can reduce the scalar potential in Eq. ( [3.671) to the following one: 



2 9 



H 2 



+ 2 



+ 



3(1 



,5| 2 4£ 2 |S| 2 
b' s )H 2 \S\ 2 + 3(1 



+ 



4£ 2 \S\ 



+ 



h 4 W 



64£ 2 |5| 4 
h 2 I, 



3(1-^^ 



20 Here, we take Hd — assuming a positive Hubble-order mass term for the Hd field. This is necessary 
in order to avoid a contamination of Hd to the LH U flat direction. We have checked both analytically 
and numerically that the contamination of the Hd field becomes relevant only after the H /i, and hence 
the following discussion is not affected much by the Hd contamination. 
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+ ^(c;|^ + H.c.)+^ 2 |0| 2 + ^sb, (3.69) 

where = — (l/4)c£. Here and hereafter, we assume 3(1 — b'^) ~ —1 for simplicity. 
(This corresponds to c# ~ 1 in the previous analysis. See Eq. (|3.13|) .) 

One may wonder why the "flat" direction field can develop a large expectation value 
in spite of the presence of the U{1)b-l D-tevm. This is because the S field shifts and 
absorbs the D-term potential: 



' - If7 + ;SW- (3.70) 



The curvature associated with this relation is also of the order of the B — L breaking 
scale. Then we get the following effective potential: 



16 



V ^ V SB + ^ 2 \<p\ 2 + ^H[c';-^ + R. C .)+^-^- + 



2^ m 2 V C 5 ' '/ 4( 2 |S| 2 VI^IV 
- # 2 |0| 2 + 3(1 - b' s )H 2 \S\ 2 + 3(1 - &'s)# 2 ^t + 3(1 - b' N )H° 



\s\ 2 v nj ^e\s\ 2 ' 

(3.71) 



In the following discussion, we assume H^v and that the relations in Eq. (|3.68| ) and 



in Eq. ( |3.70| ) are satisfied during the leptogenesis works. Although the approximation 
we have made here is crude, we will justify by a numerical calculation that the effective 
potential in Eq. ( |3.71| ) is valid at least for the following discussion. 



As we will see later, the third term in Eq. ( 3.71| ) 



V phase = ^H^j^ + E.c) , (3.72) 
which depends on the phase of 0, plays a crucial role in the D-term stopping case. We will 



call this term "Hubble A-term." (This corresponds to the Hubble A-term in Eq. ( 3.13 
in the previous analysis.) 



3.4.2 Initial amplitude and initial phase 

In this subsection, we discuss the evolution of the <fi field during the inflation. 
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3.4. 2. a initial amplitude 



First, we discuss the amplitude of the field. Let us first investigate the potential in the 
range \(j)\ ^ v. In this range, the S field in Eq. (|3.70|) is expanded into the following form: 

\S\ 2 * v 2 -\\^\ 2 + o{^Y for \<j>\<v. (3.73) 



By substituting Eq. (|3.73|) into Eq. (|3.71| ), we get the effective potential for 0^ v in the 
following form: 



v K ^ + ^^ + |JL (c y + H. c .) + 4A6 

H 2 \4>\ 2 + I (b' s - b' s ) H 2 \4>? + ■ ■ ■ , (3.74) 



where the ellipsis denotes higher order terms in \<p\ 2 /v 2 . Here, M e g = ^v/h 2 corresponds 
to the mass scale which was denoted by the same symbol in Sec. [3.2| . (The mass of the 
right-handed neutrino in the present model is given by Recall that it is related to 
the neutrino mass as m u i = (H u ) 2 /M eS . 

The scalar potential in Eq. ( |3.74| ) is nothing but the potential obtained in Eq. ( |3.12| ) 
in the previous analysis, with some modifications^ However, there is an important 
difference in the Hubble-induced SUSY breaking mass terms. Since there appears an extra 
Hubble mass term, we need the following condition in order that the field develops a 
large expectation value during the inflation: 

b's-b's<l. (3-75) 

If this is not satisfied, the field is driven toward the origin during the inflation, and the 
leptogenesis cannot work. In the following discussion, we assume that the above condition 
is satisfied. Then the field has a negative Hubble mass term at least in the range of 
\<P\<v. 

Now let us discuss the scale where the flat direction is lifted. If the balance point be- 
tween the negative Hubble mass term and the F-term potential |0| 6 / (4M e 2 ff ) in Eq. ( |3.74| ) 
is below the B — L breaking scale v , the field is stopped at this balance point, 



|0| ~ sJM^H ini <v. (3.76) 

This is the condition for the "F-term stopping case." In this case, the evolution of the 
field and the amount of the generated baryon asymmetry result in the same conclusions 
as discussed in Sec. WZ- 
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As we mentioned, the potential terms coming from thermal effects will be discussed in Sec. 3.4.4 
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On the other hand, if the condition 



M cS H hli >v, (3.77) 

is satisfied, the <p field can develop its expectation value as large as v. In this case, we 
can stop the <fi field below the balance point y/M^TT^i by using the D-term potential. 
This is very crucial to enhance the final baryon asymmetry in the .D-term stopping case. 
To see this, let us discuss what happens if the field develops its value as large as the 
B — L breaking scale v. At this scale (i.e., \<f>\ ~ v), the expansion of the S field given in 



Eq. ( 3.73 ) becomes invalid and above this scale we must use another expansion of the S 



field as follows: 

1,2 '12 , Ln ( "\ , r. ( h2 W | |2 for |^| >B (3 78) 



2 



One might wonder whether the above expansion is reliable, since Eq. ( 3.70Q is based on 



Eq. (pi), which may not be valid for \<ft\ 3> v. Here, we first derive the conditions to 
fix the field at the scale v during the inflation, assuming this expansion is effectively 
valid at least for \<f>\ ~ v. We will justify later the validity of the obtained conditions by 
numerical calculations. 



By substituting the expansion Eq. ( 3.78|) into Eq. ( 3.71|) , we obtain the following 
Hubble-induced mass termQ 

5V^ (-1 + ^(1 -b'^H^ + ol^H 2 ) for \<j>\>v. (3.79) 



4 

Therefore, the <fi field can have a positive Hubble mass term for \<j>\ £jV, if 

Therefore, if this condition is satisfied, the field cannot develop its expectation value 
above the scale v, and hence it is fixed at the B — L breaking scale v during the inflation. 

To summarize, the <fi field, and hence S and S fields as well, are stopped at the B — L 
breaking scale v during the inflation, if the following conditions are satisfied, 

1 



b' s -b' § <-, ^M cS H iQ{ >v, b' § <--. (3.81) 



These are the conditions for the D-term stopping case. 



2,2 Here, we assume h 2 /£ 2 -C 1. If this is not the case, we must include the Hubble mass term coming 



from the coupling of the right-handed Majorana neutrino to the inflaton, the last term in Eq. (3.71). 
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b's - b' s = 8/3 




— b' s = -1/3 




★ </>/<;< 0.3 




■ 0.3 < 0/« < 1 




a l < <f>/v < 3 




A 3 < <j>/v < 10 




10 < (j>/v 



b' 



Figure 3.6: The amplitudes of the field fixed during inflation which are determined 
by numerical calculations. In this calculation, we have used the full scalar potential in 
Eq. (|3.67 ) to follow the evolution of the relevant fields. Here, we have assumed that 
Hint/v = 0.1, 3(1 — 5^) = —1, h = 10~ 4 , and g = £ = A = 1, and have randomly generated 
other coupling constants by-, \cy\ in the range —2.5 < by < 2.5 and < \cy\ < 2.5, 
respectively. Various symbols denote the \<f>\ to v ratio at the end of inflation. We see 
that the <fi field is, in fact, stopped at the B — L breaking scale v if the conditions in 
Eq. ( |3.8ip are satisfied. 
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In order to check the reliability of the conditions in Eq. ( |3.81| ), (especially that of the 

/_ 
s 



last condition b's ^ —1/3,) we have numerically solved the coupled equations of motions 



dV 

Y + 3HY+ — = 0, (3.82) 
for the relevant fields Y = 0, N, S, S and X, using the full scalar potential in Eq. fl3.67| ). 



We show the result in Fig. |3Jd|, where the amplitude of the field at the end of the inflation 
is plotted in b' s -b'^ plane. Here, we have taken M e fj = £v/h 2 = 10 8 v and H m f = O.lv, and 
hence the second condition y/M eS H in{ ^ v is satisfied. 

It is found from Fig. |3.6| that the conditions in Eq. ( |3.81| ) well explain the result of this 



numerical calculation. Actually, for b' s — b'g ^ 8/3 (positive Hubble mass term for |0 
the amplitude is damped as |0| <C v, while for b'g ^ —1/3 (negative Hubble mass term for 
\4>\^v), the field overshoots the B — L breaking scale v. 

Thus, the field can stop at the B — L breaking scale v if the conditions b' s — b'g ^ 8/3 



and 6^<-l/3 are satisfied. In fact, we see from Fig. |3T0 that the amplitude of the (p field 
at the end of the inflation lies in the range 

|0| ~ (1-3) x v, (3.83) 

in most of the parameter space where the conditions in Eq. (|3.81| ) are satisfied. 



3.4. 2. b initial phase 

Now let us turn to discuss the phase of the <fi field at the end of the inflation, which is 
crucial to estimate the resultant lepton asymmetry. 

First, in the F-term stopping case, the amplitude of the field is fixed at A/M eff if inf . 
An important point here is that the curvature around the valley of the Hubble A-term 
potential in Eq. ( [3.72j ) (which we denote by the symbol Wp haso ) is of the order of the 
Hubble parameter in this case. This can be seen from the following relation: 

o 

phase 

- 9%| 2 ^ f . (3.84) 

Therefore, as discussed in Sec. [T2], the phase of the field settles down at the bottom 
of the valley of the A-term potential during the inflation, and hence the Hubble A-term 
cannot supply a torque to rotate the field. In this case, the relevant torque for the 
field only comes from the ordinary A-term potential proportional to the gravitino mass. 



a ,n 'phase 
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As stressed in Sec. |3.2.4|, this leads to a suppression of the baryon asymmetry in gauge- 



mediation models with a small gravitino mass. 

However, if the conditions in Eq. (|3.81|) are satisfied, namely, in the D-term stopping 
case, the curvature of the potential along the phase direction is smaller than the Hubble 
parameter during the inflation: 

d 2 



^"phase 



phase 

~ ^L\^^Hl ( \£ X <H? { . (3.85) 

In this case, therefore, there is no reason to expect that the field sits down at the bottom 
of the valley of the A-term potential in Eq. ( |3.72|) during the inflation. In other words, 
unless there is an accidental fine tuning of the initial phase of the field, it is generally 
displaced from the bottom of the valley of this A-term potential when the inflation ends. 

Therefore, the A-term potential kicks the field along the phase direction when the 
Hubble parameter becomes comparable to the curvature along the phase direction ^p hase . 
As we will show in the next subsection, this phase rotation caused by the Hubble A-term 
is very crucial to enhance the final baryon asymmetry. 



3.4.3 Lepton asymmetry 

Now let us discuss the dynamics of the field and calculate the resultant lepton asymmetry 
in the present model. In the F-term stopping case, the amplitude of the field is always 
smaller than the B — L breaking scale v [see Eq. (|3.76|) 1 and the evolution of the field 
is determined by the effective potential in Eq. (|3.74 ) plus thermal effects, which we have 



investigated in Sec. pT2| . Thus, the resultant baryon asymmetry is given by Eq. ( |3.41| ) (or 
Fig. 



We turn to the D-term stopping case. As discussed in the previous subsection, the 
amplitude of the field at the end of the inflation is given by 

|0| ~v, (3.86) 

and the phase of the field is generally displaced from the bottom of the valley of the 
Hubble A-term potential. 

After the end of the inflation, the field value of the field takes an almost constant 
value for a while. Then the Hubble A-term in Eq. ( |3.72| ) kicks the field along the phase 
direction when the Hubble parameter H becomes comparable with the curvature: 

H 

M~eS 



H 2 ~ f^phasc ( — 77 ) ■ (3-87) 
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Thus, the Hubble parameter at this time if p h is given by 

H^^i^r- ( 3 - 88 ) 

At this time, the phase of the field, arg(0), starts to oscillate around the bottom of the 
valley of the Hubble A-term potential. At the same time, the amplitude of the field, 
|0|, starts to decrease according to the following equation: 



M eS H, (3.89) 



eff) 



which is the balance point between the negative Hubble mass term and the |0| 6 /(4M f 
potential. 

Since the phase of the field starts its oscillation at H = H p ^, it has already had a 
large acceleration along the phase direction before it starts coherent oscillation along the 
radius direction around the origin when H = H osc . This H osc is determined by thermal 
effects, soft SUSY breaking mass term and /x-term in the same way as in Sec. [3.2| . Here, 
we have assumed that H ph > H osc , which we will discuss in Sec. |3.4.4j . 



In Fig. |3.7| and Fig. |3.8| , we show the evolution of the amplitude and the phase of the 
field obtained by numerical calculations, respectively. In these figures, t is the cosmic 
time in the matter dominated universe t = 2/(3H). We see from these figures that the 
evolution of the field is well explained by above arguments. 

We are now at the point to estimate the baryon asymmetry. From the equation of 



motion (|3.21|), the evolution of the lepton asymmetry til = (z/2)(0*0 — 0*0) is given by 



h L + 3Hn L = 2V3-^-Im (cU 4 ) , (3.90) 



M, 



eff 



where we have used \S\ ~ v for H ^ H ph . This can be easily integrated and we obtain 
the lepton number at time t as 



R 3 n L ](t) ~ |p f'dt'R'HlmU'^) . (3.91) 



Before the 0's oscillation, i.e., for H osc ^ H ^ H ph , the amplitude of the field decreases 
as |0| ~ ^/McffH and the sign of Im(c00 4 ) in Eq. ( |3.91| ) changes with 1/H time scale. (See 
Fig. |3]7] and Fig. |3.8| .) Therefore, we see from Eq. (|3.91|) that the total lepton number 



oscillates with 1/H time scale with almost constant amplitude, since R 3 H |0 4 | oc H oc t 1 
in this period. 

Then, the total lepton number is fixed when the field starts its coherent oscillation 
around the origin at H = H osc , since the amplitude of the field decreases as fast as 
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field estimated by a numerical calculation. We assume 
1(T 10 eV), v = 10 14 GeV, H M = 10 12 GeV, 



mvi 



Figure 3.7: The evolution of the 
that M cS = 3 x 10 23 GeV (i.e., 

= 10 3 GeV. We see that the field starts to move at H = H ph ~ 3 x 10 4 GeV 
according to the equation; \<p\ ~ y/M^H. Here, t = 2/(3H) is the cosmic time. The 
field starts to oscillate around the origin at H osc ~ since we neglect the thermal effects 
in this calculation. 



argQ) 
vr/4 




Figure 3.8: The evolution of the phase of the <fi field estimated by a numerical calculation. 
The parameters used in this figure are the same as in Fig. ft.7| . Here, we also defined that 
is real, and hence the valleys coming from the Hubble A-term lie along arg(0) = 7r/4, 
37t/4, 57r/4, 77r/4. We see that the phase of the field starts to oscillate around the 
bottom of the valley at H ~ i? p h- 



112 



|0| oc H and the integrand of Eq. ( |3.91| ) decreases as ~ t 3 after the oscillation. The 
lepton number density at time t osc = 2/(3H osc ) is then given by 

n L (t osc ) ~ ^^-H osc lm [cU 4 (t OS( 

~ -^=\c^\M cS H 2 J ph , (3.92) 

where 5 p h — sin(4arg0 + argc'^) is the effective CP-phase. By using Eq. ( |3.39| ), the final 
lepton asymmetry is given by 

nr Mpfl-Tp , , , „ . 
- 'c'J5 ph - (3-93) 



Here, we have used the fact that the (f) field starts its oscillation before the reheating 
process completes, i.e., H osc > T x , which was verified in Sec. |3.2.4| . Notice that the resul- 



tant lepton asymmetry is independent of the time when the field starts the oscillation 
around the origin, if osc .f^] This is because the total lepton number itself oscillates with a 
constant amplitude for H osc ^ H ^ H ph , as discussed just above. Thus, the order of the 
lepton asymmetry is already fixed just after H ~ i? p h, before the oscillation time H osc . 

By using the relation between the lepton and baryon asymmetry in the presence of 
the sphaleron effect, we obtain 

— ~ 0.35 x — 

Here, we give several comments on the obtained baryon asymmetry. First, by com- 
paring Eq. (|3.93|) and Eq. (|3.40|) , we can see that the lepton asymmetry in the D-term 
stopping case is enhanced compared with the F-term stopping case (and the case without 
^(I-)b-l), since it is no longer suppressed by the oscillation time H osc . 

Next, as stressed before, one sees that the resultant baryon asymmetry is independent 
of the gravitino mass m 3 / 2 . Therefore, we can obtain an enough baryon asymmetry even 
with a very light gravitino. This is the crucial point in this D-term stopping case. Here, 
we have to take care of the gravitino problem. For example, if the mass of the lightest 
neutrino is as small as m u i ~ 10 -10 eV, the present baryon asymmetry is explained with 
a low reheating temperature Tr ~ 10 4 GeV [see Eq. ( |3.94j )l. In such a low reheating 



temperature, we are free from the overproduction of gravitinos even in gauge-mediated 
SUSY breaking scenarios with the small gravitino mass m 3 / 2 = (9(100 keV) pUf . 



23 We have also confirmed by numerical calculations that the final lepton asymmetry / s is independent 
of the oscillation time H osc . 
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Figure 3.9: The evolution of the baryon asymmetry %/s estimated by a numerical cal- 
culation. Parameters used in this figure are the same as in Fig. [3]7| and Fig. [3.8| . Here, 
we also take Tr = 10 4 GeV. The thick line denotes the analytic estimation in Eq. fl3.94p . 
We can see that the asymmetry begins to oscillate with constant amplitude at H ~ H p ^ 
and it is fixed when the <fi field starts the oscillations around the origin. The analytical 
estimation of the baryon asymmetry agrees well with the obtained value by the numerical 
calculation. 



On the other hand, if the mass of the gravitino is large enough, say, m 3 / 2 ~ 1 TeV, we 
can avoid the cosmological gravitino problem even if the reheating temperature is rather 
high as T R ~ 10 9 GeV (see Sec. |1.5|) . In such a case, we see from Eq. ( 3.94j) that the mass 
of the lightest neutrino m v \ ~ eV is small enough to explain the present baryon 
asymmetry. 

Finally, we show the evolution of the baryon asymmetry %/s obtained by a numerical 
calculation in Fig. |3.9| . From this figure we see that the asymmetry begins to oscillate with 
almost constant amplitude at H = and it is fixed at the time when the field starts 
to oscillate around the origin, H = H osc . We also see that the amount of generated baryon 
asymmetry agrees well with our analytic estimation in Eq. ( |3.94| ). We confirm from this 
figure that the resultant asymmetry is well explained by the arguments described in this 
subsection. 
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3.4.4 Conditions to avoid an early oscillation 



So far, we have neglected the thermal effects and assumed that there is no early oscillation, 
i.e., ifph > H osc . If the <fi field starts its oscillation before H = H p ^, i.e., if H osc < H p ^, the 
resultant baryon asymmetry is strongly suppressed. (In this case, the resultant baryon 
asymmetry is inversely proportional to H% BC .) Thus, we verify the assumption if p h > H osc 
here. 

First, in order to avoid the early oscillation caused by the soft SUSY breaking mass 
and the /x-term, we need 

,,2 



P h 



M, 



(3.95) 



wr 



Next, the thermal mass term Cfc/|T 2 discussed in Sec. |3.2.3| may cause the early oscil- 
lation. This can be avoided if the following condition is satisfied: 



H 2 > E c k ftT 2 

fkW<T 



for H > K 



ph 



(3.96) 



Finally, we need the following third condition to avoid the early oscillation caused by 
the thermal log term T 4 ln(|0| 2 ), which was also explained in Sec. |3.2.3 : 

2 mA 

a g a s T 



H 2 > 



for H > H, 



P h 



(3.97) 



Here, we have used 



v for H > H, 



P h- 



From the second ( |3.9£| ) and third ( |3.97| ) conditions, the reheating temperature is 
bounded from above as 

„,2 



Tr < min 



1/2 „ ,1/2 , ,1/2 

a g ' a s MJ M c ^ j 



where 



T, 



max 



fiV 



3/2 



1/4; 



f'J 2 



1/2 . ,3/2 
G IV1 eS ■ 



(3.98) 



(3.99) 



We show the allowed regions which are free from the early oscillation in Fig. |3.10| , 
taking account the conditions in Eqs. ( |3.95|) and (|3.98|) . From this figure, we see that the 
early oscillations can be easily avoided, especially when the B — L breaking scale satisfies 
10 14 GeV. 



24 At first sight, Eq. (3.98) seems to be obtained just by rewriting the relation H osc < H p h by using 

the formula of H osc in Eq. (3.34). However, it is slightly different, since \<f>\ takes a constant value v for 

> 



H ?t -ffph in the present case, while \cj>\ ~ y/M e gH in the previous case. 
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Figure 3.10: The plot of the parameter regions which are free from the early oscillation in 
m v \ — Tr plane. The regions below the four solid lines are free from the early oscillation. 
These four lines correspond to the breaking scale of the U(1)#-l, v = 10 16 , 10 15 , 10 14 , 
10 13 GeV from left to right, respectively. The shaded region corresponds to the present 
baryon asymmetry, %/s — (0.4 — 1) x 10 -10 . We have taken £ p h = 1 and |c'J = 1 in 
this figure. We see that if we take the breaking scale of the U(1)b-l as v ^ 10 14 GeV, 
the early oscillation of the field H osc > can be avoided in most region of parameter 
space. 
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3.5 The neutrino- less double-beta decay 



In Sec. [3.2| , we have investigated the leptogenesis via LH U flat direction in detail and 
have found that the produced baryon asymmetry is determined almost independently of 
the reheating temperature Tr, because of the early oscillation of the flat direction field 
caused by thermal effects. As shown there, an ultralight neutrino m v \ ~ 10~ 9 eV is 
required to explain the empirical baryon asymmetry ub/s ~ 10~ 10 , except for the case 
where the model in Sec. [3.3.2| is adopted.^ By introducing a gauged \]{1)b-l symmetry, 
the early oscillation of the flat direction due to the thermal effect can be avoided (Sec. |3.4|) , 
but small neutrino mass m v \ <C 10~ 4 eV is still necessary if we require low enough 
reheating temperature Tr ^ 10 9 GeV to avoid the cosmological gravitino problem in a 
wide range of gravitino mass (see Fig. p.lOj ). In this section we show that such a tiny 



mass of the lightest neutrino mass m v \ <C 10~ 4 eV have an interesting implication on a 
low-energy experiment; the neutrino-less double beta (OvfiP) decay, which will be a strong 
evidence for the lepton number violation if observed. 

The crucial parameter to determine the Qvf3/3 decay rate is the v e -v e component of 
the neutrino mass matrix, which is given by 



(3.100) 



where the mixing matrix U a i is defined in Sec. |2.1| . If the mass of the lightest neutrino is 
actually so small as m v \ <C 10~ 4 eV, the contribution from m v \ to m ee can be neglected as 
long as we consider an accessible 0is(3(3 decay rate, which corresponds to \m ee \ ^ 10~ 3 eV. 
Thus, the parameter \m ee \ is written in terms of masses and mixings of two other neutrinos 



as 



\ m ee\ = \U 2 2 m u2 + U 2 3 m u3 \ , (3.101) 
which immediately leads to the following upper and lower bounds: 

|m ee |max = \U e 2\ 2 m u 2 + \U e3 \ 2 m u3 , 

U e2 \ 2 m u2 - \U e3 \ 2 m u3 . (3.102) 



m 



ee |min 



In general, the mass pattern of the neutrinos can be classified into two cases: the normal 
hierarchy and the inverted hierarchy. Let us consider these two cases in turn. 



25 See also the last comment in Sec. [3.2.4 
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normal mass hierarchy 



In the case of normal mass hierarchy, the parameters of atmospheric and solar neutrino 
oscillations are given by 



5m 



atm 



m„ 3 - m 



2 

v2 i 



5m 2 



sol 



m 



v2 5 



tan 2 6> so i = \U e2 /U e i\ 



(3.103) 



and the mixing angle constrained by CHOOZ experiment [|116|1 corresponds to sin^chooz = 
\U e s\, which is severely bounded as sin^hooz < 0.15. In terms of these oscillation param- 
eters, Eqs. ( |3.102|) become 



m. 



ee max 



Cchooz^soiV Sm soi + s chooz JSml tm + 5m 2 ol 



C chooz S sol \ I s 



2 

chooz 



Sml tm + Sm 2 



sol 



(3.104) 



cos 2 #x and s 2 x = sin 2 ^x- Here, we can safely neglect the contribution of 



where c 2 x 

Sm 2 ol in the second term s 2 



rjl ,, oz y 5m 2 tm + 5m 2 ol , since it changes |m ee | at most A|m ee | < 



chooz 



5m s 2 ol /(2 v /5m 2 tm ) < 10~ 4 eV. 



Then, the ratio of the mass parameter \m ee \ to the mass scale of atmospheric neutrino 



oscillation J 5m 2 tm is given by 



\m P 



^atm 



2 2 
^chooz sol 



\ 



S ™L , 2 
i m 2 6 chooz 5 



2 2 
^chooz^sol 



\ 



^»sol 
^ m atm 



5 chooz 



(3.105) 



Fig. RTTTJ shows the upper and lower bounds on \m ee \ in units of 



<5m 2 tm , for sin^hooz 



0.15, 0.10, 0.05 and 0. As can be seen in Fig. p.llj , the mass parameter \m ee \ is predicted 
in a narrow range. We stress that this is because of the absence of the contribution from 
the lightest neutrino, i.e., m v \ <C 10~ 4 eV. It is quite interesting to see that the predicted 
\m ee \ can be as large as \m ee \ ~ 0.001-0.01 eV, which could be accessible at future Ov/3/3 
decay experiments | 120 |. Furthermore, it is also found from Fig. |3.11| that, if the sin# c h 00 z 
becomes more constrained by future experiments, \m ee \ is predicted in a much narrower 
range. 
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Figure 3.11: The predicted value of \m ee \ in units of 
hierarchy. The lines represent the upper and lower values of |m ee 



^ m atm> i n the case of normal mass 



for sin6> cho o Z = 0.15, 



0.10, 0.05 and 0. The vertical line corresponds to the set of best fit values, sin 9 BO \ — 0.3, 
5m, 



sol 



7 x 10~ 5 eV 2 |IT3| and 5ml tm = 2.5 x 10~ 3 eV 2 [|10 
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Figure 3.12: The predicted value of \m ee \ in units of 
mass hierarchy. The lines represent the upper and lower values of \m 
corresponds to the best fit value, sin 2 6> so i 



<5m 2 tm , in the case of inverted 
The vertical line 
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inverted hierarchy 



Now let us turn to discuss the case of inverted mass hierarchy, where the parameters of 
atmospheric and solar neutrino oscillations are given by0 



5m 



atm 



m 



u3 > 



m u3 -m 



2 

v2 i 



tan 2 0. 



sol 



\u e3 /u e2 \ 2 



(3.106) 



The mixing angle constrained by CHOOZ experiment ||116| now corresponds to sin # c hooz 
\U e i\. Then from Eqs. ( |3.102| ) we obtain 



[in 



ee max 



TTlee min 



5m 2 



chooz 



r 2 

^chooz 



"sol 



\ 



1 - 



^ m atm 



atm 



"sol 



\ 



5m sol 
^ m atm 



'sol 



'sol 



(3.107) 



$mltm, for given 



In Fig. |3.12| , we show the upper and lower bounds on \m ee \ in units of 
tan 2 sol . Here, we have used sin^chooz = for simplicity, but \m ee \ changes (reduces) 
only 2% even for sm0 c h OO z = 0.15. As for the ratio r = <5m 2 ol /<5m 2 tm , we have adopted 
r = 0.008-0.15 [[10], |115|| . As seen from the figure, the \m ee \ is restricted in a very narrow 
range such as \m ee \ ~ (0.01-0.05) eV x (i75m 2 tm /0.05 eV) for tan 2 # sol < 0.7, which is 



indeed in the reach of the future 0v(3(3 decay experiments ||120|. 



The prediction of \m ee \ discussed in this section is a generic consequence of a hierar- 
chical neutrino masses m v \ <C mu2, m u3- Thus the Qv(3(3 decay experiments, combined 
with neutrino oscillation experiments, can provide a consistency test for the tiny mass of 
the lightest neutrino m v \ <C 10~ 4 eV required in the leptogenesis via LH U flat direction. 



26 The convention adopted here is different from the one in Ref. 1 103 , where ra„ 
m^i < m V 2 for inverted hierarchy. 



are defined as m v ^ < 
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3. A Z}-flat condition for the LH U flat direction 

In this appendix, we discuss the D-flat condition along the LH U direction. Let us take 
the number of generations of the lepton doublets Lj to be N in general. Then, there are 
(2 + 2N) complex scalar fields in the L-H u sector. 

u = ( 50 i \ 




H u = „o ■ ( 3 - 108 ) 



First, the freedom to make SU(2)^ gauge transformation allows us to rotate away a 
possible vacuum expectation value of one of the weak isospin components of one of the 
scalar fields. Thus, we can take H£ = without loss of generality. Then, vanishing 
.D-term potential requires the following relations: 

-\v£i\ 2 -\eii\ 2 + \H^\ 2 = forU(l)y, (3.109) 

VhftU + e£*VLi = 
vii^Li - zilvLi = for SU(2) L , 
\fii\ 2 - \e£i\ 2 - \H%\ 2 = (3.110) 

where the summations over % are suppressed for simplicity. From these equations, the 
D-flat condition is given by 

N 

£I^I 2 = I#T- (3.1H) 

It is still possible to rotate away one of the phases by using the U(l)y gauge transfor- 
mation (or third component of the SU(2)x gauge transformations^]). Therefore, this flat 
"direction" is parameterized by [(iV+1) complex fields] — [one real constraint Eq. (|3.11 I 



— [one phase] = 2N real fields. In the absence of effective superpotential Eq. ( |3.2j ), this 
"direction" would be completely flat in the SUSY limit. Furthermore, it would be isotropic 
for iV directions, v£i, ■ ■ • , vln- However, due to the scalar potential induced by the op- 
erator Eq. (p.2|), the potential is lifted and modified. Then, the flattest direction among 
the above field space is parameterized by one scalar field 0, which corresponds to the first 
family i = 1. 



27 



Notice that one U(l) group remains unbroken by the condensation of the flat direction field. 
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Chapter 4 



Conclusions and discussion 



In this thesis, we have investigated various leptogenesis mechanisms in the framework of 
SUSY. We have also discussed some aspects of those scenarios, such as constraints from 
cosmological gravitino problems and possible experimental signatures. In particular, we 
have made essential improvements compared with earlier works in the following points: 

• In the leptogenesis scenario where the lepton asymmetry is created by the decay 
of right-handed sneutrino N having dominated the early universe (Sec. |2.4j) , the 
dilution of the gravitino due to the entropy production by the iV's decay has been 
discussed. As a result, it has been shown that the cosmological gravitino problem 
is drastically ameliorated. 

• In the leptogenesis via the LH U flat direction (Chapter^), all the relevant thermal 
effects has been taken into account, and detailed analyses on the dynamics of the 
flat direction field has been done both by analytic and numerical calculations. 

In this chapter, we summarize the aspects of various leptogenesis mechanisms in- 
vestigated in this thesis. The produced baryon asymmetry in each mechanism is well 



summarized in Table 4.1 



First of all, we see that in all the mechanisms the neutrino mass plays a very crucial 
role to determine the resultant baryon asymmetry. In the mechanisms where lepton asym- 
metry is produced by the asymmetric decay of the right-handed (s)neutrino (Chapter |2p, 
the asymmetry parameter e\ is proportional to the mass of the heaviest neutrino m v z for 
fixed effective CP- violating phase 8 e g. We see that the information from the atmospheric 



neutrino oscillation [IU|, which leads to m v s — 0.05 eV as long as masses of the light 
neutrinos are not degenerate, is very important to calculate the final baryon asymmetry. 
On the other hand, when the lepton asymmetry is produced via the LH U flat direction 
(Chapter |3|), the mass of the lightest neutrino m v \ is the crucial parameter. 
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produced baryon asymmetry 
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(See Fig. ^ 
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s 



0.35 x — ;= _|c' |5 ph 



v^M 2 



G 



D-term stopping 



(5) 



1 x 10 



-10 



R 



10 8 GeV 



10" 6 eV~ 



-i 



Table 4.1: Summary of produced baryon asymmetries in various leptogenesis mechanisms. 
Necessary conditions are: (1) T R > M 1 . (2) M x > T R and m x > 2M 1 . (3) Mi > T Nl 
and T R > (for \N U \ ~ M G ). (4) negative Hubble mass for 0. (5) See Eq. fljQI 



and Fig. |3.10| . Some comments are in order: (a) Perturbative decay is assumed. For 
preheating, see Sec. [2.3.6| . (b) H in { > H osc is assumed. (See the last comment in Sec. |3.2.4| .) 
(c) This is not the case if effective mass for the right-handed neutrinos are very large during 



leptogenesis (Sec. |3.3.2|) 
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It is interesting to see that the dependences on the neutrino masses are opposite in 
those two scenarios. As for the t\ parameter of the right-handed neutrino decay, larger 
Yukawa coupling is favored to make e\ larger. When we normalize the Yukawa couplings 
by neutrino masses (Sec. |2TT| ), ei turns out to be proportional to m vi for fixed <5 e ff. On the 
other hand, for leptogenesis via LH U flat direction, flatter direction is favored to enhance 
the final lepton asymmetry. This is because, the flatter the potential is, the larger the 
amplitude of <fi becomes, and the larger the final lepton asymmetry is produced. Thus in 
this case the relevant neutrino mass is the lightest one, m v \. 

Now let us comment on each mechanism in turn. 



thermal production of the right-handed (s) neutrino 



The simplest one is the thermal production of the right-handed (s)neutrino (Sec. [2.2|) . The 
attractive point in this scenario is nothing but the fact that the production of the right- 
handed (s)neutrino requires no assumption besides high enough reheating temperature 
Tr > Mi. The production occurs just by ordinary thermal scatterings. 

On the other hand, since this scenario requires high enough reheating temperature, 
cosmological gravitino problems should be considered seriously. One solution is the stable 
(LSP) gravitino with mass m 3 / 2 ~ 100 GeV. In this case, one has to take care of the decays 
of the next-to-lightest SUSY particle (NLSP) after the big-bang nucleosynthesis f48| , 
which gives severe constraints on the lifetime and the abundance of the NLSP at the decay. 
Since both of the lifetime and the abundance of the NLSP is determined by the properties 
of the NLSP (and the gravitino mass), they might be going to be well calculated if SUSY 
particles are discovered and their properties are studied in future collider experiments. 
This is quite an interesting possibility. 

If the gravitino is not the LSP, cosmological gravitino problem is relatively severe. 
Therefore, it will be interesting and might be necessary to determine a definite lower bound 
on the reheating temperature Tr at which the leptogenesis from thermally produced right- 
handed (s)neutrino works marginally. This requires to see what happens if the reheating 
temperature is just around the mass of the right-handed neutrino, T R ~ Mi. 



• leptogenesis in inflaton decay 

The advantage of this mechanism is that the production of the right-handed (s) neutrino is 
possible for Tr < Mi. Thus the cosmological gravitino problem can be avoided in a wider 
range of the gravitino mass. Actually, we have shown that the enough lepton asymmetry 
can be produced even with a reheating temperature Tr ~ 10 6 GeV. On the other hand, 
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the produced baryon asymmetry crucially depends on the physics of the inflation. In 
particular, we have taken the branching ratio B r (x —>■ NiNx) to be B r ~ 1. In this sense 



the first one of the hybrid inflation models we have discussed in Sec. |2.3.3.a| , which has 
Yukawa couplings between the inflaton and right-handed neutrinos, is attractive since it 
automatically gives rise to B r ~ 1. 

We should mention that the current and future satellite experiments on anisotropies of 



the cosmic microwave background radiation [90] will give precious information about the 



inflation physics, especially the spectrum index n s . Thus we will be able to distinguish 
which inflation model is favored by using future data. 



By the way, from Table [4.1| , it is found that the formula of the baryon asymmetry from 
inflaton decay reproduces that of the iV-dominant universe if we take B r — 1, Tr — 
and 2Mi = m x . This corresponds to the case where the right-handed sneutrino plays a 
role of the inflaton itself [Ell, |91 



• iV-dominant universe 

The resultant baryon asymmetry in this scenario is expressed in a very simple way, as 



given in Table |4.1| . It is determined only by the decay temperature for fixed m v z and 
<5 e ff- This is because the iVi's decay produces dominant component of the radiation as 
well as the lepton asymmetry at once. 

The most attractive feature of this scenario is the dilution of the gravitino. In fact, as 
we have shown, cosmological gravitino problems can be avoided even for Tr ^> 10 11 GeV, 
which is very special among other baryogenesis scenarios. 

It is also quite an interesting possibility that the cosmic density perturbation coming 
from the iVi's fluctuation might be detected at future experiments @ . 



• LH U flat direction 

We have performed a detailed analysis on this scenario, including all the relevant thermal 
effects. In spite of the complicated potential and dynamics of the flat direction field 0, 
it was shown that the result is very simple and attractive. The baryon asymmetry is 
determined almost only by the mass of the lightest neutrino m u \. Remarkably, we can 
estimate (or predict) the m v \ from the empirical baryon asymmetry ns/s ~ (0.4-1) x 10 -10 
as 

m ul ~ (0.1 -3) x 10~ 9 eV, (4.1) 
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almost independently of the reheating temperature for Tr ^ 10 5 GeV. It is amazing that 
the baryogenesis which happened in the very early universe is directly related to the 
neutrino mass. 



• LH U flat direction + gauged \](1) B _ L , D-term stopping case 

When we introduce an gauged XJ(1)b-l symmetry the dynamics of the LH U flat direc- 
tion field cf) can be modified, if certain conditions are satisfied (D-term stopping case). 
Then the resultant baryon asymmetry is enhanced compared with the case with F-term 
stopping case, or the case without a gauged U(1)b_l symmetry. Although the "reheating 
temperature independence" disappears in this case (and m v \ is no longer predicted), this 
scenario has an advantage in the case of gauge-mediated SUSY breaking (GMSB) mod- 
els [45||, compared to the previous case. The point is that the phase rotational motion 



of the field is provided not by the usual SUSY-breaking A-term which is proportional 
to the gravitino mass, but by a Hubble-induced A-term. Thus this mechanism can work 
even with a very light gravitino as in the GMSB models. 



• The neutrinoless double beta decay 

We have also shown that a very light mass of the lightest neutrino m v \ <C 10~ 4 eV, which 
is suggested by the above leptogenesis mechanisms via LH U flat direction, leads to a high 
predictability of the rate of the neutrino-less double beta (Qv(3f3) decay. Although it is 
highly difficult to determine the mass of the lightest neutrino directly, future experiments 
on 0^/3/3 decay and neutrino oscillations will provide us precious information by which 
we can test the consistency of the very small lightest neutrino mass m v \ <C eV, as 
shown in Fig. |3.11| and |3.12| . 
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Appendix A 
Notations 



We use natural units c = h/(2ir) = 1, where c is the speed of light (in vacuum) and h is 
the Planck constant. We also take the Boltzmann constant ^ to be unity, which leads 
to 1 eV = 11605.fr. 

M G = (87rG) -1 / 2 = 2.4 x 10 18 GeV denotes the reduced Planck scale, where G is the 
Newton constant. 

Throughout this thesis, except for generic discussions in Sec. |3T], |3]^, |3]5| and |3.A| , we 
assume the existence of the right-handed neutrinos, which have superpotentials: 

W = ^MiNiNi + h ia NiL a H u , (A.l) 

where iVj (z = 1,2,3), L a (a = e, //, r) and H u denote the supermultiplets of the heavy 
right-handed neutrinos, lepton doublets and the Higgs doublet which couples to up-type 
quarks, respectively. 

We often use a basis where the mass matrix of the light neutrinos becomes diagonal, 
where the Yukawa couplings are given by hy. (The left subscript i denotes the family of 
the heavy right-handed neutrinos iVj, while the right subscript j denotes the family of the 
light neutrinos Vj in mass eigenstates. See Sec. [D] and Eq. fl2.12p .) 

As for the parameters related to the cosmology, see next Appendix [B[ 
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Appendix B 

Some notes of standard cosmology 



In this appendix, we briefly review some ingredients of the standard cosmology, in order 
to introduce our notations and relevant equations. For more details, see, for example, 
Ref. @. 



B.l Friedmann- Robertson- Walker universe 

The standard cosmology is based on the following homogeneous and isotropic metric 
(Robertson- Walker metric): 

where R(t) is the "scale factor" of the expanding universe, and t is the cosmic time. 
(t, r, 8, <fr) is called the comoving coordinates. From the above metric and the perfect- 
fluid form of the energy-momentum tensor T^ u = diag(p,p, p, p) (p(t) and p(t) denote the 
energy density and the pressure, respectively), the Einstein equations lead to the following 
two equations: 

[R\ 2 k 8ttG ,_ n , 

UJ + ^ = -T'' (R2) 

where the overdot denotes the derivative with cosmic time t, and G is the Newton constant. 



The second equation (|B.3|) , which corresponds to the conservation of the energy and 
momentum T^ v . v = 0, describes the dilution of the energy density due to the expansion: 

radiation (p rad = ip rad ) p rad oc R~ A , 
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matter (p mat = 0) p mat oc R 3 , 

vacuum (p vac = -p vac ) p vac = const . (B.4) 

Notice that the above relations hold for each component of the energy density px, as long 
as the energy-momentum conservation T x v . v = is satisfied separately for each X. 

The other equation (|B.2[), called Friedmann equation, describes the expansion of the 



universe. Notice that the curvature term k/R 2 decreases with time more slowly than 
the right-hand side of Eq. ( |B.2| ) as long as p is dominated by matter (p mat ) or radiation 
(Prad) [see Eq. (|B.4j )] . Therefore, it can be safely neglected as far as we consider the early 
universe. 

The Hubble parameter is defined as 

ff = |, (B.5) 

which means how fast the universe is expanding at a given stage of its expansion. In terms 
of this Hubble parameter H, reduced Planck mass M G = (87rG) -1 / 2 = 2.4 x 10 18 GeV, 
critical energy density p crit = 3MqH 2 and density parameter fl = p/p cr it, Eq. ( P-2|) is 
rewritten as follows: 

p i k 

Neglecting the curvature term, usually we use the following form of the Friedmann equa- 
tion: 



p = 3M^H 2 . (B.7) 



The present universe (t = to) observed is as follows [p9 j. The present Hubble expansion 
rate is H = 100 h [km sec -1 Mpc -1 ], where h ~ 0.7, and hence critical energy density is 
Pcrit(to) = 3MqHq = 1.9 h 2 x 10~ 29 [g cm -3 ]. In terms of density parameter fl x = Px/ Pcrit, 
total energy density is given by Q(to) ~ 1, which includes radiation fi ra d(to) = (2.47 + 
0.56^^) h~ 2 x 10~ 5 (N„ is the number of the generations of neutrinos which are still 
relativistic today), matter f2 ma t(to) — 0.3 (including baryon Qb^o) — 0.02/i~ 2 ) and "dark 
energy" fl vac (to) — 0.7.0 For more details and the errors included in these values, see 



Ref. [29]. The nature of the cold dark matter (dominant component of the fi ma t) as well 
as that of the "dark energy" are big puzzles in cosmology, which is beyond the scope of 
this appendix. (As for the cold dark matter, however, supersymmetry provides a good 



candidate, the lightest SUSY particle (LSP). See also Sec. |2.4.3j . 



1 The notation "il V ac" might not be correct, since the "dark energy" density could be a dynamical 
value (like "quintessence"), not a cosmological constant. 
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From the above ratio of the present energy densities, it is found that the early universe 
was dominated by radiation (for 10~ 4 -R(t )), and subsequently dominated by matter. 
Eqs. (|B.4|) and (|B.7|) lead to the following familiar relations for each epoch.[] 

radiation— dominated universe R oc t 1//2 , H = — , 

2t 
2 

matter— dominated universe R oc t 2//3 , H = — . (B.8) 

Ob 

B.2 Thermodynamics 

Here, we briefly review some basic thermodynamics in the radiation-dominated early 
universe, in which many particles are relativistic and in thermal equilibrium. First of all, 
the equilibrium density of particles of type i with momenta in a range d 3 p centered on p 
is given by 

9i^fiiP) , (B.9) 

where Qi is the number of degrees of freedom and /i(p) is the Fermi-Dirac or Bose-Einstein 
distribution function: 

/■(P)= fEl i . ■ (B.10) 
exp (^jA) ± 1 



Here, E$ is the energy = y p 2 + mf, fii is the chemical potential of the particle i, and 
the plus (minus) sign is for fermions (bosons). The number density Hi, energy density pi 
and pressure Pi of particle i are then given by the following equations: 



* = ^lMP)d 3 P, (B.ll! 



2w 



Eifi(p)d*p, (B.12) 



* = &!m m **- (R13) 



In Table [B.l| , we show these quantities for the relativistic (T ^> m^) and non-relativistic 
(T <^ mi) limits. Here, we have assumed |/^| < T and no Bose-Einstein condensation 
< rrii). 



2 If we include the change of g*(T) with time, which is defined below, the relation in Eq. ( |B.S| ) is 
slightly modified. 
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T > mi 




fermion 


boson 


n *=n 2 J exp l tJ 

Pi = mjUj 
^ = Tnj (< ^) 


pi = h (io) T4 

Pi = ~pi 


Hl ~ 9i ( vr 2 ) T ' 

= ft (io) T4 

1 

Pi = oPi 



Table B.l: The number density rii, energy density p, and pressure pi of the particle i, 
which is thermal equilibrium, in the limits of T ^> rrii and T <^ rrii. We have assumed 
|pi| < T and |pi| < m,. 



Because the energy density of a non-relativistic particle is exponentially suppressed 
compared with the relativistic one, the total energy density of the radiation p rad is given 
by the following simple form: 

Prad = |^*Cn^ 4 , (B.14) 

where 

9*(T)= E 9i + I E 9i- (B.15) 
rrii < T rrij < T 

i = boson j = fermion 

If there are particles which have different temperatures from that of the photon T, another 
factor (Tj/T) 4 should be multiplied in the above expression. (For example, at T MeV, 
neutrinos have temperature T v = (4/ll) 1//3 T for m u <C T u .) 

Notice that all leptogenesis mechanisms discussed in this thesis work at temperatures 
far above the electroweak scale T ^> 1 TeV, where all the MSSM particles are expected 
to be in thermal equilibrium. In this case, we obtain 

#, = 228.75 for MSSM. (B.16) 



In the expanding universe, it is convenient to introduce the entropy density s, which 
is defined by 

p + p 
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(Again, in the presence of particles with different temperatures, a factor of (Tj/T) 3 is 
multiplied in Eq. ( |B.17|) . In this case, the in Eq. ( B.17 ) becomes slightly different from 
the g* in Eq. ( |B.14j) .) Notice that the entropy per co moving volume sR 3 is conserved as far 
as no entropy production takes place. Thus it is quite convenient to take the ratio rix/s 
when we discuss some number density nx- For example, if some AT-number is conserved, 
the ratio of the A-number density to the entropy density takes a constant value 



s 



const , 



:b.i8) 



as long as there is no entropy production, since both nx and s scales as R~ 3 as the 
universe expands. As another example, if the A-particle is in thermal equilibrium and 
relativistic (T 3> m x ), the ratio is given by 



eq 
n X 



45C(3) g x 



x - for fermion 
4 



(B.19) 



s 2tt 4 g m (T) 

where the temperature (or time) dependence only comes from g*(T). 

Before closing this section, we calculate the relations between the particle number 
asymmetry n[ — n\ ^ and the particle's chemical potential /ij, which can be obtained by 
integrating Eq. ([B.ll ). In order to calculate the asymmetry, it is necessary to calculate 
higher order terms than those in Table |B.1|. The results are given by 



n 



(+) 



n 



(-) 



n 



(+) 



n 



(-) 



V 3 
o 

-g{T 3 
3 y 



for fermion . 
for boson , 



(B.20) 



where ellipses denote higher order terms in the expansions of rrii/T and fii/T. Here, we 
have assumed no Bose- Einstein condensation |/^| < for boson, and relativistic limit 
mi T. Eq. ( |B.20| ) is used in Sec. O, in deriving the relation between the baryon 
and lepton asymmetries in the presence of sphaleron effect. In this case, the chemical 
potentials are of order ~ 10~ 10 T, and hence the above assumptions are justified. 
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